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number  of  wheels  is  not  limited.  Two  programs  were  developed,  one  called 
WESLIQID  and  the  other  WESLAYER.  The  former  is  for  pavements  on  liquid 
foundations  and  the  latter  is  for  linear  layered  elastic  solids.  Variable 
slab  thicknesses  and  moduli  of  subgrade  reaction  are  incorporated  in  the 
WESLIQID  program,  and  any  number  of  slabs  arrayed  in  an  arbitrary  pattern 
can  be  handled.  Because  of  larger  computer  storage  requirement  and  compu¬ 
tational  complexity,  the  WESLAYER  program  is  limited  to  two  slabs^ 

The  theoretical  background  of  the  model  is  presented  first^  sind  the 
capability  and  logic  of  the  programs  are  fully  described.  A  discussion  of 
the  load  transfer  mechanism  along  the  joints  and  cracks  is  then  given. 
Results  computed  by  the  WESLIQID  program  were  compared  with  those  of  avail¬ 
able  solutions,  such  as  the  Westergaard  solution  and  Pickett  and  Ray's  in¬ 
fluence  charts,  and  the  comparisons  were  very  favorable.  The  comparisons 
with  the  discrete-element  computer  program  were  excellent,  except  that  the 
edge  stresses  computed  by  the  discrete  element  program  were  much  smaller. 
The  WESLIQID  and  WESLAYER  programs  were  used  to  analyze  test  pavements  at 
eight  U.  S.  Air  Force  bases  and  at  the  Ohio  River  Division  Laboratories. 

The  comparisons  were  centered  on  the  percent  stress  transfer  across  the 
joint,  and  the  comparisons  were  good.  Based  on  the  conclusions  of  the  com¬ 
puted  results,  the  design  implication  of  WESLIQID  in  the  rigid  pavements 
is  discussed. 
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PREFACE 


The  study  described  herein  was  sponsored  by  the  Office,  Chief  of 
Engineers,  U.  S.  Army  (OCE) ,  as  a  part  of  the  Mobility  and  Weapons  Ef¬ 
fects  Technology  RDT&E  Project  No.  4A762719AT40 ,  Work  Unit  001,  "Air¬ 
field  Pavement  Design  and  Parametric  Sensitivity  Analysis,"  and  Work 
Unit  003,  "Rigid  Airfield  Pavement  Load-Deformation  Response  Analysis." 
The  study  was  conducted  during  the  period  from  1  October  1978  to 
30  September  1979- 

This  report  is  Report  1  of  a  three-report  series  concerning  the 
computer  programs  WESLIQID  and  WESLAYER,  which  provide  for  analysis  of 
rigid  multi component  pavements  with  discontinuities  on  liquid  founda¬ 
tions  (WESLIQID)  and  on  linear  layered  elastic  solids  (WESLAYER).  This 
report  presents  the  theoretical  background  and  numerical  results  and 
discusses  the  capability  of  the  two  programs  and  their  logic.  Reports  2 
and  3  are  user's  manuals  for  WESLIQID  and  WESLAYER,  respectively. 

The  study  was  conducted  by  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  Geotechnical  Laboratory  (GL),  under  the  gen¬ 
eral  supervision  of  Dr.  Don  C.  Banks,  Acting  Chief,  GL;  Dr.  Paul  F. 
Hadala,  Assistant  Chief,  GL;  and  Mr.  Alfred  H.  Joseph,  Chief,  Pavement 
Systems  Division  (PSD),  GL.  Dr.  Yu  T.  Chou,  PSD,  was  in  charge  of  the 
study  and  is  the  author  of  the  report.  Professor  Y.  H.  Huang  of  the 
University  of  Kentucky,  who  originally  developed  the  computer  programs, 
assisted  in  the  study. 

COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were  Com¬ 
manders  and  Directors  of  the  WES  during  this  study  and  the  preparation 
of  this  report.  Mr.  Fred  R.  Brown  was  Technical  Director. 
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APPENDIX  A:  EQUATIONS  FOR  STRESSES  AND  DEFLECTIONS  UNDER 


A  POINT  LOAD  AND  UNDER  A  CIRCULAR  LOAD . A1 

APPENDIX  B:  NOTATION  .  B1 


CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 

UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  he  con¬ 
verted  to  metric  (Si)  units  as  follows: 


Multiply 

By 

To  Obtain 

Fahrenheit  degrees 

0.555 

Celsius  degrees  or  Kelvins* 

feet 

0.30U8 

metres 

inches 

2.5U 

centimetres 

pounds 

( force) 

1*.  1+1*8222 

newtons 

pounds 

inch 

(force)  per 

cubic 

0.2714 

megapascals  per  metre 

pounds 

(force)  per 

inch 

175.1268 

newtons  per  metre 

pounds 

inch 

(force)  per 

square 

6.89475? 

kilopascals 

*  To  obtain  Celsius  (c)  temperature  readings  from  Fahrenheit  (F) 
readings,  use  the  following  formula:  C  =  0.555(F  -  32).  To  obtain 
Kelvin  (K)  readings,  use:  K  =  0.555(F  -  32)  -  273.15. 
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STRUCTURAL  ANALYSIS  COMPUTER  PROGRAMS  FOR  RIGID 


MULTICOMPONENT  PAVEMENT  STRUCTURES  WITH  DIS¬ 
CONTINUITIES—' WESLIQID  AND  WESLAYER 

PROGRAM  DEVELOPMENT  AND  NUMERICAL  PRESENTATIONS 

PART  I :  INTRODUCTION 


Background 

1.  The  determination  of  stresses  and  deflections  in  concrete  pave¬ 
ments  due  to  wheel  loads  has  been  a  subject  of  major  concern  for  nearly 
half  a  century.  In  the  early  1920's,  Westergaard  (1925)  assumed  the  sub¬ 
grade  to  be  a  Winkler  foundation*  and  assumed  that  the  slab  was  infinite 
in  extent  in  all  directions  away  from  the  load,  and  used  the  theory  of 
elasticity  to  develop  a  mathematical  method  for  determining  the  stresses 
in  concrete  pavements  resulting  from  corner,  edge,  and  interior  loads. 

In  extending  the  method  to  airport  pavements,  he  later  developed  new 
formulas  (Westergaard  1939 >  19^*8)  that  give  the  stresses  and  deflections 
at  an  edge  point  far  from  any  corner  and  at  an  interior  point  far  from 
any  edge.  These  formulas  were  then  employed  by  Pickett  and  Ray  (1951) 
to  develop  influence  charts,  which  have  been  used  by  the  Portland  Cement 
Association  (1955>  1966)  for  the  design  of  highway  and  airport  pavements. 

2.  In  spite  of  their  wide  acceptance  and  usage,  the  Westergaard 
solutions  have  been  subject  to  many  criticisms,  including  the  following: 

a.  The  solutions  are  based  on  an  infinitely  large  slab,  with 
a  load  at  the  corner,  on  the  edge,  or  in  the  interior. 

They  may  not  be  applicable  to  today's  airfield  pavements 
for  aircraft  equipped  with  large  multiple-wheel  gear  loads. 

b.  The  assumption  of  a  Winkler  foundation  is  not  realistic 
because  a  Winkler  foundation  consists  of  a  series  of 
springs  in  which  the  pressure  at  any  point  between  the 


*  A  Winkler  foundation  is  also  called  a  liquid  foundation.  The  inten¬ 
sity  of  the  reaction  of  the  subgrade  is  assumed  to  be  proportional  to 
the  deflection  of  the  slab  and  to  be  vertical  only;  frictional  forces 
are  neglected. 
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slab  and  the  subgrade  is  directly  proportional  to  the 
deflection  only  at  that  point  and  not  elsewhere. 

c_.  The  slab  and  the  subgrade  may  not  always  be  in  full  con¬ 
tact  as  assumed  in  the  Westergaard  solution.  Gaps  are 
frequently  observed  in  the  subgrade  near  the  joint  because 
of  pumping  action  or  plastic  deformation.  Temperature 
warping  can  also  cause  the  slab  to  curl  up  and  lose  con¬ 
tact  with  the  subgrade. 

d.  Westergaard  solutions  are  based  on  an  infinitely  large 

slab  with  no  discontinuities  and  thus  could  not  be  applied 
to  analyze  stress  conditions  at  a  joint  or  at  a  crack. 

3.  In  the  1960's,  a  discrete  element  method  based  on  the  finite 
difference  technique  was  developed  at  the  University  of  Texas  (Hudson 
and  Matlock  1966)  to  analyze  concrete  slabs.  The  method  considers  the 
slab  to  be  an  assemblage  of  elastic  joints,  rigid  bars,  and  torsional 
bars.  This  method  of  modeling  was  helpful  in  visualizing  the  problem 
and  forming  the  solution.  It  does  give  reasonable  values  for  pavement 
deflections,  but  there  are  problems  in  achieving  accurate  stress  values 
along  the  edges.  Serious  problems  exist  in  the  analysis  of  joints, 
cracks,  and  gaps  under  the  slab  because  of  the  nature  of  the  method. 

b.  The  Corps  of  Engineers  (CE)  realizes  that  much  of  the  main¬ 
tenance  of  rigid  pavements  is  associated  with  cracks  and  joints.  The 
current  CE  rigid  pavement  design  procedures  (Department  of  the  Army  and 
the  Air  Force  1970)  have  certain  limitations  that  were  imposed  by  the 
state  of  the  art  at  the  particular  stage  of  development.  During  the 
development  of  the  procedure,  it  was  necessary  to  make  simplifying  assump¬ 
tions  and  in  many  instances  to  ignore  the  effects  of  cracks  and  joints. 
Since  the  advent  of  high-speed  computers  and  the  development  of  the  finite 
element  method,  a  more  comprehensive  investigation  of  the  state  of 
stress  at  pavement  joints,  cracks,  and  other  locations  in  multicomponent 
pavement  structures  is  now  tractable  within  the  assumptions  of  the  theory 
of  elasticity.  Consequently,  a  better  and  more  reasonable  design  method 
may  be  developed  for  rigid  pavements. 

Purpose 

5.  The  purpose  of  the  study  was  to  develop  two-dimensional  workable 
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finite  element  computer  programs  that  have  the  capability  of  analyzing 
stress  conditions  in  a  rigid  pavement  containing  cracks  and  joints  and  in 
the  supporting  subgrade  soil.  The  programs  should  be  able  to  analyze 
slabs  made  up  of  two  layers  of  materials  with  different  engineering 
properties,  and  should  be  able  to  accommodate  full  or  partial  loss  of 
subgrade  support  over  designated  regions  of  the  slabs.  The  subgrade 
soil  can  be  either  the  Winkler  foundation  or  a  layered  elastic  solid. 

The  program  should  be  easy  and  economical  to  operate. 

Scope 


6.  The  finite  element  computer  programs  originally  developed  by 
Professor  Y.  H.  Huang  (Huang  and  Wang  1973,  197^+;  Huang  197^+a,  197^8) 
of  the  University  of  Kentucky  were  modified  and  extended  to  suit  the 
purpose  of  this  study.  The  programs  were  developed  based  on  the  two- 
dimensional  plate-bending  theory.  Two  computer  programs  were  developed: 
one  named  WESLIQID  and  the  other  WESLAYER.  WESLIQID  is  developed  for 
subgrade  soil  represented  as  a  Winkler  foundation.  The  program  can 
treat  any  number  of  slabs  connected  by  steel  bars  or  other  load  trans¬ 
fer  devices  at  the  joints.  WESLIQID  can  be  applied  to  two-layer  slabs, 
either  bonded  or  unbonded.  WESLAYER  is  for  subgrade  soil  represented 

as  either  a  linear  elastic  solid  or  a  linear  elastic  layered  system, 
because  of  additional  computer  storage  space  and  other  computational 
complexity,  WESLAYER  is  limited  to  two  slabs  connected  by  load  transfer 
devices . 

7.  Report  1  of  this  series  presents  the  basic  theoretical 
development  of  the  programs.  Explanations  are  given  in  the  concept  of 
stress  transfer  along  the  joint  and  the  capability  of  the  programs. 
Numerical  results  are  presented  comparing  the  values  computed  by  the 
computer  programs  with  field  measurements  and  with  those  computed  by 
the  Westergaard  and  other  available  programs.  The  design  implication 
of  the  computed  results  are  discussed. 

8.  Reports  2  and  3  are  user's  manuals  for  WESLIQID  and  WESLAYER, 
respectively.  Descriptions  of  the  two  programs  are  presented  in  detail, 


and  the  programming  approaches  and  logics  are  also  explained.  The  flow¬ 
charts  and  input  guidance  are  presented  with  several  example  problems 
that  illustrate  the  use  of  the  input  guides. 


PART  II:  FINITE  ELEMENT  PLATE-BENDING  MODEL 


Introduction 


9-  To  analyze  the  stress  conditions  of  a  rigid  pavement  involv¬ 
ing  cracks  and  joints,  the  most  ideal  representation  of  such  a  system 
would  be  the  application  of  a  three-dimensional  finite  element  method. 
The  inherent  flexibility  of  such  an  approach  permits  the  analysis  of  a 
rigid  pavement  with  steel  bars  and  stabilized  layers  and  provides  an 
efficient  tool  for  analyzing  stress  conditions  at  the  joint.  Unfor¬ 
tunately,  such  a  procedure  would  require  a  tremendously  large  amount  of 
computer  core  space  to  solve  an  extremely  large  number  of  simultaneous 
equations.  Such  a  procedure  is  highly  uneconomical  and  impractical 
until  much  larger,  less  costly  computers  become  available.  Also,  the 
use  of  a  three-dimensional  finite  element  method  still  could  not  solve 
some  basic  problems  existing  in  a  rigid  pavement,  such  as  the  loss  of 
support  between  the  pavement  and  the  subgrade  due  to  temperature  warping 
or  other  causes.  The  difficulty  lies  in  satisfying  the  continuity  con¬ 
ditions  in  +he  three-dimensional  finite  element  method.  However,  this 
problem  does  not  exist  in  the  two-dimensional  plate-bending  model  used 
in  this  study. 

10.  Recently,  several  researchers  used  the  finite  element  plate¬ 
bending  model  for  analysis  of  concrete  pavement  with  considerable 
success.  They  are:  Eberhardt  (1973a,  1973b),  Huang  and  Wang  (Huang  and 
Wang  1973,  197**;  Huang  197**a,  197**b),  Pichumani  (1971),  and  Tabatabaie 
and  Barenberg  (1978).  The  obvious  advantage  of  the  plate-bending  model 
is  that  it  is  two-dimensional;  it  can  thus  save  greatly  on  the  computer 
core  space  and  computing  time  and  can  make  the  model  workable  and  more 
acceptable  to  the  general  users. 

11.  After  a  thorough  review  of  the  available  models,  it  was 
decided  that  the  models  developed  by  Huang  and  Wang  (1973,  197*0  were 
more  complete  than  the  others.  Mainly,  these  programs  consider  the 
partial  subgrade  contact  of  the  pavement  and  elastic  subgrade 
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foundation,  which  are  essential  considerations  in  a  rational  pavement 
design. 


12.  The  finite  element  method  employed  in  this  study  is  "based  on 
the  classical  theory  of  a  thin  plate  by  assuming  that  the  plane  before 
bending  remains  a  plane  after  bending.  Because  the  slab  is  modeled  as 
a  thin  plate,  there  is  no  variation  in  vertical  deflection  along  the 
thickness  of  the  plate;  i.e.,  the  deflection  at  the  top  of  the  plate 

is  the  same  as  that  at  the  bottom.  When  the  slab  is  divided  into 
rectangular  finite  elements,  the  division  is  made  only  in  the  longitu¬ 
dinal  and  transverse  directions;  the  vertical  direction  is  not  needed. 
The  model  is  thus  two-dimensional.  Another  advantage  of  the  plate¬ 
bending  model  is  that  the  application  of  finite  element  method  does  not 
involve  the  subgrade  soil  and  thus  saves  computer  time.  Only  the  sub¬ 
grade  reactive  forces  acting  at  the  nodes  are  important.  The  subgrade 
reactive  forces  are  evaluated  by  numerical  procedures. 

13.  The  procedure  of  the  model  can  be  found  in  many  textbooks  and 
papers,  such  as  Zienkiewicz  and  Cheung  (1967)  and  Cheung  and  Zienkiewicz 
(1965),  and  will  not  be  presented  herein.  Only  the  general  approach  is 
described. 

Slabs  on  the  Winkler  foundation 

lh.  Figure  1  shows  a  rectangular  finite  element  with  nodes*  i  , 
j  ,  k  ,  and  St  .  At  each  node,  there  are  three  fictitious  forces  and 
three  corresponding  displacements.  The  three  forces  are  a  vertical 


force  F  ; 
w 

y-axis  M 

y 

z-direction 
the  y-axis 


a  moment  about  the  x-axis  M  ;  and  a  moment  about  the 

x 

The  three  displacements  are  the  deflection  in  the 

w  ;  a  rotation  about  the  x-axis  0  ;  and  a  rotation  about 

x 

9  .  These  forces  and  displacements,  for  plates  on  a 


Winkler  foundation,  are  related  by 


*  Symbols  used  in  this  report  are  listed  and  defined  in  the  Notation 
(Appendix  B). 
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where 

[K]  =  stiffness  matrix  of  the  slab,  the  coefficients  of  which 
depend  on  the  dimensions  a  and  b  of  the  element  and 
the  Young's  modulus  and  Poisson's  ratio  of  the  slab 

6  =  displacements  in  the  slab 

6  *  =  displacements  in  the  subgrade 

k  =  modulus  of  subgrade  reaction 

At  any  given  node  i 


where  w.  is  vertical  deflection. 


2a 


FORCES  AND  CORRESPONDING  DISPLACEMENTS 
Figure  1.  Rectangular  plate  element 
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15.  The  stiffness  matrix  for  a  rectangular  element  was  tabulated 
in  Table  7.1  of  the  book  by  Zienkiewicz  and  Cheung  (1967)  and  is  used 
in  the  analysis.  The  type  of  elements  used  is  isoparametric.  For  illus¬ 
trative  purposes,  assuming  i  ,  j  ,  k  ,  and  l  are  1,  2,  3,  and  k, 
respectively.  Equation  1  can  be  expanded  into  12  simultaneous  equations: 


where 

k^  =  vertical  force  (index  l)  Fw^  at  node  i  due  to  vertical 

deflection  (index  l)  wj  at  node  j 

k^2  =  vertical  force  (index  l)  Fwi  at  node  i  due  to  rotation 

about  x-axis  (index  2)  0xj  at  node  j 

k  =  moment  about  x-axis  (index  2)  MXi  at  node  i  due  to 

rotation  about  y-axis  (index  3)  0yj  at  node  j 

k  =  moment  about  y-axis  (index  3)  Myi  at  node  i  due  to 

J  rotation  about  x-axis  (index  2)  0xj  at  node  j 


It  is 


with  similar  meanings  for  k^  ,  k ^  »  ^21  ’  ^13  * 

also  to  be  noted  in  Equation  3  that  vertical  deflection  vi  at  the 
subgrade  surface  is  implicitly  assumed  to  be  equal  to  vertical  deflec¬ 
tion  w^  of  the  slab. 

16.  By  superimposing  stiffness  matrices  over  all  elements  and 
replacing  the  assumed  fictitious  nodal  forces  with  the  statistical 
equivalent  of  the  externally  applied  loads,  a  set  of  simultaneous  equa¬ 
tions  can  be  obtained  for  solving  the  unknown  nodal  displacements. 

IT.  For  illustrative  purposes.  Figure  2  shows  a  pavement  with 


Figure  2.  Computation  of  stresses  under  a  point 
load  acting  at  the  corner  of  a  slab 

8l  nodes.  Since  each  node  has  three  unknowns,  there  are  2U3  (8l  *  3 
=  21*3)  simultaneous  equations  to  be  solved  for  the  2b3  unknown  displace¬ 
ments.  The  nodal  numbering  system  shown  in  Figure  2  (9  nodes  in  the 
vertical  direction)  indicates  the  half  bandwidth  is  (9  +  2)  x  3  =  33  , 
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so  the  number  of  coefficients  in  the  stiffness  matrix  [K]  needing  to 
be  stored  in  the  core  is  only  33  x  2^3  =  8,019  ,  instead  of  a  full 
matrix  2h 3  x  2^3  =  59 ,0k9  elements.  The  reason  that  the  coefficients 
outside  the  half  bandwidth  need  not  be  stored  is  that  they  are  all  zeros. 

18.  Once  the  nodal  displacements  are  computed,  the  nodal  moments 
are  then  computed  using  the  stress  matrix  tabulated  in  Table  7.2  of 
Zienkiewicz  and  Cheung  (1967).  The  nodal  stresses  are  then  computed 
from  the  nodal  moments.  Because  the  stresses  at  a  given  node  computed 
by  means  of  one  element  might  be  different  from  that  computed  by  means 
of  neighboring  elements,  the  stresses  in  all  adjoining  elements  were 
computed  and  their  average  values  obtained. 

Slabs  on  an  elastic  foundation 

19-  Similar  to  Equation  1  for  the  Winkler  foundation,  the  rela¬ 
tionship  for  the  forces  and  displacements  can  be  written  as: 

{F}  =  ([K]  +  [H]){6}  (5) 

in  which 


and 

{F}  =  externally  applied  nodal  forces 

[K]  =  stiffness  matrix  of  the  slab,  the  coefficients  of  which 

depend  on  the  finite  element  configuration  and  the  flexural 
rigidity  of  the  slab 

[H]  =  stiffness  matrix  of  the  subgrade,  the  coefficients  of  which 
depend  on  the  nodal  spacings  and  the  Young's  modulus  and 
Poisson's  ratio  of  the  subgrade 

{6}  =  nodal  displacements,  each  consisting  of  a  vertical  deflec¬ 
tion  and  two  rotations 

20.  The  characteristics  of  Equation  5  are  such  that  the  stiffness 
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matrix  of  the  slab  [K]  is  banded,  but  the  stiffness  of  the  subgrade 
matrix  [H]  is  not  banded.  When  the  two  matrices  are  added,  the 
composite  matrix  [P]  ,  i.e.,  [p]  =  [K]  +  [H]  ,  is  not  banded.  An 

iterative  scheme  was  developed  by  Huang  (I97^b)  to  make  [P]  a  banded 
matrix  to  save  computer  storage  space. 

21.  To  illustrate  the  iterative  scheme,  consider  a  simple  example 
for  a  slab  divided  into  only  two  finite  elements  with  a  total  of  six 
nodes.  Because  each  node  has  three  unknown  displacements,  there  are 
18  simultaneous  equations: 


K11  K12 

K13 

Kll 

0 

0 

K22 

K23 

K2U 
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K33 

K3U 

K35 

K36 
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F5 

H66 
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L  6J 

where 

K. .  =  3  by  3  submatrix 
1 J 

-  3  by  3  subraatrix 
i  =  nodal  numbers  from  1  to  6 
j  =  nodal  numbers  from  1  to  6 


(6) 


and 
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where 


I 


8  .  =  rotation  about  x-axis 

xi 

0  .  =  rotation  about  y-axis 
yi 

F  ^  =  equivalent  vertical  force  determined  by  statics 
As  no  angular  continuity  is  assumed  between  the  slab  and  the  subgrade 


"h .  0  0 

ij 

0  0  0 

0  0  0 


(8 


where  h.,  is  vertical  force  at  node  i  due  to  vertical  displacement 
at  node  j  . 

22.  The  developed  scheme  is  to  transfer  part  of  the  f H ]  matrix 
to  the  right  side  of  Equation  6,  or 
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(9 


First,  assume  the  displacement  {6}  on  the  right  side  of  Equation  9 
as  zero  and  compute  a  set  of  nodal  displacements  {6}  by  solving 
Equation  9-  Enter  the  displacements  thus  computed  into  the  right  side 
of  Equation  9,  and  find  a  new  set  of  displacements.  The  process  is 
repeated  until  the  displacements  converge  to  a  specified  tolerance. 

23-  In  separating  the  [ H ]  matrix,  any  half  bandwidth  may  be 
used.  The  coefficients  within  the  half  bandwidth  are  placed  on  the 
left  of  Equation  9,  whereas  those  outside  the  half  bandwidth  are 
transferred  to  the  right.  It  was  found  that  when  the  number  of  equa¬ 
tions  is  large  while  the  half  bandwidth  is  small,  the  displacements  may 
not  converge,  and  a  larger  bandwidth  should  be  used. 

2h.  In  solving  Equation  9,  each  slab  can  be  considered  separately. 
First  consider  the  left  slab  only,  and  the  first  set  of  nodal  displace¬ 
ments  is  determined  from  Equation  8.  Next,  consider  the  right  slab. 

The  existence  of  the  left  slab  has  two  effects  on  the  right  slab: 

(a)  the  vertical  deflections  along  the  joint  in  the  right  slab  must  be 
set  equal  to  those  in  the  left  slab;  and  (b)  the  vertical  nodal  forces 
in  the  right  slab  must  include  those  due  to  the  deflections  of  the  left 
slab.  Because  the  deflections  of  the  left  slab  have  been  previously 
determined,  the  nodal  displacements  of  the  right  slab  can  be  computed. 
Then  return  to  the  left  slab  again.  The  existence  of  the  right  slab 
also  has  two  effects  on  the  left  slab:  (a)  the  displacements  of  the 
right  slab  will  induce  a  set  of  nodal  forces  along  the  joint,  which 
must  be  transferred  to  the  left  slab;  and  (b)  the  deflections  of  the 
right  slab  will  induce  vertical  reactive  forces  in  the  left  slab.  By 
taking  these  effects  into  consideration,  a  new  set  of  displacements  for 
the  left  slab  is  determined.  The  process  is  repeated  until  the  nodal 
displacements  converge. 

Basic  differences  between  the  liquid 
(Winkler)  foundation  and  the  elastic  foundation 

25.  The  basic  difference  between  the  liquid  foundation  and  the 
elastic  foundation  is  that  in  a  liquid  foundation,  the  deflection  at 
a  given  node  depends  only  on  the  forces  at  the  node  and  does  not  depend 
on  forces  or  deflections  at  any  other  nodes.  In  an  elastic  foundation. 
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however,  the  deflection  at  a  given  node  depends  not  only  on  the  forces 
at  the  node,  but  also  on  the  forces  or  deflections  at  other  nodes. 

26.  In  both  the  WESLIQID  and  WESLAYER  computer  programs,  separate 
nodal  numbers  are  used  at  nodes  at  either  side  of  the  joint.  The  dis¬ 
tance  between  the  two  nodes  is  zero;  i.e.,  these  two  nodes  are  physi¬ 
cally  one  node.  The  use  of  separate  nodal  numbers  is  necessary  because 
the  stresses  and  deflections  in  the  two  slabs  across  the  joint  may  not 
be  the  same.  The  presence  of  separate  nodal  points  across  a  joint  is 
not  a  problem  in  a  liquid  foundation  because  in  this  case  the  deflection 
at  a  nodal  point  along  the  joint  does  not  depend  on  deflections  else¬ 
where.  In  an  elastic  foundation,  however,  the  problem  is  more  involved. 
For  instance,  a  deflection  at  a  node  away  from  the  joint  will  produce  a 
certain  subgrade  reactive  force  at  nodes  along  the  joint. 

Description  and  Capability  of  the  Programs 


Program  logic 

27.  The  programming  approaches  for  programs  WESLIQID  and  WESLAYER 
are  presented  separately  in  Reports  2  and  3  of  this  series,  respectively. 
For  convenience  of  discussion,  the  basic  logic  of  the  programs  is  pre¬ 
sented.  Two  cycles  of  iterations  are  involved  in  the  programs.  One  is 
for  checking  the  subgrade  contact  condition,  and  the  other  is  for  check¬ 
ing  pavement  shear  forces  for  a  deflection  convergence.  At  the  outset 
of  the  computation,  the  program  first  assumes  a  full  contact  between  the 
slab  and  the  subgrade,  except  at  the  nodes  where  gaps  are  preassigned. 

The  gaps  may  be  a  result  of  pumping  or  plastic  deformation  of  the  sub¬ 
grade.  Deflections  are  computed  sequentially  for  each  slab  by  succes¬ 
sive  approximations  until  the  deflection  convergence  criterion  is  met. 

28.  The  relationships  between  the  slabs  along  the  joints  are 
such  that  after  the  deflections  are  computed  for  slab  i  ,  the  deflec¬ 
tions  are  superimposed  to  the  adjacent  (i  +  l)t}l  slab  through  the 
Joint,  and  when  the  iteration  returns  to  computing  the  deflections  of 
slab  i  ,  shear  forces  exist  at  slab  i  along  the  joint  between  slabs 
i  and  (i  +  l)  induced  by  the  deflections  of  slab  (i  +  1)  . 
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29.  When  the  deflections  have  all  converged,  the  deflection  at 
each  node  is  checked  as  to  whether  the  contact  condition  has  changed. 

If  the  condition  has  not  changed,  the  computed  deflections  are  the  final 
values;  if  it  has  changed,  a  new  subgrade  contact  condition  is  assigned 
and  the  iteration  for  deflection  computation  starts  again  until  the  sub- 
grade  contact  condition  ceases  to  change.  To  assure  proper  convergence, 
an  under relaxat ion  factor  is  used  in  the  programs  to  change  the  relaxa¬ 
tion  factor  automatically. 

Subgrade  types 

30.  As  was  explained  earlier,  the  application  of  the  two- 
dimensional  finite  element  process  does  not  involve  the  subgrade  soil. 

Only  the  subgrade  reactive  forces  between  the  subgrade  and  the  slab  at 
each  node  are  important.  The  subgrade  reactive  forces  can  be  due  either 
to  deflection  of  the  subgrade  at  the  node  or  to  deflections  of  an  adja¬ 
cent  slab  transferred  through  the  joint.  The  subgrade  reactive  forces 
at  nodal  points  are  combined  with  the  externally  applied  forces  when  the 
displacements  are  being  solved  in  the  simultaneous  equations  shown  in 
Equation  1.  The  subgrade  reactive  forces  can  be  evaluated  readily  in 
the  case  of  Winkler  foundation  but  are  more  laborious  in  the  case  of  the 
elastic  foundation.  These  forces  are  explained  in  the  following  paragraphs. 

31.  Winkler  foundation.  The  reactive  force  between  the  subgrade 
and  the  slab  at  each  node  equals  the  product  of  the  modulus  of  the  sub¬ 
grade  reaction  k  and  the  deflection  w  at  the  node.  For  reactive 
forces  at  nodes  along  the  joint  that  are  induced  by  the  deflections  of 
adjacent  slabs,  the  forces  are  computed  through  the  stiffness  matrix  of 
the  elements  adjacent  to  the  joint. 

32.  Elastic  foundation.  Boussinesq's  solution  and  Burmister's 
layered  elastic  solution  are  used  to  compute  subgrade  surface  deflec¬ 
tions  for  the  cases  of  a  homogeneous  elastic  foundation  and  a  layered 
elastic  foundation,  respectively.  Once  the  flexibility  matrix  is  formed, 
a  matrix  inversion  subroutine  is  used  to  invert  the  flexibility  matrix 

to  the  stiffness  matrix.  The  subgrade  stiffness  matrix  is  not  banded, 
and  at  each  node  there  is  only  a  vertical  component. 
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Stresses,  strains,  and  deflections 
in  the  supporting  subgrade  soil 

33.  Once  the  subgrade  reactive  forces  between  the  subgrade  and 
the  slab  at  each  node  are  determined,  stresses  and  strains  in  the  sup¬ 
porting  subgrade  soil  can  be  computed.  The  stresses  and  strains  are  in¬ 
duced  by  the  nodal  reactive  forces,  but  the  forces  are  acting  in  the 
direction  opposite  to  those  when  the  stress  conditions  in  the  slab  were 
computed.  When  the  subgrade  soil  is  represented  by  the  Winkler  founda¬ 
tion  (WESLIQID  program),  the  Boussinesq's  equations  can  be  used  to  com¬ 
pute  the  stresses  and  strains  induced  by  the  concentrated  nodal  forces. 

In  order  to  use  the  equations,  an  equivalent  elastic  modulus  E  corre¬ 
sponding  to  the  modulus  of  subgrade  reaction  k  (used  in  the  program  to 
compute  stress  condition  in  the  slab)  should  be  selected.  When  the  sub¬ 
grade  soil  is  represented  by  a  layered  elastic  foundation,  Burmister's 
layered  elastic  solution  is  used  for  computation.  Since  the  stresses 
and  strains  in  the  subgrade  soil  under  the  concrete  slabs  are  very  small, 
the  principle  of  superposition  is  valid  and  is  used  to  compute  the 
stresses  and  strains  in  the  soil  induced  by  all  the  nodal  forces.  It 
should  be  pointed  out  that  at  the  subgrade  surface,  the  deflection  at 
the  subgrade  at  a  node  is  the  same  as  that  of  the  concrete  slab  at  the 
same  node;  the  vertical  stress  at  the  node  is  equal  to  the  reactive 
force  acting  at  the  node  divided  by  the  affected  area. 

3^.  In  the  WESLIQID  computer  program,  the  Boussinesq's  equations 
for  a  point  load  (Harr  1966)  are  used  to  compute  the  stresses  mid  deflec¬ 
tion  in  the  supporting  elastic  subgrade  soil.  In  using  the  equations, 
however,  the  stresses  and  deflections  become  infinitely  large  or  inde¬ 
terminate  at  the  surface  directly  under  the  point  load  and  at  locations 
very  close  to  the  point  load.  However,  closed- form  solutions  for  uni¬ 
formly  applied  circular  loads  are  available  only  for  the  vertical  stress 
and  vertical  deflections  w  directly  under  the  center  of  the  cir¬ 
cular  load.  These  equations  are  presented  in  Appendix  A.  For  computed 
locations  in  the  subgrade  soil  not  directly  under  a  node,  the  computa¬ 
tions  are  made  based  on  the  point  loads  acting  at  the  nodes.  Since  the 
computed  values  at  locations  close  to  a  point  load  may  be  erroneous,  the 
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computations  are  not  made  at  shallow  depths  (less  than  1  in.*  in  the 

program),  except  for  the  vertical  stress  o  and  vertical  deflection 

z 

w  at  computed  locations  directly  under  a  node.  This  problem  does  not 
exist  in  the  layered  elastic  subgrade  soil  (WESLAYER  program)  in  which 
uniformly  applied  circular  loads  are  used. 

35*  When  the  subgrade  soil  is  represented  as  an  elastic  solid 
(WESLAYER  program) ,  the  nodal  reactive  forces  are  larger  at  the  slab 
edges  than  at  the  slab  center  where  the  load  is  applied.  The  distribu¬ 
tion  of  contact  pressures  under  a  rigid  footing  can  be  found  in  the 
textbook  by  Terzaghi  and  Peck  (1962). 

Subgrade  contact  options 

36.  The  concept  of  full  and  partial  subgrade  contact  was  origi¬ 
nally  developed  by  Huang  (1971).  Most  materials  presented  in  this  sec¬ 
tion  are  taken  from  this  source. 

37-  Discussion  on  the  Westergaard  solution.  Complete  subgrade 
contact  condition  was  assumed  in  the  Westergaard  solution.  The  slab 
always  has  full  contact  with  the  subgrade  soil ,  and  gaps  are  not  allowed 
between  the  slab  and  subgrade,  no  matter  how  much  the  slab  has  warped 
upward  due  to  temperature  change  or  to  the  applied  load.  In  other  words, 
the  slab  is  supported  by  a  group  of  springs,  and  the  springs  are  always 
connected  to  the  slab.  In  reality,  the  pavement  can  lose  subgrade  sup¬ 
port  at  some  parts  due  to  temperature  warping,  pumping,  and  plastic  de¬ 
formation  of  the  subgrade.  Results  from  the  Arlington  test  (Teller  and 
Sutherland  1935 »  1936,  19^2)  indicated  that  the  pavement  and  the  sub¬ 
grade  were  not  in  full  contact  even  when  the  slab  was  flat  and  there  was 
no  temperature  differential  between  the  top  and  the  bottom.  It  was  also 
found  that  the  stresses  in  concrete  pavements  due  to  corner  loading  de¬ 
pended  strongly  on  the  condition  of  warping.  When  the  corner  was  warped 
down  and  the  slab  and  subgrade  were  in  full  contact ,  the  observed  corner 
stresses  checked  favorably  with  Westergaard' s  solutions.  However,  the 
observed  stresses  were  **0  to  50  percent  greater  when  the  corner  was 
warped  up.  Consequently,  Westergaard' s  equation  for  corner  loading  was 

*  A  table  of  factors  for  converting  U.  S.  customary  to  metric  (Si)  units 
of  measurement  is  presented  on  page  b. 
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modified  by  Bradbury  (1938),  Kelley  (1989),  Spangler  (19^2),  and  Pickett 
(1951)  to  account  for  the  loss  of  subgrade  contact  due  to  temperature 
warping,  pumping,  and  plastic  deformation  of  the  subgrade.  These 
modifications  were  based  on  empirical  results,  and  no  theoretical 
methods  to  the  author's  knowledge  have  been  developed  so  far.  With  the 
advent  of  high-speed  computers  and  the  finite  element  method,  it  is  now 
possible  to  analyze  concrete  pavements  subjected  to  warping  and  loading 
without  assuming  that  the  slab  and  the  subgrade  are  in  full  contact. 

38.  Stres  ses  and  deflections  due  to  temperature  warping.  Fig¬ 
ure  3  shows,  in  an  exaggerated  scale,  a  thin  slab  subjected  to  a 
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Figure  3.  Warping  of  a  concrete  slab 

temperature  differential  AT  between  the  top  and  the  bottom.  If  the 
slab  is  weightless  and  unrestrained,  it  will  form  a  spherical  surface 

with  a  radius  R  .  Because  the  slab  is  only  slightly  curved,  the  length 

of  the  arc  on  the  upper  surface  is  practically  the  same  as  that  on  the 
lower  surface,  so  the  length  L  of  the  upper  surface  is  shown  as  the 
length  of  the  lower  surface.  The  length  is  actually  greater  at  the 
bottom  than  at  the  top  by  aLAT  ,  where  a  is  the  coefficient  of 

thermal  expansion.  Since  the  radius  R  is  much  greater  than  the  thick¬ 

ness  h  and  L  is  much  greater  than  aLAT  ,  it  can  be  easily  shown 
from  geometry  that 


and 


c 


(10b) 


where 

c  =  initial  curling  of  a  weightless  and  unrestrained  slab  due  to  a 
temperature  differential  between  the  top  and  the  bottom 

d  =  distance  to  the  center  of  slab  where  curling  is  zero 
Substituting  Equation  10a  into  Equation  10b  gives 
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Note  that  AT  is  positive  when  the  slab  is  warped  up  with  a  temperature 
at  the  top  smaller  than  that  at  the  bottom  and  negative  when  it  is 
warped  down. 

39-  The  general  formulation  involving  warping  is  similar  to  that 
for  loading.  After  the  stiffness  matrix  is  superimposed  over  all 
elements  and  the  nodal  forces  are  replaced  with  the  statistical  equiva¬ 
lent  of  the  externally  applied  loads,  the  following  simultaneous  equa¬ 
tions  can  be  obtained  for  solving  the  nodal  displacements: 


(K]{6)  =  {F}  +  k[A]{6'} 


(12a) 


where 

[A]  =  diagonal  matrix  representing  the  area  over  which  subgrade 
reaction  is  distributed 

{6'}  =  subgrade  displacements 

Note  that  the  second  term  on  the  right  side  of  Equation  12a  represents 
the  nodal  forces  due  to  the  subgrade  reaction.  If  the  slab  has  a  total 
of _  n  nodes,  then 
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where 

w  =  vertical  deflection,  downward  positive 

F  =  vertical  force  due  to  externally  applied  load,  downward 
positive 

Note  that  c  =  0  when  there  is  no  warping.  The  reason  that  F  and 
<5^  contain  only  one  nonzero  element  is  that  the  nodal  forces  are  deter¬ 
mined  by  statics  and  only  vertical  loads  and  reactions  are  involved. 

ho.  Analysis  based  on  full  contact.  The  concept  of  full  contact 
between  the  slab  and  subgrade  can  best  be  explained  by  the  spring  anal¬ 
ogy  shown  in  Figure  U.  Figure  ka  shows  a  foundation  consisting  of  a 
series  of  springs,  each  representing  a  nodal  point  in  the  finite-element 
analysis.  When  a  slab  is  placed  on  the  foundation,  the  weight  of  the 
slab  will  cause  a  precompression  of  the  springs,  as  shown  in  Figure  kb. 
Because  the  slab  is  uniform  in  thickness,  each  spring  will  deform  the 
same  amount,  and  no  stresses  will  be  induced  in  the  slab.  The  amount 
of  precompression  can  be  determined  directly  by  dividing  the  weight 
of  the  slab  per  unit  area  by  the  modulus  of  subgrade  reaction.  When 
the  temperature  is  colder  at  the  top  of  the  slab  than  at  the  bottom, 
as  is  usually  the  case  at  night,  part  of  the  slab  will  deflect  upward, 
as  shown  in  Figure  kc.  However,  the  slab  and  the  springs  still  remain 
in  contact  because  the  upward  deflections  are  smaller  than  the 
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Figure  1* .  Fpring  analogy  for  full  contact 


precompression.  The  deflection  of  the  slab  due  to  warping  can  be  deter¬ 
mined  by  subtracting  the  precompression  due  to  the  weight  of  slab  from 
the  deflection  due  to  the  weight  and  the  warping  combined,  as  indicated 
by  the  shaded  area  in  Figure  l*c.  The  result  is  exactly  the  same  as  when 
the  warping  alone  is  considered.  The  same  is  true  when  a  load  is  applied 
to  a  warped  slab,  as  shown  in  Figure  Ud.  Therefore,  when  the  slab  and 
the  subgrade  are  in  full  contact,  the  principle  of  superposition  applies. 
The  stresses  and  deflections  due  to  warping  and  loading  can  be  deter¬ 
mined  separately,  one  independent  of  the  other,  disregarding  the  weight 
of  the  slab.  This  principle  forms  the  basis  of  Westergaard' s  analysis. 

1*1.  The  assumption  that  the  slab  remains  in  contact  with  the  sub¬ 
grade  implies  that  the  subgrade  reaction  always  exists  no  matter  how  the 
slab  is  warped.  If  the  slab  is  warped  up,  the  subgrade  will  2 ull  the 
slab  down,  and  a  deflection  w  is  obtained,  as  shown  in  Figure  3.  The 
displacement  of  the  subgrade  is  thus  c  -  w  ,  as  indicated  by  Equa¬ 
tion  12c.  If  w  within  {6'}  in  the  second  term  on  the  right  side  of 
Equation  12a  is  moved  to  the  left  and  combined  with  w  on  the  left  and 
c  is  combined  with  { F }  ,  Equation  12a  becomes 

t K] { 6 }  =  {?}  (13) 
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where 


[K]  =  composite  stiffness  matrix  of  the  system 
(F}  =  composite  nodal  forces 

If  there  is  no  warping,  then  c  =  0  ,  and  {F}  =  {F}  .  Equation  13  is 
used  in  the  computer  programs  to  compute  the  displacements. 

1*2.  The  above  derivation  for  upward  warping  also  ap^xied  to  down¬ 
ward  warping.  When  the  slab  is  warped  down,  the  temperature  differential 
is  negative.  If  the  temperature  differential  is  the  same  for  downward 
warping  as  for  upward  warping,  the  stresses  and  deflections  will  be  the 
same  in  magnitude  but  opposite  in  sign. 

1*3.  Analysis  based  on  partial  contact.  The  major  difference  in  pro¬ 
cedure  between  full  and  partial  contact  is  that  it  is  not  necessary  to  con¬ 
sider  weight  of  the  slab  in  case  of  full  contact,  but  for  partial  contact 
weight  of  the  slab  must  be  considered.  The  latter  case  involves  two  steps. 
First,  gaps  and  precompressions  of  the  subgrade  due  to  weight  of  the  slab 
and  warping  combined  are  determined.  These  gaps  and  precompressions  are 
then  used  to  determine  stresses  and  deflections  due  to  applied  loads. 

1*1*.  It  should  be  noted  that  full  contact  is  a  special  case  of  par¬ 
tial  contact.  Every  problem  in  partial  contact  is  analyzed  first  by  as¬ 
suming  that  the  slab  and  subgrade  are  in  full  contact.  If  it  turns  out 
that  they  actually  are  in  full  contact,  no  iterations  are  needed.  If  some 
points  are  found  out  of  contact,  the  reactive  force  at  those  points  is  set 
to  zero.  The  process  is  repeated  until  the  same  contact  condition 
is  obtained. 

1*5-  Partial  contact  without  initial  gaps.  This  case  applied  to 
new  pavements  not  subjected  to  a  significant  amount  of  traffic.  Each 
spring  in  the  Winkler  foundation  is  in  good  condition  and  if  the  slab  is 
removed,  will  rebound  to  the  same  elevation  with  no  initial  gaps,  as 
shown  in  Figure  5a.  Under  the  weight  of  the  slab,  each  spring  is 
subjected  to  a  precompression,  as  shown  in  Figure  5b.  If  the  slab  is 
warped  up,  gaps  will  form  at  the  exterior  springs,  as  indicated  by  a 
positive  s  in  Figure  5c,  and  precompressions  will  form  at  the  interior 
springs,  as  indicated  by  a  negative  s  .  If  the  slab  is  warped  down, 
all  springs  will  be  under  precompression,  as  shown  in  Figure  5b,  except 
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(a)  Initial  subgrade  condition 
without  initial  gap 


(b)  Precompression  due  to 
weight  of  slab 


LOAD 


)  dap  and  preeompress ion 
ue  to  weight  and  warping 


(d)  Deflection  due  to 
loading  (shaded  area) 


Figure  5*  Partial  contact  without  initial  gaps 
that  the  precompressions  are  not  equal.  The  displacements  due  to  the 
weight  of  the  slab  and  warping  combined  can  be  determined  from 
Equation  12a,  except  that  the  subgrade  displacements  are  expressed  as 


and 
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when  w .  >  c . 

l  l 
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when  w.  <  c. 
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Note  that  Equation  lUa  is  exactly  thi  same  as  Equation  12c  for  full  con¬ 
tact  and  is  used  to  start  the  iteration.  After  each  iteratior ,  a  check 
is  made  on  each  nodal  point  to  find  out  whether  any  contact  exists.  If 
the  deflection  w  is  smaller  than  ihe  initial  curling  c  ,  the  slab  is 
not  in  contact  with  the  subgrade,  and  the  subgrade  displacement  is  set 
to  zero,  as  indicated  in  Equation  lUb.  Thus  after  each  iteration,  a  new 
set  of  simultaneous  equations  is  established.  The  process  is  repeated 
until  the  same  equations  are  obtained.  In  most  cases,  this  can  be 


achieved  by  five  or  six  iterations.  After  the  deflections  due  to  the 
weight  and  warping  are  determined,  the  gaps  and  precompressions  can  be 
computed  and  used  later  for  computing  the  stresses  and  deflections  due 
to  the  load  alone. 

**6.  To  determine  the  stresses  and  deflections  due  to  the  load 
alone  requires  that  the  gaps  and  precompressions  shown  in  Figures  5b  and 
5c,  depending  on  whether  warping  exists,  first  be  determined.  When 
these  gaps  and  precompress  ions  are  used  as  s  ,  the  deflections  due  to 
the  load  alone  (Figure  5d)  can  be  determined  from  Equation  12,  except 
that  the  subgrade  displacements  are  expressed  as 


a: 

l 


when  w .  <  s . 

l  l 


(15a) 
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and  s .  >  0 
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when  w.  >  s.  and  s.  0 
ii  i 
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When  w  is  checked  with  s  ,  downward  deflection  is  considered  positive, 
upward  deflection  is  considered  negative,  gap  is  considered  positive, 
and  precompression  is  considered  negative.  First,  assume  that  the  slab 
and  the  subgrade  are  in  full  contact  and  the  deflections  of  tiie  slab  due 
to  the  applied  load  are  determined.  Then  check  the  deflections  with  s 
and  form  a  new  set  of  equations  based  on  Equation  l1'.  The  process  is 
repeated  until  the  same  equations  are  obtained. 

1*7.  When  the  slab  and  the  subgrade  are  in  partial  contact,  the 
principle  of  superposition  no  longer  applies.  To  ietermine  the  stresses 
and  deflections  due  to  an  applied  load  requires  that  the  deformed  shape 
of  the  slab  immediately  before  the  application  of  the  load  be  computed 
first.  Since  the  deformed  shape  depends  strongly  on  the  condition  of 
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warping,  the  stresses  and  deflections  due  to  loading  are  affected  apprec¬ 
iably  by  warping.  This  fact  was  borne  out  in  both  the  Maryland  (High¬ 
way  Research  Board  1952)  and  the  AASHO  (Highway  Research  Board  1962) 
road  tests. 

U8.  In  the  method  presented  here,  the  stresses  and  deflections 
due  to  weight  and  warping  are  computed  separately  from  those  due  to 
loading.  This  is  desirable  because  the  modulus  of  subgrade  reaction 
under  the  sustained  action  of  weight  and  warping  is  much  smaller  than 
that  under  the  transient  load  of  traffic.  If  the  same  modulus  of  sub¬ 
grade  reaction  is  used,  the  stresses  and  deflections  due  to  the  combined 
effect  of  weight,  warping,  and  loading  can  be  computed  in  the  same  way 
as  those  due  to  weight  and  warping,  except  that  additional  nodal  forces 
are  needed  to  account  for  the  applied  loads. 

1*9.  Partial  contact  with  initial  gaps.  This  case  applies  to  pave¬ 
ments  subjected  to  a  high  intensity  of  traffic,  such  as  the  traffic 
loops  in  the  AASHO  road  test  (Highway  Research  Board  1962).  Because  of 
pumping  or  plastic  deformation  of  the  subgrade,  some  springs  in  the 
Winkler  foundation  become  defective  and,  if  the  slab  is  removed,  will 
not  return  to  the  original  elevation.  Thus,  initial  gaps  are  formed,  as 
indicated  by  the  two  exterior  springs  shown  in  Figure  6a.  These  gaps  s 
must  be  assumed  before  an  analysis  can  be  made. 

50.  The  displacements  due  to  the  weight  of  the  slab,  as  shown  in 
Figure  6b,  can  be  determined  from  Equation  12a,  except  that  the  subgrade 
displacements  must  be  expressed  as 


s .  -  w. 
1  1 


when  w.  >  s. 

1  1 


(l6a) 


6'  =  <0  >  when  w^  < 


(16b) 


First,  assume  that  the  slab  and  the  subgrade  are  in  full  contact.  The 
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Figure  6.  Partial  contact  with  initial  gaps 


vertical  deflections  of  the  slab  are  determined  from  Equation  12a. 

Then  check  the  deflection  at  each  node  against  the  gap  s  .  If  the  de¬ 
flection  is  smaller  than  the  gap,  as  shown  by  the  left  spring  in  Figure 
6b,  Equation  l6b  is  used.  If  the  deflection  is  greater  than  the  gap, 
as  shown  by  the  other  springs  in  Figure  6b,  Equation  l6a  is  used.  The 
process  is  repeated  until  the  same  equations  are  obtained.  After  the 
deflections  are  obtained,  the  gaps  and  precompressions  can  be  computed 
and  used  later  for  computing  the  stresses  and  deflections  due  to  lead¬ 
ing,  if  no  warping  exists. 

51.  If  the  springs  are  of  the  same  length,  as  shown  in  Figure 
the  weight  of  the  slab  will  result  in  a  uniform  precompression,  and  no 
stresses  will  be  set  up  in  the  slab.  However,  if  the  springs  are  of 
unequal  lengths,  the  deflections  will  no  longer  be  uniform,  and  stress¬ 
ing  of  the  slab  will  occur. 

52.  Figure  6c  shows  the  combined  effect  of  weight  and  warping 

when  the  slab  is  warped  down.  The  reason  that  downward  warping  is  con¬ 
sidered  here,  instead  of  the  upward  warping,  is  that  the  case  of  upward 
warping  is  similar  to  that  shown  in  Figure  5c  except  that  the  gaps  are 
measured  from  the  top  of  the  defective  springs.  The  method  is  applic¬ 
able  to  both  upward  and  downward  warping,  but  downward  warping  is  used 
for  an  illustration.  The  procedure  for  determining  the  deflections  is 
similar  to  that  involving  the  weight  of  slab  alone  except  that  the  ini¬ 
tial  curling  of  the  slab,  as  indicated  by  Equation  11,  is  added  to  the 
gap  shown  in  Figure  6a  to  form  the  total  gap  and  precompression  s  for 

use  in  Equation  l6.  Since  the  gap  is  either  positive  or  zero  and  the 

initial  curling  may  be  positive  or  negative,  depending  on  whether  the 
slab  is  warped  up  or  down,  s  may  be  positive  or  negative.  After  the 
deflections  of  the  slab  are  obtained,  the  gaps  and  precompressions,  as 
shown  in  Figure  6c,  can  be  determined.  These  gaps  and  precompressions 
are  used  for  computing  the  stresses  and  deflections  due  to  the  load 
alone,  as  shown  in  Figure  6d. 

Symmetry 

53.  The  application  of  the  finite  element  method  for  analyzing 

rigid  pavements  involves  solving  a  large  set  of  simultaneous  equations. 
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However,  due  to  symmetry  the  number  of  simultaneous  equations  could 
be  greatly  reduced  by  considering  only  one  half  or  one  quarter  of  the 
slabs.  Consequently,  the  computer  time  and  storage  necessary  can  be 
greatly  reduced  and  the  results  yielded  are  the  same  as  for  the  use  of 
the  whole  slab. 

Stress  Transfer  Along  the  Joints  and  Clacks 

5h.  Joints  are  placed  in  rigid  pavements  to  control  cracking  and 
provide  enough  space  and  freedom  for  movement.  Load  is  transferred 
across  a  joint  or  crack  principally  by  shear  forces  and  in  some  cases  by 
moment  transfer.  Shear  force  is  provided  either  by  dowel  bars,  key 
joint,  or  aggregate  interlock.  Moment  transfer,  on  the  other  hand,  is 
provided  by  the  strength  of  the  concrete  slab  and/or  in-plane  thrust 
(which  is  ignored  in  this  analysis)  but  is  produced  by  heating  of  the 
slab.  When  a  joint  or  a  crack  has  a  visible  opening,  however,  the 
transfer  of  moment  across  the  joint  or  crack  becomes  negligible.  It  is 
therefore  justified  to  assume  there  is  no  moment  transfer  across  a  joint 
or  a  crack,  except  in  cases  such  as  a  tied  joint  where  some  moment  trans¬ 
fer  may  be  expected  if  the  joint  remains  tightly  closed. 

55.  In  a  continuously  reinforced  concrete  pavement  (CRCP),  a 
large  number  of  closely  spaced  cracks  may  form  soon  after  placing  the 
concrete.  The  tightly  closed  cracks  can  transfer  a  large  percentage  of 
moment  during  the  early  life  of  the  pavement,  but  as  traffic  load  is 
applied,  the  openings  of  the  cracks  will  increase  and  the  ability  of 
moment  transfer  across  the  cracks  will  decrease. 

56.  If  moment  transfer  across  a  joint  is  neglected,  the  amount 
of  stress  transfer  at  a  joint  is  governed  by  the  difference  in  deflec¬ 
tion  between  the  two  slabs  along  the  joint.  In  other  words,  the  shear 
transfer  is  100  percent  if  deflections  at  both  slabs  are  equal.  This 
difference  in  deflection  depends  on  the  shear  deformation  of  the  dowel 
bar  and  the  dowel-concrete  interaction.  The  analyses  presented  later 
indicate  that  the  effect  of  dowel-concrete  interaction  is  more  dominant 
than  the  shear  deformation  of  the  dowel  bar.  Neglecting  the  deformation 
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of  concrete  surrounding  the  dowel  bar  in  the  program  can  make  the  bar 
more  effective  than  in  reality.  Field  measurements  conducted  by  the 
Corps  of  Engineers  in  many  military  airfields  (Ohio  River  Division  Lab¬ 
oratories  1959)  indicated  that  the  dowel  bars  were  not  very  effective; 
the  average  stress  transfer  across  a  joint  was  only  about  25  percent. 
Detailed  discussion  in  this  respect  is  presented  in  Part  III. 

Shear  and  moment  trans- 
fer  across  joints  and  cracks 

57-  The  program  provides  three  options  for  specifying  shear  trans¬ 
fer,  but  only  one  for  moment  transfer.  The  three  options  for  shear 
transfer  are:  (a)  efficiency  of  shear  transfer,  (b)  spring  constant, 
and  (c)  diameter  and  spacing  of  dowels.  The  only  option  for  moment 
transfer  is  to  assume  an  efficiency  of  moment  transfer  across  the  joint. 
The  advantages  and  disadvantages  of  each  option  are  discussed  in  the 
next  section. 

58.  Efficiency  of  shear  transfer.  The  efficiency  of  shear  trans¬ 
fer  is  defined  as  the  ratio  of  vertical  deflections  along  the  joint  be¬ 
tween  the  unloaded,  or  less  heavily  loaded,  slab  and  the  adjacent  more 
heavily  loaded  slab.  This  is  the  easiest  method  to  specify  shear  trans¬ 
fer.  By  assigning  an  efficiency  between  0.0  (no  shear  transfer)  and 
1.0  (complete  shear  transfer),  reasonable  results  can  be  obtained  with 
a  minimum  number  of  iterations.  The  efficiency  can  be  easily  checked 
on  the  printout  by  comparing  each  pair  of  deflections  along  the  joint. 
However,  the  use  of  a  given  efficiency  for  all  nodes  along  a  joint  is 
not  realistic  because  the  deflection  ratios  in  an  actual  pavement  should 
vary  along  the  joint,  with  the  smallest  ratio  at  a  point  where  the  de¬ 
flection  is  the  largest.  To  determine  the  efficiency  in  the  field  by 
measuring  the  deflection  of  both  slabs  along  their  common  joint,  the 
ratio  at  the  point  of  largest  deflection  should  be  used,  thus  giving  a 
more  conservative  estimate  of  the  efficiency.  The  method  has  the  further 
disadvantage  that  the  slabs  must  be  numbered  according  to  the  magnitude 
of  load.  This  aspect  is  discussed  in  Report  2  of  this  series.  In  the 
computer  programs,  the  slabs  are  numbered  according  to  the  magnitude 
of  load. 
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59-  Spring  constant.  The  use  of  imaginary  shear  transfer  springs 
along  the  joint  between  two  slabs  to  determine  the  difference  in  deflec¬ 
tion  is  more  realistic  than  the  use  of  the  efficiency  of  load  transfer 
because  the  use  of  imaginary  springs  takes  into  consideration  the  shear 
force  at  the  joint.  The  spring  constant  is  defined  as  the  force  in 
pounds  per  linear  inch  to  cause  a  difference  in  deflection  of  1  in.;  so 
the  unit  of  the  spring  constant  is  in  pounts  per  inch  per  inch.  This 
option  can  be  used  either  for  key  joints  or  joints  with  aggregate  inter¬ 
lock  (such  as  tightly  closed  cracks)  or  joints  using  both  dowel  bars 
and  aggregate  interlock  for  shear  transfer.  The  spring  constant  can  be 
determined  in  deflections  at  a  number  of  points  along  the  joint.  The 
ratio  between  the  average  load  (applied  along  the  joint)  per  unit  width 
and  the  average  difference  in  deflection  is  the  spring  constant  of  the 
joint.  A  disadvantage  of  the  method  is  that,  without  any  test  data,  it 
is  very  difficult  to  assume  a  proper  spring  constant.  The  use  of  an 
improper  spring  constant  may  result  in  an  unreasonably  large  difference 
in  deflection  and  thus  require  a  large  number  of  iterations  to  obtain  a 
convergent  solution.  It  should  be  emphasized  that  the  spring  constant 
should  not  be  determined  based  on  test  data  of  a  single  wheel  load.  The 
ideal  test  procedure  is  to  place  a  long  piece  of  steel  beam  adjacent  to 
the  joint  and  to  apply  the  load  to  the  beam  through  a  series  of  wheel 
loads,  preferably  four  to  six  wheels. 

60.  Once  the  value  of  the  spring  constant  is  determined  for  a 
certain  joint  in  a  rigid  pavement,  the  difference  in  deflections  at 
each  nodal  point  across  the  joint  can  be  determined  based  on  the  shear 
forces  at  each  node  computed  at  the  particular  stage  during  the  itera¬ 
tion  cycles. 

61 .  Diameter  and  spacing  of  dowels.  This  method  is  most  straight¬ 
forward  and  does  not  require  a  field  test  to  determine  either  the  effi¬ 
ciency  of  shear  transfer  or  the  spring  constant.  This  option  takes  into 
consideration  the  diameter  and  spacing  of  dowels  or  steel  bars.  While 
this  option  should  yield  results  far  superior  to  the  other  two  options, 
the  method  has  some  disadvantages,  which  are  described  below. 

a^.  This  method  is  applicable  only  when  steel  bars  are  the 
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sole  means  of  shear  transfer.  In  a  tied  joint  that 
remains  closed,  the  shear  transfer  is  principally  pro¬ 
vided  by  the  granular  interlock,  not  by  the  thin  tie 
bars.  Therefore,  the  use  of  this  method  to  specify  the 
shear  transfer  of  a  tied  joint  would  result  in  an  effi¬ 
ciency  less  than  the  actual  value. 

b_.  This  method  requires  an  estimate  of  the  modulus  of  dowel 
support,  which  may  vary  considerably  depending  on  the 
type  of  dowel,  strength  of  concrete,  and  method  of  con¬ 
struction.  The  condition  of  the  dowel  bars  in  the  joint, 
i.e.,  the  degree  of  looseness  of  dowels  in  the  concrete, 
can  affect  the  modul  is  value  of  dowel  support.  In  Part  IV, 
computed  results  of  stress  transfer  across  the  joint  for 
many  military  airfield  pavements  are  presented.  The  re¬ 
sults  indicate  that  the  joint  performance  varies  greatly 
with  the  modulus  of  dowel  support  used  in  the  computation. 
Field  tests  conducted  in  a  number  of  military  airfields 
(Ohio  River  Division  Laboratories  1959)  indicated  that  in 
many  airfields  the  dowel  bars  in  the  concrete  were  loose 
(i.e.,  excessive  amounts  of  space  around  the  dowels  in 
the  concrete). 

62.  Efficiency  of  moment  transfer.  Analysis  of  moment  transfer 
across  a  joint  or  a  crack  is  more  involved  than  analysis  of  the  transfer 
of  shear  force.  The  amount  of  moment  transfer  depends  on  the  width  of 
the  crack,  thickness  of  the  slab,  amount  of  reinforcing  steel,  and  many 
other  factors,  and  is  difficult  to  analyze.  It  is  believed  that  for  a 
crack  with  a  visible  opening,  the  transfer  of  moment  is  negligible. 

63.  While  the  method  of  efficiency  of  shear  transfer  works  well 
in  specifying  shear  transfer  across  a  joint,  the  efficiency  of  moment 
transfer,  if  defined  as  the  rotation  ratio  between  the  unloaded  and 
loaded  slabs,  is  not  applicable  because  an  efficiency  of  zero  indicates 
a  zero  rotation  of  the  unloaded  slab  and  a  zero  rotation  is  not  realis¬ 
tic.  The  shear  forces  in  a  joint  not  only  cause  the  slabs  to  deflect, 
but  also  cause  the  slabs  to  rotate.  In  the  case  of  a  crack  with  a  vis¬ 
ible  opening,  the  rotation  at  the  unloaded  slab  is  not  zero;  the  joint 
acts  as  a  hinge  with  large  rotational  movements  at  the  joint.  Unlike 

a  zero  deflection,  a  zero  rotation  actually  requires  the  addition  of  a 
very  large  moment  in  the  unloaded  slab,  which  contradicts  the  defini¬ 
tion  of  zero  moment  transfer.  Because  the  moments  are  very  sensitive 
to  the  rotations,  and  could  cause  problems  in  solution  convergence,  it 
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is  thus  impractical  to  use  a  rotational  constant  similar  to  the  spring 
constant  to  determine  the  difference  in  rotation.  The  following  method 
for  moment  transfer  was  developed  and  found  to  be  satisfactory. 

6U.  Moment  transfer  across  the  joint  (or  crack)  is  specified  by 
the  efficiency  of  moment  transfer.  One  hundred  percent  moment  transfer 
is  defined  by  equal  rotations  at  nodal  points  at  both  sides  of  the  joint 
with  the  moments  computed  accordingly.  Zero  percent  moment  transfer  is 
defined  to  be  that  the  moments  of  nodal  points  along  the  joint  are  all 
zero  while  the  rotations  are  not  required  to  be  zero.  The  efficiency  is 
not  defined  as  the  rotation  ratio  between  the  unloaded  and  loaded  slabs, 
but  as  a  fraction  of  the  full  moment,  which  is  determined  by  assuming 
that  the  rotations  on  both  sides  of  the  joint  are  the  same.  Unless  the 
efficiency  of  moment  transfer  for  all  cracks  is  either  0.0  or  1.0,  it  is 
necessary  to  analyze  the  problem  twice.  First,  an  efficiency  of  1.0  is 
assumed  for  all  cracks  having  an  efficiency  other  than  zero,  and  the 
moments  at  each  node  along  the  crack  are  computed.  These  full  moments 
are  then  multiplied  by  the  efficiency  of  moment  transfer  at  the  corres¬ 
ponding  joint  to  determine  the  moments  that  actually  exist.  These 
moments  are  then  assigned  for  each  slab  edge,  as  externally  applied 
moments,  and  a  second  analysis  of  the  slabs  is  made. 

65.  It  can  be  seen  that  the  efficiency  of  moment  transfer  is  de¬ 
fined  differently  than  the  efficiency  of  shear  transfer.  The  efficiency 
of  shear  transfer  is  based  on  vertical  deflections  instead  of  vertical 
forces,  whereas  the  efficiency  of  moment  transfer  is  based  on  moments, 
instead  of  rotations.  It  is  possible  to  define  shear  transfer  on  the 
basis  of  vertical  forces,  in  the  same  way  as  for  moment  transfer,  and 
analyze  the  problem  twice.  However,  this  is  not  warranted  because  many 
of  the  practical  problems  involve  shear  transfer  only  and  can  be  solved 
in  one  analysis  when  shear  transfer  is  defined  by  vertical  deflections. 
Also,  in  the  efficiency  of  shear  transfer,  the  vertical  deflections  on 
both  sides  of  the  slab  are  different,  but  in  the  efficiency  of  moment 
transfer,  the  moments  on  both  sides  of  the  siabs  are  the  same. 
Computations  for  dowel  bars 

66.  For  stress  transfer  using  the  method  of  diameter  and  spacing 
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of  the  dowel  bars,  the  derivation  of  equations  is  presented  in  the  fol¬ 
lowing  paragraphs . 

67.  Figure  7a  shows  an  exaggerated  view  of  a  doweled  joint  with 
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H  K- 


(a)  EXAGGERATED  CROSS-  (b)  VERTICAL  SHEAR  FORCE 

SECTION  VIEW  OF  DOWEL  TRANSMITTED  TO  DOWEL 

JOINT 


Ic)  ILLUSTRATION  OF  DIFFERENCES  IN  DEFLECTION 

Figure  7.  Difference  in  deflection  between  two  slabs 


a  width  d  .  In  the  finite  element  method,  the  slabs  are  divided  into 
a  number  of  elements.  Suppose  that  1  and  2  are  a  pair  of  the  many  nodal 
points.  If  a  vertical  shear  force  P  exists  at  node  1,  it  will  be 
transmitted  to  the  dowel,  as  shown  in  Figure  7b.  After  deformation,  the 
deflections  at  nodes  1  and  2  will  not  be  the  same  and  the  difference 
in  deflection  A  between  nodes  1  and  2,  as  shown  in  Figure  7c,  is 
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(IT) 


A  =  A  +  2A 
s  c 


where 

A  =  shear  deformation  of  the  dowel 
s 

A  =  shear  deformation  of  the  concrete  due  to  the  bending  of  the 
dowel  as  a  beam  on  an  elastic  (concrete)  foundation 

68.  The  shear  deformation  of  the  dowel  can  be  determined  approx¬ 


imately  by 


(18) 


in  which  G  and  A  are  the  shear  modulus  and  the  area  of  the  dowel, 
respectively.  By  considering  the  dowel  as  a  beam  on  an  elastic  founda¬ 
tion,  A^  can  be  determined  by  (Timoshenko  and  Lessels  1925,  Yoder  and 
Witczak  1975) 


A  =  — | —  (2  +  8d ) 
C  he3EI 


(19a) 


(19b) 


where 

E  =  modulus  of  elasticity  of  the  dowel 
I  =  moment  of  inertia  of  the  dowel 
K  =  modulus  of  dowel  support,  pci 
b  =  diameter  of  the  dowel 

By  using  realistic  values  for  the  parameters  in  Equations  18,  19a,  and 
19b,  it  can  be  shown  that  Ag  is  one  or  more  orders  of  magnitude  smaller 
than  Ac  and  can  usually  be  neglected;  thus,  the  difference  in  deflec¬ 
tion  depends  principally  on  the  dowel-concrete  interaction.  Because  the 
amount  of  stress  transfer  across  the  joint  is  governed  by  the  difference 
in  deflection,  the  negligence  of  the  dowel-concretQ  interaction  will  re¬ 
sult  in  a  stress  transfer  that  is  too  large. 

69.  It  should  be  pointed  out  that  the  dowel  bars  are  not  modeled 
as  bar  elements  in  the  programs.  The  joints  or  cracks  in  a  rigid  pave¬ 
ment  generally  have  a  width  of  l/l6  to  1/32  in.,  which  is  much  smaller 
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than  the  bar  diameter.  A  dowel  bar  in  a  joint  should  therefore  serve 
merely  as  a  key  to  transmit  shear  forces.  If  it  is  desirable  to  model 
dowels  as  bar  elements,  meaningful  results  can  be  obtained  only  if  the 
theory  of  the  deep  beam  is  employed  in  the  analysis,  which,  of  course, 
is  very  complex.  Tabatabaie  and  Barenberg  (1978)  modeled  dowel  bars  as 
bar  elements  in  their  work;  a  detailed  discussion  in  this  aspect  was 
made  by  Huang  and  Chou  (1978). 

Modulus  of  dowel  support 

70.  The  stresses  in  dowel  bars  result  from  shear,  bending,  and 
bearing  forces.  These  stresses  can  be  analyzed  to  determine  factors 
that  affect  the  load  transfer  characteristics.  The  stress  analysis  of 
dowels  is  based  upon  work  by  Timoshenko  (Timoshenko  and  Lessels  1925); 
Timonshenko  modeled  a  dowel  bar  encased  in  concrete  as  a  beam  on  a  Wink¬ 
ler  foundation.  The  ratio  between  the  bearing  pressure  and  the  deflec¬ 
tion  of  the  dowel  bar  was  termed  as  modulus  of  dowel  support  K  . 

Table  1  shows  a  wide  range  of  moduli  of  dowel  support  produced  by  many 
investigators.  The  information  was  gathered  by  Finney  (1956).  Finney 
reported  that  while  testing  procedures  varied  among  investigators,  K 
values  also  varied  between  specimens  for  a  given  test  procedure.  A 
study  of  these  investigations  seems  to  indicate  that  K  is  not  a  con¬ 
stant  quantity,  but  varies  with  the  concrete  properties,  dowel  bar  di¬ 
ameter  and  length,  slab  thickness,  and  the  degree  of  dowel  looseness  in 
the  concrete.*  Yoder  and  Witczak  (1975)  suggested  that  values  of  K 
range  between  300,000  and  1,500,000  pci  and  that  the  use  of  1,500,000  pci 
appears  to  be  warranted.  It  was  found  in  this  study  that  while  large 
changes  in  the  modulus  do  not  affect  £he  stress  calculations  in  the 

slab  greatly,  they  do  affect  the  deflection  of  the  dowel  which  in  turn 
can  cause  the  change  in  the  stress  transfer  across  a  joint. 

71.  To  account  in  the  computer  program  for  the  possible  loose¬ 
ness  of  dowels,  two  different  moduli  of  dowel  support  K  can  be  speci¬ 
fied.  (This  is  done  only  in  the  WESLIQID  program.)  When  A  <  ts'  , 

s  —  c 

K  =  ;  when  A^  >  A'  ,  K  =  Kg  ;  where  A'  =  an  input  parameter 

*  Slight  looseness  on  one  side  of  the  dowel  is  often  intentionally 
built  into  the  system  in  construction  practice  of  rigid  pavements. 
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specifying  the  deformation  of  concrete  below  which  modulus  is 

used  and  above  which  modulus  is  used.  If  a  very  large  A'  is 

assigned  in  the  input,  only  will  be  used  in  the  computation. 

72.  If  Ac  ,  as  computed  from  Equation  19a  and  based  on  K  = 
is  not  greater  than  A'  ,  it  will  be  used  in  Equation  17  to  determine 
the  difference  in  deflection.  Otherwise,  the  following  equations  are 
used: 


4B3EIA' 

_ c 

2  +  Bd 


(20) 


in  which  P"  =  shear  force  on  dowel  to  affect  A'  ,  and  B  can  be 
determined  from  Equation  19b  with  K  =  . 

A  =  A'  +  F-~  P-  (2  +  Bd) 
c  c  Ub3EI 

In  Equation  21,  K 2  is  used  to  determine  B  . 


(21) 


PART  III:  PRESENTATION  OF  NUMERICAL  RESULTS 
FOR  THE  WESLIQID  PROGRAM 


Comparison  with  Available  Solutions 

73.  To  check  the  accuracy  of  the  finite  element  method  and  the 
correctness  of  the  computer  program,  it  is  desirable  to  compare  the 
finite  element  solutions  with  other  theoretical  solutions  available, 
especially  with  those  involving  discontinuities  such  as  the  stress  at 
the  free  edge  of  a  slab.  Westergaard' s  original  work  and  Pickett  and 
Ray's  influence  charts  can  be  used  for  such  purposes.  Because  these 
available  solutions  are  based  on  an  infinite  slab,  a  very  large  slab 
was  used  in  the  finite  element  analysis. 

Westergaard1 s  solutions 

Jh.  The  finite  element  solutions  were  obtained  by  using  a  large 
slab,  201  long  by  10£  wide,  where  £  is  the  radius  of  relative  stiff¬ 
ness. 


where 

E  =  Young's  modulus  of  the  pavement 
t  =  thickness  of  the  pavement 
v  =  Poisson's  ratio  of  the  pavement 
k  =  modulus  of  subgrade  reaction 

Because  the  problem  is  symmetrical  with  respect  to  the  y-axis,  only  one 
half  of  the  slab  was  considered.  The  slab  was  divided  into  rectangular 
finite  elements  as  shown  in  Figure  8.  Both  the  x  and  y  coordinates  are 
0  ,  TT  1/8  ,  T\l/h  ,  TT8./2  ,  3TT5./14  ,  U  l  ,  $1  ,  6i  ,  7 1  ,  8  S.  ,  9«.  , 
and  10 l  .  (Note  the  elements  are  not  equally  spaced  for  x  and  y 
smaller  than  hi  .)  The  Poisson's  ratio  of  the  concrete  is  0.25.  The 
finite  elements  are  so  divided  because  the  Westergaard 's  solution  for 
the  problem  at  these  coordinates  is  available  in  Westergaard  (1925). 

75.  Figures  8  and  9  show  the  comparison  between  Westergaard' s 
exact  solutions  and  the  WESLIQID  finite  element  solution  for  deflection 

hi 


Distances  along  X-aris 

Figure  9.  Comparison  of  WESLIQID  finite  element  solution 
with  Westergaard' s  solution,  moments 


and  moment,  respectively.  The  former  is  for  an  infinite  slab  with  a 
point  load  P  acting  at  one  edge  far  from  any  corner.  The  Westergaard 
solutions  are  indicated  by  the  solid  curves,  and  the  finite  element 
solutions  are  indicated  by  the  small  circles.  The  bonding  moments  M 
about  the  y-axis  are  computed  along  the  edge  of  the  slab,  and  the  de¬ 
flections  w  are  along  the  edge  and  at  distances  vi./k  and  ttJI/2 
from  the  edge.  The  Westergaard  solutions  for  moment  and  deflection  are 
obtained  from  those  shown  in  Figures  11  and  8,  respectively,  of  Wester¬ 
gaard  (1925).  It  can  be  seen  that  the  finite  element  solutions  check 
very  closely  with  Westergaard' s  results. 

76.  Westergaard  computed  critical  stresses  in  concrete  pavements 
under  three  types  of  loading  conditions.  Type  1  is  a  wheel  load  acting 
close  to  a  rectangular  corner  of  a  la:  ge  panel  of  the  slab.  The  criti¬ 
cal  stress  is  a  tension  at  the  top  of  the  slab.  Type  2  is  a  wheel  load 
acting  at  a  considerable  distance  from  the  edges  and  is  generally  called 
a  center  load.  The  critical  tension  occurs  at  the  bottom  of  the  slab 
under  the  center  of  the  load.  Type  3  is  a  wheel  load  acting  at  the 
edge  of  the  slab  but  at  a  considerable  distance  from  any  corner;  it  is 
generally  called  an  edge  load.  The  critical  stress  is  a  tension  at  the 
bottom  under  the  center  of  the  loaded  circle. 

77-  The  WESLIQID  finite  element  program  was  used  to  compute 
stresses  in  concrete  pavements  for  these  three  types  of  load.  Point 
loads  were  used  because  they  are  easier  to  work  with  in  the  finite  ele¬ 
ment  method.  The  pavement  was  assumed  to  be  10  in.  thick  in  all  calcula¬ 
tions,  and  the  modulus  of  subgrade  reaction  had  three  different  values. 
The  computed  maximum  stresses  are  tabulated  in  Table  2  together  with 
those  computed  values  obtained  by  the  Westergaard  method  (1925). 

78.  In  the  cases  of  interior  and  edge  load,  the  stresses  com¬ 
puted  by  the  finite  element  method  are  slightly  higher  than  those  of 
the  Westergaard  solution.  This  could  be  attributed  to  the  fact  that 
the  Westergaard  solution  uses  a  semi-inifinite  slab,  while  the  finite 
element  method  employs  a  finite-size  slab.  It  is  believed  that  if  the 
slab  size  is  increased  in  the  finite  element  method,  the  computed 
stresses  will  be  reduced  slightly.  In  the  case  of  corner  load,  the 


stresses  computed  by  the  finite  element  method  are  considerably  less 
than  the  300  psi  computed  by  the  Westergaard  solution.  The  explanation 
of  this  discrepancy  is  that  the  corner  stress  computed  by  the  Wester¬ 
gaard  solution  is  the  maximum  stress  under  the  point  load,  but  the 
corner  stresses  computed  by  the  finite  element  method  as  given  in 
Table  2  are  not  the  maximum  values  but  the  stresses  computed  at  node  11 
in  Figure  2.  More  discussion  in  this  respect  is  given  in  the  next 
paragraph. 

T9-  Figure  2  shows  the  finite  element  layout  of  a  concrete  slab 
subjected  to  a  point  load  acting  at  the  corner  of  the  slab.  The  stresses 
computed  at  the  nodes  surrounding  the  load  are  presented  in  Table  3. 

It  is  seen  that  node  11  has  the  highest  stress  among  those  nodes  where 
stresses  were  computed,  but  the  stress  at  node  11  is  not  necessarily 
the  maximum  in  the  slab  under  the  load.  The  maximum  stress  in  the  slab 
lies  possibly  between  nodes  1  and  11  (or  between  11  and  12,  21,  or  20), 
depending  on  the  modulus  of  subgrade  reaction,  and  it  can  be  determined 
only  by  further  dividing  the  finite  element  grid.  Since  the  stresses 
under  the  corner  load  computed  by  the  finite  element  method  in  Table  2 
are  not  the  maximum  stresses,  they  are  therefore  smaller  than  the  maxi¬ 
mum  stresses  computed  by  the  Westergaard  solution.  The  discrepancy  is 
greater  for  a  greater  modulus  of  subgrade  reaction.  In  Figure  2,  node 
11  is  28.3  in.  away  from  the  corner;  it  is  believed  that  when  the  mod¬ 
ulus  of  the  subgrade  reaction  increases,  the  location  of  the  maximum 
stress  in  the  slab  would  be  closer  to  the  corner  where  the  load  is 
applied. 

80.  Influence  charts  by  Pickett  and  Ray  (1951)  were  not  used  di¬ 
rectly  to  check  the  finite  element  results.  Instead,  the  results  of  the 
WESLIQID  program  were  checked  with  those  computed  by  the  H-51  computer 
program.  H-51  was  developed  by  the  General  Dynamics  Corporation  to  com¬ 
pute  edge  stress  in  a  concrete  pavement  based  on  the  Westergaard  solu¬ 
tion.  The  program  has  been  used  by  the  Corps  of  Engineers  for  several 
years  to  determine  edge  stresses  under  multiple-wheel  loads.  It  was 
found  that  the  edge  stresses  computed  by  H-51  compare  very  closely  with 
those  from  Pickett  and  Ray's  influence  charts  (1951). 


8l.  Figure  10  shows  the  layouts  of  the  finite  element  and  H-51 
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Figure  10.  Computation  of  edge  stresses  using  WESLIQID 
finite  element  and  H-51  program 

programs  to  compute  the  edge  stresses  under  a  50,000-11)  wheel  load.  In 
the  finite  element  method,  a  20-ft  by  20- ft  slab  was  used,  and  because 
of  symmetry  only  one  half  of  the  slab  and  a  25, 000- lb  load  were  used  in 
the  computation.  The  maximum  radial  tensile  stress  computed  by  the 
H-51  program  was  607.3  psi,  and  that  computed  by  the  finite  element  pro¬ 
gram  was  624.3  psi,  a  difference  of  2.7  percent. 

Discrete  element  solution 

82.  The  finite  element  solutions  were  also  compared  with  the 
discrete  element  solution  developed  at  the  University  of  Texas.  The 
program  is  called  SLAB30  (Hudson  and  Matlock  1965).  The  discrete  ele¬ 
ment  approach  is  mathematically  equivalent  to  the  finite  difference 
method.  The  pavement  slab  is  represented  as  a  combination  of  elastic 
blocks,  rigid  bars,  and  torsion  bars.  Figure  11  shows  the  layout  of 


the  finite  element  solution.  The  slab  has  four  edges  with  four  concen¬ 
trated  loads  acting  av  nodes  18,  19*  20,  and  21.  The  same  slab  with 
the  same  loading  condition  was  used  for  the  CLAB30  program,  except  that 
the  slab  was  divided  with  equal  increments  of  1  ft  in  both  x-  and  y- 
directions.  Table  h  presents  the  comparisons  of  the  computed  stresses 
ana  deflections  at  the  nodal  points  designated  in  Figure  11.  It  can  be 


Figure  11.  Computations  of  stresses  and  deflections  using 
the  WESLIQID  finite  element  program,  single  slab 

seen  that  excellent  agreement  was  obtained  on  deflections  but  that  the 
stresses  agree  only  at  the  interior  nodes.  For  nodes  along  the  edge, 
the  absolute  values  of  the  stresses  computed  by  the  discrete-element 
program  are  generally  much  smaller  than  those  computed  by  the  finite- 
element  program.  Theoretically,  the  stresses  at  the  corner  of  a  slab 
should  be  zero,  but  fairly  large  stresses  are  computed  at  the  corner 
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node  (node  ll)  by  the  discrete-element  program.  The  stresses  com¬ 
puted  at  the  corner  and  edge  nodes  by  the  discrete  element  program  are 
believed  to  be  wrong.  The  errors  are  possibly  caused  by  the  practice 
of  reducing  the  bending  stiffness  of  the  edge  and  corner  nodes  to  one- 
half  and  one-fourth  of  the  full  bending  stiffness. 

83.  To  further  demonstrate  the  characteristics  of  stress  dis¬ 
tribution  computed  by  the  SLAB30  program,  comp\itations  are  made  for  a 
26-ft  by  26-ft  concrete  slab  subjected  to  two  concentrated  loads  placed 
at  the  edge  but  away  from  the  corners  of  the  slab.  The  maximum  princi¬ 
pal  stresses  are  plotted  along  the  center  line  of  the  slab,  as  shown  in 
Figure  12.  Theoretically,  maximum  tensile  stress  should  occur  at  the 
loaded  edge  of  the  slab  and  in  the  direction  parallel  to  pavement  edge, 
and  gradually  reduce  its  magnitude  toward  the  interior  of  the  slab. 

The  results  plotted  in  Figure  12  show  that  the  stress  at  the  edge 
(365  psi)  is  smaller  than  the  stress  at  the  node  1  ft  away  from  the 
edge  (^98  psi),  which  obviously  is  not  correct. 


Figure  12.  Stress  distribution  in  a  concrete  pavement  under  two 
concentrated  loads  computed  by  SLAB30  program 


Comparison  with  Experimental  Results 


AASHO  road  tests 

8U.  Efforts  were  made  to  compare  the  finite  element  solution 
with  the  strain  measurements  from  the  (Highway  Research  Board  1962 ) 

AASHO  road  tests.  Because  a  similar  comparison  had  been  made  by 
Tabatabaie  and  Barenberg  (1978) ,  also  using  the  two-dimensional  finite 
element  plate-bending  model,  the  same  test  section,  material  properties, 
and  loads  were  thus  used  in  the  computation.  The  following  information 
is  excerpted  from  Tabatabaie  and  Barenberg: 

Tests  were  conducted  on  the  main  traffic  loops 
where  the  strain  due  to  moving  traffic  12  to  22  in. 
from  the  edge  was  measured  at  the  slab  edge.  The 
lengths  of  slabs  were  15-ft  nonreinforced  sections 
and  UO-ft  reinforced  slabs.  Slab  thickness  ranged 
from  5  to  12.5  in.  The  measured  dynamic  modulus 
of  elasticity  and  Poisson's  ratio  of  concrete  were 
found  to  be  6.25  *  ic£  psi  and  0.28,  respectively. 

The  modulus  of  subgrade  reaction  (k  values)  on  the 
subbase  obtained  by  the  plate-bearing  tests  varied 
from  approximately  85  to  200  pci  over  all  of  the 
loops  throughout  the  two-year  test  period.  An 
average  of  150  pci  was  used  for  modulus  of  subgrade 
reaction  in  the  finite  element  analysis. 

85.  Figure  13  shows  the  comparison  of  the  WESLIQID  finite  ele¬ 
ment  solutions  with  AASHO  experimental  results  and  the  calculations  by 
Tabatabaie  and  Barenberg  (1978).  Except  for  the  results  computed  by  the 
WESLIQID  program,  the  plotted  results  shown  in  Figure  13  are  taken  di¬ 
rectly  from  Tabatabaie  and  Barenberg.  It  is  seen  that  excellent  agree¬ 
ment  is  obtained  in  the  theoretical  solutions  computed  by  the  program 
prepared  by  Tabatabaie  and  Barenberg  and  the  WESLIQID  program  prepared 
by  the  WES.  This  agreement  is  not  surprising  because  both  programs 
were  developed  based  on  the  plate-bending  theory.  The  computed  results 
should  be  very  close  when  the  problem  is  dealing  with  a  single  slab 
that  has  full  contact  with  the  subgrade.  The  computed  results  may 
differ  when  joint  and  partial  contact  are  considered. 
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Figure  13.  Comparison  of  edge  stresses  computed  with 
WESLIQID  finite  element  program  and  those  measured  at 
the  AASHO  road  test 


Strain  measurements  from 
the  Corps  of  Engineers  tests 

86.  Description  of  strain  measurements.  During  the  1950' s,  a 
series  of  full-scale  static  loading  tests  on  rigid  pavements  was  con¬ 
ducted  at  eight  Air  Force  bases  and  at  the  Ohio  River  Division  Labora¬ 
tories  (ORDL)  (1959).  The  purpose  of  such  tests  was  to  evaluate  dif¬ 
ferent  methods  of  dowel  installation  in  terms  of  the  load  transfer 


achieved  at  the  joint  and  to  check  whether  any  of  the  methods  failed  to 
provide  a  joint  meeting  the  minimum  requirements  of  25  percent  load 
transfer  as  assumed  in  the  design  and  evaluation  of  rigid  pavements. 

87.  Doweled  joints  investigated  for  the  study  were  constructed 
hy  five  different  methods,  which  included  (a)  cast  in  place,  (b)  remove 
and  replace,  (c)  split  dowels  (screw-type  only),  (d)  dummy  half-dowel, 
and  (e)  oversize  dummy  half-dowel  (cement  grout  only).  Aircraft  used 
in  the  loading  tests  included  B-U7,  B-52,  and  KC-97  aircraft.  The  test 
cart  used  on  the  ORDL  experimental  pavement  had  a  four-wheel  twin- 
tandem  configuration  similar  to  that  used  as  the  main  gear  for  the  B-36 
aircraft.  Joint  performance  was  evaluated  on  the  basis  of  strain  mea¬ 
surements  at  the  pavement  surface  in  the  direction  parallel  to  the  joint. 
Test  procedures  included  the  no  load ,  load,  and  rebound  measurements. 

88.  The  description  of  fc st  sites,  instrumentation,  test  proce¬ 
dures,  and  test  results  of  the  study  are  documented  in  Ohio  River  Divi¬ 
sion  Laboratories  (1959).  Table  5  is  a  summary  of  pavement,  joint,  and 
subgrade  characteristics  pertinent  to  the  evaluation.  Table  6  presents 
a  summary  of  the  average  strain  measurements  obtained  from  the  loading 
tests  at  each  of  the  field  test  sites.  Comparative  average  strains  are 
shown  for  the  load  wheels  on  the  bonded  and  unbonded  sides  of  the  dowel, 
as  well  as  for  the  load  positioned  either  directly  over  a  dowel  or  mid¬ 
way  between  adjacent  dowels.  At  ORDL  Test  Track  "A,"  no  mine  detector 
was  available  for  locating  the  dowels.  Therefore,  the  data  for  this 
test  area  were  the  average  strains  for  all  measurements  without  regard 
to  the  relative  positioning  of  the  strain  gages  and  the  dowels.  Table  7 
presents  a  summary  of  the  average  amount  of  load  transfer  measured  at 
each  of  the  nine  test  areas.  The  load  transfer  is  defined  as  the  ratio 
of  the  strain  on  the  unloaded  side  of  joint  to  that  of  the  total  strain 
(the  sum  of  the  strains  on  both  unloaded  and  loaded  sides)  expressed  as 
a  percentage.  Further  discussion  on  the  adequacy  of  the  definition  of 
"load  transfer"  is  given  in  paragraph  93. 

89.  It  was  pointed  out  in  Ohio  River  Division  Laboratories 
(1959)  that  based  on  Westergaard ’ s  theoretical  analysis  of  rigid  pave¬ 
ments,  the  sum  of  the  strains  on  the  loaded  and  unloaded  sides  of  a 
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doweled  joint  is  equal  to  the  strain  produced  by  loading  tangentially 
to  a  free  edge  of  the  pavement.  It  was  thus  valid  to  measure  strains 
simultaneously  on  both  sides  of  the  joint  to  permit  a  reasonable  deter¬ 
mination  of  the  load  transfer  characteristics  of  doweled  joints.  This 
is  verified  by  the  results  computed  by  the  WESLIQID  program  that  the 
sum  of  the  strains  on  both  sides  of  the  joint  is  nearly  a  constant  and 
is  equal  to  the  strain  at  the  free  edge. 

90.  Table  7  indicates  that  the  load  transfer  produced  in  the 
nine  test  areas  ranges  from  a  minimum  of  11. 5  percent  to  a  maximum  of 
k0.2  percent.  The  average  load  transfer  for  all  test  siwes  was  28.0 
percent . 

91.  It  was  suggested  by  the  Ohio  River  Division  Laboratories 
(1959)  that  the  most  improtant  single  factor  affecting  the  load  trans¬ 
fer  capacity  of  the  joint  is  the  amount  of  looseness  that  exists  between 
the  dowel  and  the  pavement.  The  initial  looseness  is  probably  related 
to  the  particular  construction  procedure  used  in  the  installation  of  the 
dowel.  During  the  service  life  of  the  pavement,  some  additional  loose¬ 
ness  will  develop  as  a  result  of  the  gradual  enlarging  of  the  dowel 
socket  under  the  action  of  repetitive  loading.  Further  loss  of  load 
transfer  capacity  in  a  doweled  joint  can  also  result  from  a  reduction 

in  the  effective  cross-sectional  area  of  the  dowels  due  to  corrosion. 

The  study  further  concluded  that  the  load  transfer  characteristics  of 
doweled  longitudinal  construction  joints  are  not  affected,  regardless 
of  which  side  of  the  joint  is  loaded. 

92.  Computations  by  the  WESLIQID  program.  The  WESLIQID  program 
was  used  to  analyze  the  nine  test  sites  listed  in  Table  5.  Pavement 
and  dowel  information  given  in  Table  t  was  not  sufficient  to  carry  out 
the  computations;  assumed  values  for  some  variables  had  to  be  used. 

The  modulus  of  the  concrete  E  ,  joint  spacing  d  ,  and  modulus  of 
dowel  support  K  were  at  first  assumed  to  be  6, 000,000  psi ,  1/32  in., 
and  1,500,000  pci,  respectively.  However,  all  these  parameters  were 
varied  in  the  case  of  the  computation  given  in  Table  8.  The  values  of 
the  dowel  support  K  take  into  account  the  different  construction 
method  of  the  doweled  joint  and  the  condition  of  the  joint.  The  computed 


51 


strains  and  load  transfer  were  generally  close  to  those  measured  values. 
It  was  found  also  that  the  computed  values  were  influenced  greatly  with 
changes  of  E  and  K  ,  but  insignificantly  with  the  change  of  joint 
spacing  d  .  Increasing  the  E  value  decreases  the  strains  on  both 
sides  of  the  joint  and  vice  versa.  Since  the  changes  are  nearly  the 
same  on  both  sides  of  the  joint,  the  load  transfer  is  not  influenced 
significantly  by  the  changes  of  E  values.  On  the  other  hand,  increas¬ 
ing  the  modulus  of  dowel  support  K  makes  the  concrete  surrounding  the 
dowel  bars  stiffer,  reduces  the  strains  at  the  loaded  side,  and  in¬ 
creases  the  strains  at  the  unloaded  side  of  the  joint.  The  combined 
effect  significantly  increases  the  load  transfer  of  the  joint. 

93.  Table  8  shows  the  measured  and  computed  strains  and  load 
transfers  for  the  nine  test  sites.  The  computed  stresses  in  the 

direction  parallel  to  the  joint  are  also  shown  in  the  table.  It  should 
be  noted  that  in  the  theoretical  computations,  the  load  is  placed  di¬ 
rectly  next  to  the  joint  and  the  strains  are  computed  at  the  joint.  In 
the  actual  field  measurements,  however,  the  strain  gages  were  placed  a 
small  distance  away  from  the  joint,  and  the  load  is  then  placed  a  small 
distance  from  the  gage.  Therefore,  the  actual  strains  measured  should 
be  slightly  less  than  those  computed.  Also  in  Table  8,  the  term  "load 
transfer"  is  changed  to  "stress  transfer."  Current  design  criteria  for 
rigid  pavements  for  military  airfields  are  based  on  either  the  premise 
that  all  types  of  joints  transfer  a  minimum  of  25  percent  of  the  load 
applied  along  the  joint  to  the  adjacent  slab.  The  stresses  along  the 
joints  with  adequate  dowel  construction  are  calculated  as  75  percent 
of  the  free  edge  stresses  with  the  load  applied  at  the  edge  of  the 
pavement.  It  is  believed  that  the  term  "stress  transfer"  describes  the 
situation  better  than  the  term  "load  transfer." 

9^.  The  comparisons  between  the  computed  and  measured  values  in 
Table  8  are  discussed  separately  for  each  test  site  as  follows: 

a_.  Lockbourne  AFB,  Ohio,  Entry  1-a  shows  the  measured 

strains  and  stress  transfer.  The  values  shown  in  entry 
1-b  are  computed  with  E  ,  K  ,  and  d  equal  to 
6,000,000  psi,  1,500,000  pci,  and  1/32  in.,  respectively. 
The  stress  transfer  is  6  percent  higher  than  the  measured 
value,  and  the  computed  strains  are  much  higher  than 


52 


those  measured.  Entry  1-c  shows  the  values  computed 
with  E  =  8,000,000  psi.  It  is  seen  that  the  strain 
values  are  reduced  considerably  and  are  very  close  to 
those  measured  values.  Also,  the  stress  transfer  value 
is  reduced  by  about  2  percent.  Entry  1-d  shows  the  com¬ 
puted  values  by  keeping  the  value  of  E  to  6,000,000  psi 
but  reducing  the  modulus  of  dowel  support  K  to  300,000 
pci.  It  is  seen  that  while  the  stress  transfer  is  further 
reduced  (2  percent  less  than  measured  value),  the  strain 
values  increased  considerably,  particularly  on  the  loaded 
side.  In  other  words,  when  the  modulus  of  dowel  support 
is  reduced,  the  strains  in  the  concrete  increase;  and 
since  the  strain  in  the  loaded  side  increases  much  more 
than  that  in  the  unloaded  side,  the  percentage  of  stress 
transfer  is  thus  decreased.  To  check  the  effect  of  the 
crack  opening  d  on  the  computed  values,*  a  value  of 
1/2  in.  was  used  with  E  and  K  equal  to  6,000,000  psi 
and  1,500,000  pci,  respectively.  The  computed  values  are 
tabulated  in  entry  1-e.  It  is  seen  that  the  computed 
values  are  insignificantly  different  from  those  in  entry 
1-b,  indicating  that  the  crack  opening  has  very  little 
influence  on  the  joint  performance. 

b.  Lincoln  AFB,  Nebraska.  Entry  2-a  shows  the  measured 

strains  and  stress  transfer.  Entry  2-b  shows  the  values 
computed  with  the  E  ,  K  ,  and  d  values  equal  to 
6,000,000  psi,  1,500,000  pci,  and  1/32  in.,  respectively. 
It  is  seen  that  the  difference  in  total  strain  is  very 
small  between  the  measured  and  computed  values,  but  there 
is  a  9-percent  difference  in  the  stress  transfer,  result¬ 
ing  from  too  large  a  computed  strain  in  the  loaded  side 
and  too  small  a  computed  strain  in  the  unloaded  side. 

By  increasing  the  modulus  of  dowel  support  K  from 
1,500,000  to  9,000,000  pci  (i.e.,  the  concrete  surrounding 
the  dowel  becomes  stiffer  or  the  looseness  in  the  doweled 
joint  is  reduced),  the  strain  in  the  loaded  side  decreases 
noticeably  and  the  strain  in  the  unloaded  side  increases 
slightly.  The  combined  effect  results  in  a  computed 
stress  transfer  very  close  to  the  measured  value,  although 
the  computed  total  strain  becomes  slightly  less  than  the 
measured.  The  good  comparison  between  the  computed 
(entry  2-c )  and  measured  values  suggests  that  the  dowel 
bars  had  a  very  good  and  tight  fit  in  the  joint  in  the 
rigid  pavements  tested  at  Lincoln  APB. 

£_•  Hunter  AFB,  Georgia.  Entry  3-a  shows  the  measured 
strains  and  stress  transfer.  With  E  ,  K  ,  and  d 
values  of  6,000,000  psi,  1,500,000  pci,  and  1/32  in., 


*  In  this  report,  the  crack  opening  is  treated  the  same  way  as  the 
joint  spacing  in  the  computer  programs. 
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respectively,  the  computed  values  shown  in  entry  3-b 
are  very  close  to  measured  values.  Similar  to  the  ex¬ 
planation  Riven  for  entry  2-c,  it  is  believed  that  by 
increasing  slightly  the  modulus  of  dowel  support  K  , 
the  computed  values  can  be  closer  to  those  measured. 

d. .  McCoy  AFB,  Florida.  The  values  shown  in  entry  L-b  are 

computed  with  E  ,  K  ,  and  d  equal  to  6,000,000  psi, 
1,500,000  pci,  and  1/32  in.,  respectively.  It  is  seen 
that  the  computed  stress  transfer  is  very  close  to  the 
measured  value,  but  the  computed  strains  are  slightly 
higher  than  those  measured.  The  differences  can  be 
readily  corrected  by  increasing  slightly  the  modulus  of 
concrete  E  .  However,  this  small  difference  may  be 
merely  due  to  the  test  method  explained  earlier  that 
while  the  strains  are  computed  exactly  at  the  joint  with 
the  load  placed  directly  next  to  the  joint,  the  strain 
gages  were  actually  placed  at  a  small  distance  away  from 
the  joint  and  the  load  is  then  placed  behind  the  gages. 

e.  ORDL  Test  Track  "A”. 

(1)  Three  different  airfield  pavements  were  tested  at 
the  ORDL  test  track.  Pavement  and  joint  information 
can  be  found  in  Table  5.  The  measured  strains  shown 
in  entry  5-a  are  the  average  strains  for  all  mea¬ 
surements.  The  strains  computed  for  the  three  test 
pavements  are  shorn  separately  in  entries  5-b  to 
5-d.  The  average  computed  strains  and  stress  trans¬ 
fer  are  listed  in  entry  5-e.  While  the  computed 
strain  values  vary  among  the  three  different  pave¬ 
ments,  the  computed  stress  transfers  are  nearly  the 
same.  This  is  reasonable  because  the  test  pavements 
were  designed  and  constructed  in  accordance  with  the 
criteria,  so  the  computed  stress  transfers  across 
the  joints  for  the  three  pavements  are  nearly  iden¬ 
tical.  Also,  under  a  given  loading  condition,  the 
strain  decreases  as  the  thickness  of  the  concrete 
pavement  increases.  This  explains  why  the  computed 
strain  varies  among  the  three  pavements. 

(2)  The  average  strains  and  stress  transfer  shown  in 
entry  5-e  indicate  that  the  computed  total  'strain 
and  stress  transfer  are  considerably  larger  than  the 
measured  values.  The  strain  at  the  loaded  side  of 
the  joint  is  slightly  less  than  the  measured,  but 
the  strain  at  the  unloaded  side  is  much  greater  than 
the  measured,  resulting  in  a  much  greater  computed 
stress  transfer.  With  such  strong  dowel  bars  in  the 
joints,  i.e.,  2,  3,  and  k  in.  in  diameter,  it  is 
hard  to  believe  that  the  average  transfer  is  only 
28.1  percent.  The  only  explanation  for  the  discrep¬ 
ancies  is  that  the  pavements  at  the  ORDL  test  tracks 
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used  the  "remove  and  replace"  construction  method  of 
installing  the  doweled  joints.  It  was  possible  that 
the  joints  were  very  loose. 

(3)  To  take  into  account  the  looseness  of  the  doweled 
joint  in  the  computation,  the  modulus  of  dowel  sup¬ 
port  K  can  be  reduced.  In  doing  so,  as  is  the 
case  in  the  computed  values  shown  in  entry  1-d,  the 
computed  strain  at  the  loaded  side  will  increase 
slightly,  and  the  strain  at  the  unloaded  side  will 
decrease  considerably.  The  combined  effect  will  re¬ 
duce  both  the  computed  total  strain  and  stress 
transfer.  However,  computation  was  not  made  for 
smaller  value  of  the  modulus  of  dowel  support  K  . 

Ellsworth  AFB.  South  Dakota.  The  values  shown  in  entry 
6-b  are  computed  with  E  ,  K  ,  and  d  equal  to 
6,000,000  psi,  1,500,000  pci,  and  1/32  in.,  respective¬ 
ly.  The  computed  strains  and  stress  transfer  are  much 
smaller  than  the  measured  values.  To  increase  the 
strains,  the  modulus  of  the  concrete  can  be  reduced; 
and  to  increase  the  stress  transfer,  the  modulus  of 
dowel  support  can  be  increased.  Entry  6-c  shows  the 
computed  values  for  E  and  K  equal  to  4,000,000  psi 
and  14,000,000  pci,  respectively.  The  computed  strains 
and  stress  transfer  become  quite  close  to  the  measured 
values . 

Beale  AFB,  California.  Entry  7-b  shows  the  values  com¬ 
puted  with  E  ,  fC  ,  and  d  equal  to  6,000,000  psi, 
1,500,000  pci,  and  1/32  in.,  respectively.  The  computed 
strains  and  stress  transfer  are  smaller  than  the  measured 
values.  The  differences  can  be  reduced  by  using  a 
greater  value  of  the  modulus  of  dowel  support  K  . 

March  AFB.  California.  Entry  8-b  shows  the  values  com¬ 
puted  with  E  ,  K  ,  and  d  equal  to  6,000,000  psi, 
1,500,000  pci,  and  1/32  in.,  respectively.  The  computed 
strains  and  stress  transfer  are  very  close  to  measured 
values . 


Dow  AFB,  Maine. 

(l)  The  measured  values  are  shown  in  entry  9-a.  The 
11.5  percent  stress  transfer  indicates  that  the  ef¬ 
ficiency  of  load  transfer  across  the  joint  in  the 
pavements  at  Dow  AFB  must  be  very  low.  Computation 
was  first  made  with  E  ,  K  ,  and  d  equal  to 
6,000,000  psi,  1,500,000  pci,  and  1/32  in.  As  ex¬ 
pected,  the  computed  strains  and  stress  transfer 
were  too  large.  Computation  was  made  by  reducing 
the  modulus  of  dowel  support  from  1,500,000  to 
300,000  psi.  The  stress  transfer  was  reduced  con¬ 
siderably,  but  considerable  difference  still  exists 
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in  the  strains.  The  computed  values  are  shown  in 
entry  9-c. 

(2)  When  the  dowel  bars  are  loose  in  the  joint,  the 
modulus  of  dowel  support  K  must  be  low.  It  is 
possible,  however,  that  the  K  value  is  low  only 
when  the  condition  of  looseness  prevails.  When  the 
looseness  between  the  dowel  bars  and  the  concrete 
is  eliminated  by  the  applying  load,  a  much  larger 
K  value  should  be  used.  Computation  was  made  for 
such  a  case  for  the  pavement  at  Dow  AFB.  A  very 
small  K  value  of  50,000  pci  was  first  used  in  the 
computation  when  the  deformation  of  the  concrete  is 
less  than  0.01  in.;  i.e.,  the  looseness  between  the 
concrete  and  the  dowel  bars  is  assumed  to  be  0.01 
in.  When  the  deformation  exceeds  this  value,  a  K 
value  of  1,500,000  pci  was  used.  The  computed  value 
based  on  such  a  combination  is  given  in  entry  9-d. 
Because  a  very  low  initial  K  value  was  used,  the 
computed  stress  transfer  became  very  small. 

95-  The  comparison  between  the  measured  and  computed  performance 
presented  in  Table  8  indicates  that  the  WESLIQID  computer  program  can 
yield  satisfactory  results  to  predict  joint  performance  for  a  rigid 
pavement.  Generally,  6,000,000  psi  and  1,500,000  pci  can  be  used  for 
the  modulus  of  the  concrete  E  and  the  modulus  of  the  dowel  support 
K  in  most  cases.  In  some  cases,  other  values  of  E  and  K  have  to 
be  used.  Unfortunately,  it  is  very  difficult,  if  not  impossible,  to 
evaluate  the  values  of  E  and  K  for  either  a  new  or  existing  pave¬ 
ment.  The  E  value  depends  on  the  physical  properties  at  any  given 
point  in  time  of  the  concrete,  and  the  modulus  dowel  support  K  depends 
on  the  construction  method,  the  concrete  properties,  and  the  condition 
of  the  concrete  surrounding  the  dowel  bars. 

Factors  affecting  the 

stress  transfer  across  a  joint  1 

96.  The  computed  results  presented  in  Table  8  indicate  that  the 
stress  (or  load)  transfer  across  a  joint  as  calculated  in  program  WESLIQID 
is  influenced  by  many  factors.  A  series  of  computations  were  made  for 
the  pavement  at  Lockbourne  AFB  in  Ohio  subject  to  the  same  loads  to  de¬ 
termine  the  stress  transfer  across  the  joint.  The  variables  were  the 
modulus  of  dowel  support  K  ,  the  modulus  of  concrete  E  ,  the  joint 
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(or  crack)  spacing  d  ,  pavement  thickness  t  ,  modulus  of  subgrade 
reaction  k  ,  and  dowel  bar  diameter  b  .  The  computed  results  are  pre¬ 
sented  in  Figure  1^. 

97.  Figure  lUa  shows  that  the  stress  transfer  across  a  joint  is 
greatly  influenced  by  the  modulus  of  dowel  support  K  ,  particularly  in 
the  range  from  0  to  1,500,000  pci.  When  the  K  value  is  increased, 
i.e.,  the  dowel  bars  have  a  tighter  fit  in  the  joint,  more  load  can  be 
transferred  from  one  slab  to  the  other  across  the  joint. 

98.  Figure  l*+b  shows  that  within  the  possible  range  of  the  mod¬ 
ulus  of  concrete  E  ,  the  stress  transfer  is  influenced  by  E  to  some 
extent.  It  was  discussed  earlier  that  increasing  E  would  greatly  re¬ 
duce  the  strains  in  the  concrete  and  that  since  the  decrease  is  slightly 
more  in  the  loaded  side  than  the  unloaded  side  of  the  slab,  the  percent 
stress  transfer  is  consequently  reduced  when  E  is  increased  as  shown 
in  Figure  lUb. 

99.  Figure  lUc  shows  the  effect  of  joint  (or  crack)  spacing  d 
on  the  load  transfer  mechanism  across  the  joint.  It  is  seen  that  the 
stress  transfer  is  changed  merely  2  percent  when  the  spacing  d  is 
increased  from  1/32  to  1/2  in.  This  is  reasonable  because  the  stress 
transfer  across  a  joint  is  principally  due  to  the  shear  force  and  not 
due  to  the  moment  transfer  of  the  steel  bars. 

100.  Figure  lUd  shows  the  effect  of  pavement  thickness  t  on  the 
capability  of  stress  transfer  across  the  joint.  The  comparison  is  based 
on  the  same  loading  condition  and  other  variables.  When  the  pavement 
thickness  is  increased,  the  percent  stress  transfer  is  reduced  and  vice 
versa.  Evidently  under  a  given  load  when  the  pavement  thickness  is 
increased,  not  only  are  the  stresses  and  the  deflections  reduced,  but 
also  disproportionately  lesser  stresses  are  transferred  to  the  other 
slab.  If  it  is  desired  to  design  rigid  pavements  for  different  perform¬ 
ance  levels  subjected  to  the  same  aircraft  load,  it  is  expected  that 
thinner  pavements  will  have  greater  stress  transfer  capacity  across  the 
joint,  although  the  stresses  in  the  thinner  pavement  are  greater  and 
the  joint  condition  is  the  same  for  the  pavements. 

101.  It  is  interesting  to  review  the  test  pavements  in  the  ORDL 
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Figure  1^.  Variations  of  stress  (or  load)  transfer  across  a  joint  (Sheet  1  of  2) 
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test  track  in  connection  with  the  stress  transfer  capability.  The  three 
test  pavements  listed  in  Table  5  have  thicknesses  of  2b ,  28,  and  32  in., 
and  their  corresponding  dowel  diameters  and  spacings  are  2,  3,  and  b  in. 
and  IT,  IT,  and  2T  in.,  respectively.  The  numbers  were  determined 
based  on  the  Corps  of  Engineers  design  criteria  for  rigid  pavement.  In 
other  words,  the  stress  transfers  for  the  three  test  pavements  were  ex¬ 
pected  to  be  the  same.  The  measured  stress  transfers  for  the  three 
pavements  were  all  close,  and  the  computed  stress  transfers  shown  in 
Table  8  were  nearly  identical,  even  though  the  measured  and  computed 
values  for  each  pavement  do  not  agree  with  each  other. 

102.  Figure  l4e  shows  the  effect  of  the  modulus  of  subgrade  re¬ 
action  k  on  the  stress  transfer  capability  of  the  joint.  As  the  k 
value  is  increased,  the  stress  transfer  is  decreased.  Similar  to  the 
discussion  about  the  relationship  between  the  pavement  thickness  t 
and  the  stress  transfer  shown  in  Figure  l4d,  when  the  strength  of  the 
subgrade  is  increased,  the  stresses  and  deflections  are  reduced,  but 
disproportionately  lesser  stresses  are  transferred  to  the  other  slab. 
Therefore,  if  the  subgrade  soil  of  a  rigid  pavement  is  stabilized  and 
other  conditions  are  kept  identical,  the  stress  transfer  across  the 
joint  can  be  expected  to  reduce  slightly,  although  the  stresses  in  the 
pavement  with  stabilized  subgrade  are  smaller. 

103.  Figure  l4f  shows  the  effect  of  dowel  bar  diameter  b  on  the 
stress  transfer  capability.  The  computations  were  made  for  different 
diameters  of  dowel  bars  at  a  constant  6-in.  dowel  spacing.  As  was  ex¬ 
pected,  the  amount  of  stress  transfer  is  dominated  by  the  number  and 
size  of  dowel  bars  used  in  the  joint.  It  should  be  pointed  out  that 
according  to  the  definition  of  percent  stress  transfer  used  in  this 
report,  a  50  percent  stress  transfer  is  a  total  load  transfer,  i.e., 
the  stresses  on  both  sides  of  the  joint  are  identical.  For  a  4-in. 
dowel  bar  spaced  6  in.  apart,  the  stress  transfer  is  still  less  than 

a  total  transfer,  thus  indicating  that  the  dowel  bar  is  not  a  very 
efficient  load  transfer  device.  This  aspect  will  be  discussed  later  in 
the  report. 
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10U.  The  computed  results  presented  in  this  paper  indicate  that 
the  stress  (or  load)  transfer  across  a  joint  in  a  rigid  pavement  is 
affected  by  many  factors.  These  factors,  such  as  bar  diameter  and 
spacing,  subgrade  modulus,  and  pavement  thickness,  can  be  readily 
measured  and  thus  appropriately  accounted  for.  However,  the  factors, 
such  as  modulus  of  concrete  E  and  modulus  of  dowel  support  K  ,  are 
hard  to  measure  and  estimate,  but  their  effects  on  the  stress  transfer 
are  very  large.  It  is  also  known  that  the  value  of  K  in  a  rigid 
pavement  can  decrease  with  time  because  of  the  repetitive  action  of 
aircraft  loads. 
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PART  IV:  PRESENTATION  OF  NUMERICAL  RESULTS 
FOR  THE  WESLAYER  PROGRAM 


Introduction 


105.  In  using  the  WESLAYER  computer  program  to  analyze  field  mea¬ 
surements  of  test  sections,  the  primary  difficulty  lies  in  the  selection 
of  proper  modulus  values  F,  and  Poisson's  ratio  v  to  represent  the 
material  properties  of  the  layered  elastic  subgrade  because  most  test 
sections  measured  k  (modulus  of  subgrade  reaction)  values  for  sub¬ 
grade  soils.  Efforts  were  made  to  establish  correlations  between  k 
and  E  values  for  a  given  value  of  v  for  the  case  of  a  single  layer 
under  the  pavement.  This  was  done  by  computing  the  stresses  and  deflec¬ 
tions  for  a  given  pavement  with  a  range  of  subgrade  k  values  (using 
the  WEULIQID  program)  and  for  the  same  pavement  with  a  r^nge  of  subgrade 
E  values  (using  the  WEU LAYER  program  and  assuming  the  i  ■  ograde  is 
homogeneous  and  v  =  0.fi).  The  corresponding  E  and  k  values  were 
determined  by  matching  the  maximum  (or  near  maximum  as  will  be  explained 
later)  stress  or  deflection  at  a  critical  nodal  point  in  two  pavements 
that  have  the  same  loading  and  geometrical  conditions  but  have  different 
subgrade  conditions.  The  reason  for  selecting  maximum  values  is  that 
the  shapes  of  the  stress  and  deflection  basins  under  the  liquid  founda¬ 
tion  (represented  by  the  WESLIQID  program)  and  the  elastic  foundation 
(represented  by  the  WESLAYER  program)  are  different.  In  the  computation 
the  elastic  modulus  of  the  concrete  E  was  assumed  to  be  6,000,000  psi. 

Stress  and  Deflection  Basins 


106.  The  computer!  results  indicate  that  the  correlations  between 
the  k  and  K  values  are  not  unique.  The  correlations  are  presented 
in  Fii"ure  15  for  three  ,)oint  conditions  for  stress-based  computations.* 


*  The  equation  relating  k  and  E  values  developed  by  Vesic  and 
Saxena  (19T0)  was  compared  with  the  relations  shown  in  Figure  15  and 
was  found  unsatisfactory  in  some  cases. 
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MODULUS  OF  SUBGRADE  REACTION  k,  PCI 


Figure  15-  Correlations  between  subgrade  E  and  k  for 
the  Lockbourne  AFB  test  pavement;  the  subgrade  is  assumed 
to  be  one  layer  and  homogeneous  and  the  Poisson's  ratio 

v  assumed  to  be  0.4 


The  pavement  and  loading  conditions  are  the  same  as  those  used  in  the 
Lockbourne  AFB  test  pavement  (see  Table  5).  Figure  15  shows  that  the 
differences  among  the  three  different  joint  conditions  are  large,  espe¬ 
cially  for  the  dowel  bar  case.  Computations  were  also  made  for  other 
pavements.  It  was  found  that  the  correlations  were  different  for  dif¬ 
ferent  pavements  and  different  loading  positions. 

107.  As  explained  earlier,  the  correlations  of  E  and  v  ,  shown 
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in  Figure  15,  were  established  by  matching  the  maximum  stress  or  deflec¬ 
tions  at  a  certain  nodal  point  in  the  pavement.  The  reason  for  select¬ 
ing  the  maximum  (or  near  maximum)  values  is  that  the  shapes  of  the 
stress  and  deflection  basins  under  the  liquid  foundation  (represented 
by  the  WESLIQID  program)  and  the  elastic  foundation  (represented  by  the 
WESLAYER  program)  are  different  due  to  the  nature  of  the  assumption. 

In  a  liquid  subgrade,  the  deflection  at  a  node  depends  solely  on  the 
modulus  of  subgrade  reaction  k  at  the  node  and  not  elsewhere.  In  an 
elastic  subgrade,  however,  the  deflection  at  a  node  depends  not  only  on 
the  elastic  modulus  E  of  the  subgrade,  but  also  on  the  deflections  at 
other  nodes.  Consequently,  the  deflection  basin  should  be  steeper  in 
the  liquid  subgrade  than  in  the  elastic  subgrade. 

108.  The  stress  and  deflection  basins  for  the  Lockbourne  AFB  test 
pavement  for  two  different  joint  conditions  are  presented  in  Figure  16. 
The  stresses  (a  )  and  deflections  are  normalized  with  respect  to  the 
values  at  node  31  (see  Part  V)  at  the  joint.  These  basins  show  the 
cross  sections  along  the  X-axis  for  nodes  1,  6,  11,  l6,  21,  26,  and  31. 
The  computations  were  based  on  the  assumption  that  the  subgrade  elastic 
modulus  E  was  equal,  to  10,000  psi .  The  corresponding  values  of  modu¬ 
lus  of  subgrade  reaction  k  used  in  the  computations  are  determined 
from  Figure  15  for  different  cases. 

109-  Figures  l6a  and  l6b  show  the  normalized  deflection  basins 
for  the  conditions  of  100  and  50  percent  shear  transfers  across  the 
joint,  respectively.  It  can  be  seen  that  the  basins  are  much  flatter 
when  the  subgrade  is  elastic  than  when  it  is  liquid  and  when  the  joint 
has  better  load  transfer  capability. 

110.  Figures  l6c  and  l6d  show  the  shapes  of  the  stress  basins  for 
the  two  joint  conditions.  The  curvatures  of  the  stress  basins  for  the 
elastic  and  liquid  subgrades  are  nearly  the  same.  It  is  interesting  to 
note  the  hump  near  the  wheel  load  at  node  21  in  the  stress  basin  for  the 
liquid  subgrade  shown  in  Figure  l6c .  It  is  believed  that  the  hump 
(sharp  change  of  stress  magnitude)  is  attributable  to  the  assumption  of 
liquid  foundation;  i.e.,  the  large  deflection  under  the  applied  wheel 
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load  is  primarily  supported  by  the  spring  under  the  load  and  causes 
larger'  stresses  in  the  slab  near  the  load.  Tn  the  elastic  subgrade 
condition,  however,  the  wheel  load  is  supported  by  the  entire  elastic 
soil  and  the  distribution  of  slab  stress  around  the  load  is  more  uni¬ 
form  and  thus  the  hump  is  not  observed  in  the  elastic  subgrade  condi¬ 
tion.  Nevertheless,  the  basins  shown  in  Figure  l6  indicate  that  the  use 
of  elastic  subgrade  soil  (WESLAYER  program)  can  better  represent  the 
response  of  the  slab  to  the  load  than  does  the  liquid  subgrade  soil 
(WESLIQID  program). 

111.  The  question  may  arise  as  to  why  humps  shown  in  Figures  l6c 
and  l6d  are  not  observed  in  deflection  basins  in  Figures  l6a  and  l6b 
for  the  liquid  subgrade  case.  The  answer  may  be  that  although  the  sub¬ 
grade  springs  near  the  wheel  load  are  affected  much  more  than  those  else¬ 
where,  the  portion  of  the  slab  near  the  load  cannot  deflect  much  more 
than  elsewhere  because  of  the  rigidity  of  the  concrete  slab. 

112.  It  should  be  pointed  out  here  that  although  maximum  stress 
occurs  at  node  26,  the  equivalent  subgrade  modulus  was  determined  by 
matching  the  stresses  at  node  31. 

Comparisons  with  Strain  Measurements 
from  the  Corps  of  Engineers 

113.  The  test  pavements  shown  in  Table  5  were  analyzed  by  the 
WESLIQID  program  and  the  results  are  presented  earlier  in  Part  III  of 
this  report.  The  same  pavements  were  also  analyzed  by  the  WESLAYER 
program  and  the  results  are  presented  in  the  following  paragraphs.  In 
the  computations,  the  E  values  for  each  pavement  were  determined  from 
the  k  versus  E  relations  for  the  dowel  bar  case.  The  finite  element 
layout  is  shown  in  Part  V. 

11*4.  The  percentages  of  stress  transfer  across  the  joint  are 
computed  based  on  stress  o  at  nodes  31  and  36.  Table  9  shows  the 
computed  results  for  five  test  pavements.  All  of  the  pavements  were 
tested  under  the  B-U7  aircraft  loads.  For  completeness,  the  percentages 
of  stress  transfer  computed  for  liquid  subgrade  using  the  WESLIQID 
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program  are  also  included  in  the  table.  It  is  seen  that  the  computed 
values  are  very  close  to  the  measured  ones  and  are  also  close  to  those 
computed  by  the  WESLIQID  program. 

115.  It  should  be  pointed  out  that  the  correlations  shown  in 
Figure  15  are  established  only  for  the  Loekbourne  AFB  test  pavement 
and  are  strictly  not  suitable  for  determining  E  values  for  other  pave¬ 
ments.  However,  it  is  believed  that  the  difference  in  the  computed 
percentage  of  shear  transfer  may  not  be  large.  This  is  discussed  in  the 
following  paragraphs. 

Effect  of  Subgrade  Elastic  Modulus  E 
on  Stress  Transfer  Across  a  Joint 


116.  In  Part  II  of  this  report,  the  effects  of  a  number  of  impor¬ 
tant  parameters  on  stress  transfer  across  a  joint  were  studied  (by  the 
use  of  WESLIQID  program  for  pavements  on  liquid  subgrade);  the  results 
are  presented  in  Figure  ll).  The  parameters  studied  were  joint  spacing 
d  ,  modulus  of  concrete  E  ,  modulus  of  dowel  support  K  ,  dowel  bar 
diameter  b  ,  concrete  pavement  thickness  t  ,  and  modulus  of  subgrade 
reaction  k  .  It  is  believed  that  the  conclusions  derived  are  also 
applicable  to  pavements  on  elastic  subgrade  soils.  Efforts  were  thus 
made  to  investigate  only  the  effect  of  subgrade  elastic  modulus  E  on 
stress  transfer  across  a  joint.  The  results  computed  for  the  Loekbourne 
AFB  test  pavement  are  presented  in  Figure  17. 

117.  The  results  presented  in  Figure  17  for  elastic  subgrade  soils 
are  very  similar  to  those  presented  in  Figure  lUe  for  liquid  subgrade. 
The  percent  stress  transfer  decreases  as  the  subgrade  modulus  increases. 
This  is  reasonable  because  the  load  is  transferred  to  the  other  slab 
through  dowel  bars  when  the  loaded  slab  is  deflected.  In  the  extreme 
case  when  the  subgrade  modulus  is  increased  until  infinitely  large, 

the  loaded  slab  ceases  to  deflect  and  the  load  is  not  transferred  to 
the  other  slab  through  the  joint;  the  percent  of  stress  transfer  should 
thus  approach  zero. 

118.  Figure  17  shows  that  for  E  values  ranging  from  14,000  to 
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1-IN.  DOWEL  BARS,  15  IN.  C  TO  C 
PAVEMENT  THICKNESS  t  =  12  IN. 

SUBGRADE  IS  ASSUMED  TO  BE 
ONE  LAVER  AND  HOMOGENEOUS 


STRESS  TRANSFER  ACROSS  THE  JOINT  IN  PERCENT 

Figure  IT-  Relationship  between  stress  transfer  across 
a  joint  and  subgrade  elastic  modulus,  Lockbourne  AFB 

test  pavement 

200,000  psi,  the  percent  of  stress  transfer  varies  from  12.9  to  29-b. 
However,  the  subgrade  elastic  modulus  of  the  five  test  pavements  ana¬ 
lyzed  (Table  5)  ranges  only  from  7,^00  to  15,000  psi,  and  the  corre¬ 
sponding  variation  of  stress  transfer  is  only  3  percent,  shown  in 
Figure  17,  which  is  not  significant  at  all. 


L 


119-  It  should  be  pointed  out  that  if  the  option  of  percent  shear 
transfer  across  the  joint  is  used  rather  than  the  option  of  dowel  bars 
as  in  Figure  17,  the  percent  of  stress  transfer  across  the  joint  is  in¬ 
dependent  of  subgrade  elastic  modulus  E  (or  modulus  of  subgrade  reac¬ 
tion  k  ) . 


PART  V:  DESIGN  IMPLICATIONS 


120.  The  results  and  discussions  presented  in  this  part  of  the 
report  are  limited  to  the  use  of  WESLIQID  program,  i.e.,  pavements  on 
a  liquid  subgrade.  For  illustrative  and  comparative  purposes,  some 
example  problems  were  computed  using  either  heavy  loads,  or  thin  con¬ 
crete  slab,  or  weak  subgrade  soil,  or  extreme  temperature  differential; 
computed  stresses  were  thus  too  large  and  had  exceeded  the  elastic 
limit  of  the  concrete.  Also,  in  some  cases  the  load  intensities  were 
chosen  arbitrarily  and  had  resulted  in  very  large  stresses.  It  should 
be  noted,  however,  the  finite  element  method  employed  in  this  study  is 
limited  to  the  linear  elastic  theory;  i.e.,  the  stresses  and  displace¬ 
ments  are  linearly  proportional  to  the  load.  If  the  computed  stresses 
are  beyond  the  elastic  limit  of  the  concrete,  the  computed  stresses  can 
be  reduced  to  whatever  values  desired  when  the  load  is  proportionally 
reduced. 


Efficiency  of  Load  Transfer  by  Dowel  Bars 

121.  The  shear  forces  across  a  joint  in  a  rigid  pavement  are 
transferred  by  dowel  bars  and  concrete  interlock.  It  is  difficult  to 
evaluate  the  amount  of  shear  force  transferred  by  the  concrete  inter¬ 
lock,  but  the  amount  of  shear  force  transferred  by  the  dowel  bars  can 

be  evaluated  using  the  WESLIQID  program.  Using  the  pavement  and  loading 
condition  at  Lockbourne  AFB  (Table  5),  stresses  in  the  pavement  are  com¬ 
puted  for  different  sizes  and  spacings  of  dowel  bars,  assuming  either 
100  percent  or  zero  percent  moment  transfer  across  the  joint.  The  com¬ 
puted  stresses  at  selected  nodes  are  presented  in  Table  10  for  the  con¬ 
dition  of  100  percent  moment  transfer.  The  finite  element  layout  show¬ 
ing  the  nodal  numbers  is  presented  in  Figure  18.  Special  attention 
should  be  given  to  nodes  31  and  36  where  maximum  stresses  and  deflections 
occur.  Node  31  is  in  the  loaded  slab  and  node  36  is  the  corresponding 
node  in  the  unloaded  slab  across  the  joint. 

122.  When  the  efficiency  of  shear  and  moment  transfer  across  a 
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Figure  18.  Finite  element  layout  showing  nodal  numbers 
for  computations  presented  in  Table  10 


joint  are  both  100  percent,  the  joint  actually  does  not  exist.  The 

stresses  and  deflections  computed  in  such  a  case  are  tabulated  in  entry 

A  of  Table  10.  It  is  seen  that  at  nodes  31  and  36,  the  stresses  and 

deflections  are  practically  the  same.  At  nodes  6,  11,  l6,  21,  26,  31, 

36,  1+1,  1+6,  and  51,  the  shear  stresses  are  zeros  because  these  nodes 

lie  on  the  axis  of  symmetry.  Edge  stress  under  the  same  loading  was 

computed  but  is  not  presented  in  Table  10.  The  stresses  were  computed 

by  removing  slab  2  in  Figure  18  and  thus  treating  the  joint  as  a  free 

edge.  The  maximum  stress  was  1169-2  psi;  that  is,  more  than  twice  the 

maximum  stress  computed  at  node  31,  indicating  that  in  a  rigid  pavement 

edge  stress  is  most  critical.  It  is  to  be  note\l  that  the  sum  of  the 

stress,  either  0  or  a  at  nodes  31  and  36  are  nearly  equal  to  the 
y  max  . 

maximum  edge  stress. 

123.  Entry  B  presents  the  stresses  when  the  joint  is  connected  by 
1-in.  dowel  bars  spaced  15  in.  center  to  center.  It  is  seen  that  the 
presence  of  a  joint  in  the  pavement  does  not  change  very  much  the  stress 


71 


c»x  in  the  direction  perpendicular  to  the  joint,  but  changes  consider¬ 
ably  the  stress  o  near  the  load  in  the  direction  parallel  to  the 
joint.  The  reason  that  the  stresses  at  both  sides  of  the  joint 

are  equal  and  are  close  to  those  values  computed  before  the  joint  was 
created  is  due  to  the  assumption  of  100  percent  moment  transfer  across 
the  joint.  It  is  noted  that  the  stress  increases  considerably 

near  the  load  in  the  loaded  slab  but  decreases  considerably  in  the 
unloaded  side.  At  nodes  31  and  36  across  the  joint,  the  stress  (or 
load)  transfer  computed  from  a  has  a  value  of  265.6/(265.6  +  75^.3)  = 
26  percent  ,  and  the  shear  transfer  computed  from  the  deflections  is 
0. 05^/0. 065  =  83  percent  .  It  should  be  noted  that  the  shear  forces 
transferred  by  the  concrete  interlock  are  neglected  in  the  computation; 
otherwise,  the  stress  (or  load)  and  shear  transfers  may  be  greater. 

12U.  Comparing  the  values  between  entries  A  and  B,  deflections 
are  increased  in  the  loaded  slab ,  but  are  decreased  in  the  unloaded 
slab  when  a  (doweled)  joint  is  present.  Increasing  pavement  deflection 
means  greater  subgrade  stresses. 


125.  It  is  rather  interesting  to  note  that  when  a  joint  is  con¬ 
nected  by  dowel  bars,  the  sum  of  the  stresses  at  both  sides  of  the  joint 
is  approximately  equal  to  the  free  edge  stress.  At  nodes  31  and  36  in 
entry  B,  the  sum  of  o  is  (75^.3  +  265.6)  =  1019-9  psi,  and  the  sum 

v 

of  the  maximum  principal  stress  o  is  (75^.3  +  379.2)  =  1133.8  psi, 

max 

which  is  very  close  to  the  ll69.2-psi  free  edge  stress  computed  at 

node  31.  Note  at  the  free  edge  (node  31 ),  the  shear  stress  t  and 

xy 

the  stress  perpendicular  to  the  edge  are  all  zero;  the  maximum 

principal  stress  o  equals  the  stress  o 
max  y 

126.  Entry  C  presents  the  computed  values  for  the  case  of  2-in. 
dowel  bars  spaced  6  in.  center  to  center.  As  the  reinforcement  across 


the  joint  is  increased,  while  the  stresses  still  remain  nearly 

the  same  due  to  the  assumption  of  100  percent  moment  transfer,  the  de¬ 


flections  and  stresses  are  reduced  in  the  loaded  slab  but  are  in¬ 

creased  in  the  unloaded  one.  The  stress  (o  )  and  shear  (deflection) 
transfers  at  nodes  31  and  36  across  the  Joint  are  increased  from  26  to 
39  percent  and  from  83  to  97  percent,  respectively.  At  other  nodes  the 


72 


stresses  o 


and  deflections  are  reduced  in  the  loaded  slab  and  are 


max 

increased  in  the  unloaded  slab  when  the  number  and  size  of  the  dowel 
bars  are  increased.  Note  chat  in  this  case  the  stress  o  at  node  31 

y 

(6l8.3  psi)  is  considerably  larger  than  the  corresponding  value  shown 
in  entry  A  (505.6  psi)  where  100  percent  shear  and  moment  transfers  are 
assumed,  indicating  that  the  2-in.  dowel  bars  spaced  6  in.  apart  at  the 
.joint  still  cannot  transfer  the  full  stresses  from  the  loaded  slab  to 
the  unloaded  one;  i.e.,  the  stresses  on  both  sides  of  the  joint  are 
still  not  equal.  The  computed  values  shown  in  entry  C  also  demonstrate 
that  deflection  across  a  .joint  may  not  be  a  very  descriptive  parameter 
for  the  .joint  performance.  At  nodes  31  and  36 ,  the  shear  transfer  com¬ 
puted  based  on  deflections  is  97  percent,  but  the  stress  transfer  com¬ 
puted  based  on  stress  is  only  39  percent,  although  the  method  of 

calculating  the  percent  transfer  is  different  for  these  two  parameters 
(the  numerical  computations  for  these  values  can  be  found  in  paragraph 
123). 

127.  Similar  to  the  results  presented  in  entry  B,  entry  C  also 
shows  that  the  sum  of  the  stresses  at  both  sides  of  a  doweled  joint  is 
nearly  equal  to  the  free  edge  stress. 

128.  Entry  D  shows  the  computed  values  for  the  case  of  extremely 
heavy  reinforcement  across  the  joint.  The  dowel  bars  are  8  in.  in 
diameter  and  are  spaced  9  in.  center  to  center;  i.e.,  the  spacing  be¬ 
tween  the  dowels  is  only  1  in.  The  stresses  and  deflections  are  nearly 
identical  to  those  shown  in  entry  A  in  which  100  percent  of  shea  and 
moment  are  transferred.  This  is  reasonable  because  with  the  extremely 
heavy  reinforcement  across  the  joint  and  with  the  assumption  of  a  100 
percent  moment  transfer,  the  joint  becomes  so  rigid  that  the  slabs  act 
as  if  the  joint  does  not  exist. 

129.  Table  10  presents  the  computed  results  under  the  assumption 
that  the  joint  has  100  percent  moment  transfer.  Similar  computations 
were  also  made  for  the  assumption  that  the  joint  is  not  capable  of 
transferring  moments  at  all  (Table  11 ).  The  purposes  of  the  computa¬ 
tions  were  twofold.  The  first  was  to  check  the  efficiency  of  dowel  bars 
in  transferring  the  stresses  across  a  joint  in  which  the  efficiency  of 


moment  transfer  is  not  assumed.  The  second  was  to  compare  the  stress 
and  deflection  distributions  in  the  pavement  as  a  function  of  the  effi¬ 
ciency  of  moment  transfer  across  the  joint. 


130.  When  the  computed  values  in  entry  A  of  Table  11  are  compared 
with  those  in  entry  A  of  Table  10,  it  is  seen  that  the  stresses  ar.d  de¬ 
flections  change  drastically  along  the  joint  when  the  efficiency  of  mo¬ 
ment  transfer  is  reduced  to  zero.  (Note  that  100  percent  shear  transfer 
is  not  changed.)  The  stresses  (perpendicular  to  the  joint)  at  no¬ 

dal  points  along  the  joint  are  all  diminished,  and  the  deflections  are 

increased  greatly.  Stress  a  and  o  ,  however,  are  changed  insig- 

y  max 

nificantly.  Greater  deflections  mean  greater  stresses  in  the  subgrade. 
The  significance  of  this  will  be  discussed  later.  At  nodes  away  from 
the  edge,  however,  the  deflections  are  less  and  the  stresses  are  slightly 
less  in  the  loaded  slab  and  are  generally  slightly  more  in  the  unloaded 
slab,  when  the  efficiency  of  moment  transfer  is  reduced  to  zero.  Because 
the  shear  transfer  is  assumed  to  be  100  percent  efficient,  the  stress 
transfer  at  nodes  31  and  36  across  the  joint  has  a  full  transfer,  i.e., 

50  percent. 


131.  Entry  B  shows  the  computed  values  for  1-in.  bars  spaced  15 
in.  apart,  as  compared  with  the  100  percent  shear  transfer  in  entry  A. 
Because  of  the  weaker  shear  transfer  across  the  joint,  the  deflections 
at  nodes  along  the  joint  increase  in  the  loaded  slab  and  decrease  in 
the  unloaded  slab.  This  is  also  true  for  stress  o  at  nodes  under  the 

y 

load  that  result  in  a  stress  (or  load)  transfer  of  289.7/(289.7  + 

812.7)  =  26  percent  and  a  shear  transfer  of  0.058/0.071  =  82  percent  . 
These  values  are  identical  to  those  computed  and  presented  in  entry  B 
of  Table  10  (or  see  paragraph  123)  in  which  100  percent  moment  transfer 
is  assumed.  In  other  words,  the  efficiency  of  moment  transfer  across 
a  joint  does  not  affect  the  values  of  the  stress  (or  load)  and  shear 
transfers,  although  the  values  of  deflections  and  stresses  a  along 

■  —  —  .  ..  1  .... - -  ■-  - -  —  ^  —  —  ■  ■  —  — ,  .  JL m 

the  joint  are  generally  greater  in  the  case  when  zero  percent  moment 

transfer  is  assumed. 

132.  At  nodes  away  from  the  joint,  the  stresses  and  deflections 


lh 


are  also  increased  slightly  in  the  loaded  slab  and  decreased  slightly  in 
the  unloaded  one  when  the  efficiency  of  shear  transfer  across  the  joint 
is  reduced. 


133*  Entry  C  shows  the  computed  values  for  the  case  when  2-in. 
dowel  bars  spaced  6  in.  apart  are  installed  at  the  joint.  Similar  to 


the  conclusions  derived  from  the  computed  values  shown  in  entry  C  of 
Table  10  (100  percent  moment  transfer  across  the  joint),  the  deflections 
and  stresses  along  the  joint  are  decreased  in  the  loaded  slab  and 

increased  in  the  unloaded  slab  when  the  number  and  size  of  dowel  bars 


are  increased.  The  computed  stress  (or  load)  and  shear  transfers  across 
the  joint  at  nodes  31  and  36  are  E4l.2/(UHl.2  +  658.3)  =  percent 
and  0.063/0.065  =  97  percent  ,  respectively.  Again,  these  values  are 
identical  with  those  computed  in  entry  C  of  Table  10  in  which  the  same 
number  of  dowel  bars  is  used  but  the  efficiency  of  moment  transfer  is 
assumed  to  be  100  percent. 

131*.  At  nodes  away  from  the  joint,  the  deflections  and  stresses 
decrease  in  the  loaded  slab,  but  increase  in  the  unloaded  slab,  as  the 
number  of  dowels  across  the  joint  increases. 

135-  Several  conclusions  can  be  drawn  from  the  computations  pre¬ 
sented  in  Tables  10  and  11.  They  are  discussed  as  follows: 

a.  While  the  conventional  dowel  bars  cannot  transfer  moment 
at  all,  they  are  not  a  very  effective  device  for  trans¬ 
ferring  stresses  (or  loads)  across  a  joint  either.  For 
instance,  for  1-in.  dowel  bars  spaced  15  in.  apart,  the 
stress  transfer  is  only  26  percent;  i.e.,  the  stress 

o  in  the  loaded  slab  is  about  three  times  as  much  as 
t^e  stress  in  the  unloaded  slab  across  the  joint.  The 
stress  transfer  is  increased  to  39  percent  for  2-in. 
dowel  bars  spaced  6  in.  apart;  i.e.,  the  stress  o  in 
the  loaded  slab  is  one  and  one-half  times  as  much  Ss 
the  stress  in  the  unloaded  slab.  The  computed  values  are 
the  same  for  either  100  or  zero  percent  efficiency  of 
moment  transfer  across  the  joint.  As  was  previously 
stated,  100  percent  moment  transfer  is  defined  by  equal 
rotations  at  nodal  points  at  both  sides  of  the  joint 
with  the  moments  computed  accordingly.  Zero  percent 
moment  transfer  is  defined  to  be  that  the  moments  at 
nodal  points  along  the  joint  are  all  zero  while  the  ro¬ 
tations  are  not  required  to  be  zero. 

b.  For  a  given  number  and  size  of  dowel  bars  at  the  joint 


and  when  the  efficiency  of  moment  transfer  is  reduced 
to  zero,  i.e.,  the  opening  of  the  joint  becomes  visi¬ 
ble,  the  computed  results  show  that  (l)  the  deflections 
along  the  joint  increase  and  the  stresses  cx  perpen¬ 
dicular  to  the  joint  diminish,  (2)  the  stress  ay  paral¬ 
lel  to  the  joint  and  the  deflection  near  the  load  increase 
slightly  in  both  the  loaded  and  the  unloaded  slab,  but 
the  values  of  the  load  and  shear  transfers  are  generally 
not  changed,  and  (3)  at  nodal  points  away  from  the  joint, 
the  deflections  are  reduced  and  the  stresses  are  reduced 
slightly  in  the  loaded  slab  but  are  generally  increased 
slightly  in  the  unloaded  slab.  Greater  pavement  deflec¬ 
tion  implies  greater  subgrade  stress  and  vice  versa. 

£.  While  the  vertical  deflections  are  easy  to  measure  in  a 
rigid  pavement,  one  should  be  cautious  in  using  deflec¬ 
tions  to  estimate  the  efficiency  of  stress  transfer 
across  a  joint  because  deflections  do  not  change  as 
much  as  the  stress  a y  along  the  joint.  For  instance, 
entry  C  of  Tables  10  and  11  shows  that  for  2-in.  dowel 
bars  spaced  6  in.  apart,  the  stress  at  the  node  at 

the  loaded  slab  is  one  and  one-half  times  greater  than 
that  at  the  corresponding  node  at  the  unloaded  slab, 
but  the  difference  in  the  deflections  between  the  two 
nodes  is  only  3  percent,  which  can  be  difficult  to 
measure. 


Effect  of  Joint  Conditions  on  Stress  j  and  Deflections  for 
Center  and  Joint  Loading  Conditions 


136.  It  is  known  that  the  condition  of  stress  transfer  at  the 
joint  can  affect  the  stress  distributions  in  the  pavement  and  therefore 
influence  the  thickness  design.  The  computed  results  by  WESLIQID  found 
that  this  is  true  only  if  the  load  is  placed  next  to  the  joint.  When 
the  load  is  placed  at  the  pavement's  center,  the  joint  condition  has  a 
minimal  effect  on  the  stress  distribution  in  the  pavement.  Figure  19 
shows  the  finite  element  layouts  of  a  two-slab  system  for  two  loading 
positions.  One  loading  position  is  the  internal  load  where  the  load  is 
placed  at  the  center  of  the  slab,  and  the  other  position  is  the  joint 
load  where  the  load  is  placed  next  to  and  at  the  center  of  the  joint. 
Because  of  symmetry,  only  half  01  +he  slabs  are  used  for  the  computations 
as  shown  in  Figure  19.  Note  that  a  shorter  length  of  unloaded  slab  is 
used,  which  is  permissible.  In  the  computation,  the  slab  thickness. 
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Figure  19.  Finite  element  layouts  showing  nodal  numbers 
for  computations  presented  in  Table  12  for  a  center  load 
and  in  Table  13  for  a  load  next  to  a  joint 

elastic  modulus,  and  Poisson's  ratio  of  the  concrete,  and  the  modulus 
of  subgrade  reaction  are  8  in.,  6,000,000  psi,  0.2  and  100  pci,  respec¬ 
tively.  The  uniformly  applied  load  has  a  pressure  of  125  psi. 

137.  Tables  12  and  13  show  the  computed  stresses  and  defiections 
at  selected  nodal  points  for  the  center  and  the  joint  loading  conditions, 
respectively.  Assuming  the  efficiency  of  moment  transfer  to  be  zero, 
the  computations  were  made  for  three  different  efficiencies  of  shear 
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transfer.  The  definition  of  100  percent  shear  transfer  is  that  the  de¬ 
flections  at  both  sides  of  the  joint  are  equal,  and  for  zero  percent  of 
shear  and  moment  transfers,  the  load  transfer  across  the  joint  is  zero. 

138.  Table  12  shows  that  when  the  load  is  far  away  from  the  joint, 
the  condition  of  the  joint  has  no  effect  at  all  on  the  stresses  and  de¬ 
flections  in  the  slab,  except  at  nodes  near  the  joint  where  the  stresses 
and  deflections  increase  as  the  efficiency  of  the  joint  decreases.  Since 
the  stresses  and  deflections  are  so  small  near  the  joint  when  the  load  is 
far  away  from  the  joint,  the  change  of  stresses  and  deflections  has  no 
effect  on  the  overall  pavement  design.  Unlike  the  results  presented  in 
Table  12,  results  presented  in  Table  13  show  that  when  the  load  is  placed 
next  to  the  joint,  the  stresses  and  deflections  at  nodes  near  the  joint 
are  affected  significantly  by  the  joint  efficiency.  For  instance,  at 
node  57  where  maximum  stresses  and  deflections  occur,  the  stress  a  is 

y 

23^0.9  psi  when  load  transfer  at  the  joint  is  zero,  i.e.,  the  edge  load 
case.  The  stress  a^r  is  reduced  to  II60.8  psi  when  the  efficiency  of 
shear  transfer  is  increased  to  100  percent,  a  50  percent  reduction  from 
the  edge  load  condition.  Since  in  an  actual  airfield  pavement  with  1- 
to  1-1/2-in.  dowel  bars  spaced  12  in.  apart  at  the  joint,  the  measured 
shear  transfer  could  be  about  90  percent,  the  reduction  of  stress  in  the 
pavement  should  thus  be  very  significant,  possibly  a  reduction  of  U5  per¬ 
cent  from  the  edge  load.  The  results  presented  in  Table  13  indicate  that 
the  efficiency  of  stress  transfer  at  the  joint  affects  the  pavement  de¬ 
sign  when  the  load  is  placed  near  the  joint. 

Effect  of  Loading  Position  on  Stresses  and 
Deflections  in  Jointed  Pavements 


139«  The  results  presented  in  Table  13  indicate  that  when  the 
load  is  placed  close  to  the  joint,  stress  can  be  reduced  with  an  effi¬ 
cient  stress  transfer  device  at  the  joint.  A  question  arises  concerning 
whether  the  pavement  thickness  can  be  reduced  because  of  the  existence 


of  efficient  dowel  bars  in  the  transverse  joint.  The  WESLIQID  program 
was  used  to  analyze  such  conditions,  and  the  results  are  presented  in 


Tables  lit  and  15  for  edge  and  center  loading,  respectively. 
Edge  loads 


lUO.  Figure  20  shows  the  four  edge  loading  positions  in  the  pave¬ 
ment.  For  clarity,  the  finite  element  layout  is  shown  only  in  Figure 
20c.  The  stresses  are  computed  at  the  nodes  where  the  nodal  numbers 
are  shown.  The  loading  used  in  the  computation  consists  of  a  twin- 
tandem  truck  load,  with  two  concentrated  loads  representing  one-wheel 
load.  The  use  of  concentrated  loads  rather  than  uniformly  applied  loads 
is  for  the  convenience  of  preparing  input  data  in  the  finite  element 
programs.  The  computation  was  first  made  under  the  assumption  that  the 
transfer  joint  does  not  exist,  as  shown  in  Figure  20a.  The  computations 
were  then  followed  by  the  cases  that  a  transverse  joint  is  present  when 
(a)  the  twin-tandem  load  is  next  to  the  joint  (Figure  20b),  (b)  the 
twin-tandem  load  is  b  ft  away  from  the  joint  (Figure  20c),  (c)  the  front 
axle  load  is  next  to  the  joint  (Figure  20d)  ,  and  (d)  the  rear  axle  load 
is  U  ft  away  from  the  joint  (Figure  20e).  The  joint  condition  was  as¬ 
sumed  to  have  a  100  percent  efficiency  of  shear  transfer  but  a  zero  per¬ 
cent  efficiency  for  the  moment  transfer.  The  purpose  of  the  computa¬ 
tions  C  and  D  was  to  determine  which  axle  load  produces  higher  stresses 
in  the  pavement.  The  computed  stresses  and  deflections  are  tabulated 
in  Table  lU  for  the  five  cases.  Different  stress  components  and  deflec¬ 
tions  are  discussed  separately.  The  letters  a,  b,  c,  d,  and  e  in  the 
table  refer  to  the  loading  positions  and  the  pavement  condition  shown 
in  Figure  20. 


a. 


The  stresses  perpendicular  to  the  joint,  .  Table 

lU  shows  that  the  maximum  o  occurs  when  the  trans¬ 
verse  joint  does  not  exist  (column  a).  When  a  trans¬ 
verse  joint  is  present  and  when  the  twin-tandem  load  is 
placed  next  to  the  joint  (column  b),  the  stresses  o 
are  reduced  drastically  (about  50  percent),  and  the  maxi¬ 
mum  stress  occurs  along  the  edge  of  the  rear  axle  load 
away  from  the  joint.  The  maximum  o  under  the  front 
axle  load  (386.6  psi)  is  only  55  percent  of  that  under 
the  rear  axle  load  (709  psi).  Column  c  shows  the  dis¬ 


tribution  of  o  when  the 
b  ft  away  from  ?he  joint. 


twin-tandem  load  is  placed 
Under  node  h9,  the  maximum, 


0^  occurs,  having  a  magnitude  of  1307. ^  psi,  which  is 
very  close  to  the  maximum  stress  (1359.2  psi  in  column  a) 
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in  which  the  transverse  joint  is  not  present  and  the 
load  is  a  full  twin-tandem  load.  This  small  discrepancy 
is  believed  to  be  caused  by  the  small  dimensions  of  the 
slab  length.  In  other  words,  if  the  slab  length  were 
increased  and  the  load  were  placed  farther  away  from  the 
joint,  the  maximum  stress  in  column  c  would  be  increased 
and  reach  the  same  value  as  though  the  joint  did  not 
exist  (column  a).  Column  d  shows  the  distribution  of 
o_  under  the  front  axle  load  alone.  It  is  interesting 
to  note  that  the  maximum  o  does  not  occur  under  any 
load,  but  occurs  at  node  *+9»  which  is  5  ft  away  from  the 
joint.  Column  e  shows  the  uistribution  of  o  under 
the  rear  axle  load  alone.  It  is  also  interesting  to  note 
that  the  maximum  a  has  a  magnitude  of  1259-6  psi, 
which  is  nearly  equal  to  the  maximum  a  (1307.**  psi  in 
column  c )  under  the  twin-tandem  wheel  load  when  the  load 
is  1*  ft  away  from  the  joint,  but  is  much  greater  than 
the  maximum  (709  psi  in  column  b)  when  the  twin- 

tandem  load  is  placed  next  to  the  joint.  This  subject 
will  be  further  discussed  later. 

b.  The  stresses  parallel  to  the  joint,  o  .  Table  lU 

. .  .  - —  .  -  —  -tf 

shows  that  when  the  loads  are  placed  along  the  edge  of 
the  pavement ,  the  stresses  a  are  much  smaller  than 
the  stresses  a  in  all  the  ¥ive  cases,  i.e.,  a,  b,  c, 
d,  and  e  in  Tabfe  1**.  Similar  to  the  distribution  of 
a  ,  maximum  a  occurs  when  the  twin-tandem  load  is 
pfaced  **  ft  awa^  from  the  joint  (column  c).  Unlike  the 
stresses  a  ,  however,  the  stresses  o  induced  by  the 
rear  axle  load  (column  e)  are  smaller  t^an  those  induced 
by  the  front  axle  load  (column  d)  and  are  much  smaller 
than  those  induced  by  the  twin-tandem  load  (column  b). 

c.  The  maximum  principal  stress,  o  .  Since  the  stresses 

—  _ _ _ _ _  max 

are  much  greater  than  the  stresses  o  ,  the  magni¬ 
tudes  of  the  computed  maximum  principal  stresses  are  very 
close  to  those  of  o  .  The  discussions  presented  with 
respect  to  o  are  applicable  to  the  maximum  principal 
stresses  o 

max 

d_.  Vertical  deflections.  Table  1*4  shows  that  the  maximum 
deflection  occurs  at  the  edge  of  the  joint  (node  77) 
when  the  twin-tandem  load  is  placed  next  to  the  joint 
(column  b),  although  the  stresses  are  the  smallest  at 
this  loading  position.  The  maximum  deflection  at  the 
joint  is  0.3*+2  in.,  as  compared  with  the  0.212  and  0.17** 
in.  when  there  is  no  transverse  joint  (column  a)  and 
when  the  load  is  placed  away  from  the  joints  (column  c), 
respectively. 

1**1.  Some  results  presented  in  Table  1**  are  plotted  in  Figure  21. 
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Maxinum  principal  stress 


Figure  21.  Distributions  of  stresses  and  deflections 
along  the  pavement's  edge  under  edge  loads  in  a  two- 
slab  pavement  system,  WESLIQID  program 


The  plotted  results  are  limited  to  nodal  points  U9,  56,  63,  TO,  and  77 
(see  Figure  20c)  along  the  pavement  edge  where  maximum  stresses  and 
deflections  occur.  Critical  stresses  occur  in  the  pavement  near  nodal 
point  56  in  loading  positions  a,  c,  and  e.  Maximum  deflection  occurs 
at  the  pavement  corner  (node  77)  in  the  loading  position  b. 

Center  loads 

lU2.  Results  presented  in  Table  1^  are  for  loads  traveling  along 
the  edge  of  the  pavement.  Computations  were  also  made  for  loads  moving 
along  the  center  line  of  the  pavement.  Figure  22  shows  the  finite  ele¬ 
ment  layout  and  the  two  loading  positions  for  the  computation.  Because 
of  the  symmetry  of  the  loading,  only  one  half  of  the  pavement  was  used 
in  computations.  The  stresses  are  computed  at  the  nodes  where  the 
nodal  numbers  are  shown.  The  efficiency  of  the  joint  is  assumed  to  be 
100  percent  shear  transfer  and  zero  percent  moment  transfer  as  used  for 
edge  load  cases.  The  loading  used  in  the  computation  also  consists  of 
a  twin-tandem  truck  load  with  two  concentrated  loads  representing  a 
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c.  NODAL  NUMBERS  AND  FINITE  ELEMENT  LAYOUT 


Figure  22.  Finite  element  layout  showing  loading  positions  and  nodal 
numbers  for  computations  presented  in  Table  15,  a  two-slab  pavement 
system  with  center  loads,  WESLIQID  program 


one-wheel  load.  The  computed  results  are  presented  in  Table  15  in 
columns  a,  b,  and  c.  The  values  in  columns  a  and  b  refer  to  loading 
positions  a  and  b,  respectively,  in  Figure  22.  The  values  in  column  c 
are  similar  to  those  in  column  b  except  that  the  transverse  joint  is 
removed.  The  significance  of  the  values  tabulated  in  columns  d  and  e 
will  be  explained  later. 

1^3.  Table  15  shows  that  when  the  twin-tandem  load  is  traveling 
along  the  center  of  the  pavement,  the  stresses  perpendicular  to  the 

joint  when  the  load  is  U  ft  away  from  the  joint  (701.0  psi  in  column  b) 
are  107  percent  greater  than  that  when  the  load  is  next  to  the  joint 
( 338 . 14  psi  in  column  a).  Similar  to  the  edge  loading  condition  in  Table 
lU,  the  maximum  occurs  at  a  distance  h  ft  away  from  the  joint  (node 

37)  when  the  twin-tandem  load  is  placed  next  to  the  joint  (Figure  22a). 
The  stresses  in  columns  c  and  d  will  be  discussed  later. 

lUU.  The  maximum  stresses  ,  in  the  direction  parallel  to  the 
joint,  are  nearly  the  same  when  the  load  is  next  to  the  joint  (653.2  psi 
in  column  a)  and  when  the  load  is  away  from  the  joint  (692.9  psi  in 
column  b).  Table  15  shows  that  the  maximum  principal  stress,  which  is  a 
combination  of  °x  >  °y  ’  and  s*iear  stress  1  >  is  slightly  greater 
when  the  load  is  away  from  the  joint  (Tl1*^  psi  in  column  b)  than  when 
the  load  is  next  to  the  joint  (65^.0  psi  in  column  a).  In  other  words, 
unlike  the  case  when  the  load  travels  along  the  pavement  edge,  the  mag- 
nitude  of  stresses  (mainly  the  maximum  principal  stress,  °max  )  in  the 

pavement  changes  only  slightly  when  the  load  is  traveling  along  the  cen¬ 
ter  of  the  pavement.  However,  this  same  conclusion  does  not  hold  true 
for  deflections.  Table  15  shows  that  when  the  load  is  placed  next  to 
the  joint,  the  maximum  deflection  (0.116  in.  in  column  a)  is  nearly 
twice  as  much  as  the  deflection  when  the  twin-tandem  load  is  b  ft  away 
from  the  joint  (0.068  in.  in  column  b). 

1U5.  The  results  tabulated  in  columns  a  and  b  of  Table  15  are 
plotted  in  Figure  23  at  the  longitudinal  line  along  the  nodal  points  12, 
17,  22,  27,  32,  37,  b2,  hi,  52,  and  57-  It  is  seen  that  the  maximum 
stresses  are  nearly  the  same  under  the  two  different  loading  positions, 
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NOTE:  Curve  A:  Loads  placed  next  to  the  joint 

Curve  B:  Loads  are  4  ft  away  from  the  joint 
Loading  position  and  nodal  points  are  shown 
in  Figure  22. 

Figure  23.  Distributions  of  stresses  and  deflections  near 
the  center  line  of  the  pavement  under  center  loads  in  a 
two-slab  pavement  system,  WESLIQID  program 

but  the  deflections  are  much  greater  when  the  loads  are  placed  next  to 
the  joint.  As  stated  before,  greater  deflection  along  the  joint  would 
cause  greater  subgrade  stress  and  consequently  more  damage  in  the  pave¬ 
ment  along  the  joint. 

lh6.  The  computed  values  shown  in  column  c  are  similar  to  those 
shown  in  column  b  except  the  transverse  joint  is  removed.  It  is  seen 


that  the  stresses  are  generally  reduced,  but  very  slightly,  indicating 
that  if  the  load  is  far  enough  away  from  the  joint,  the  stresses  near 
the  load  are  the  same  as  if  the  pavement  joint  does  not  exist. 

lUT.  Column  d  in  Table  15  shows  the  computed  stresses  and  deflec¬ 
tions  with  butted  joints  (i.e.,  thickened  edge  at  the  joint).  The  thick¬ 
ened  section  starts  from  the  midpoint  of  the  slab  and  increases  gradu¬ 
ally  to  10  in.  at  the  joint — an  increase  of  25  percent  in  thickness. 

The  twin-tandem  load  is  placed  only  next  to  the  joint.  Comparing  the 
values  in  columns  a  and  c,  it  is  seen  that  the  stresses  are  reduced 
about  20  percent  because  of  the  increased  pavement  thickness,  but  the 
pavement  deflections  are  reduced  insignificantly. 

148.  Column  e  in  Table  15  shows  the  computed  values  when  the 
modulus  of  subgrade  reaction  k  along  the  same  Joint  as  in  a,  b,  and  d 
is  increased  from  100  to  U00  pci.  The  subgrade  soil  involved  has  a 
width  of  18  in.  at  each  side  of  the  joint.  The  twin-tandem  load  is 
placed  next  to  the  joint.  Compared  to  values  in  column  a,  it  is  seen 
that  the  maximum  principal  stresses  have  changed  insignificantly,  but 
the  deflections  are  reduced  drastically.  It  is  also  noted  that  because 
the  strengthening  of  subgrade  soil  is  limited  to  the  nodal  points  along 
the  joint,  the  maximum  deflections  occur  at  nodal  points  36  and  Ul, 
which  are  about  b  ft  away  from  the  joint. 

1U9.  It  should  be  pointed  out  that  when  the  load  is  traveling 
along  the  center  of  the  pavement,  although  the  maximum  principal  stress¬ 
es  a  are  nearly  the  same  as  when  the  load  is  next  to  or  away  from 
max 

the  transverse  joint,  the  stresses  perpendicular  to  the  joint  are 

nearly  twice  as  large  as  when  the  load  is  away  from  the  transverse 
joint  than  those  when  the  load  is  next  to  the  joint. 

150.  The  discussions  presented  concerning  the  computed  results  in 
Tables  lU  and  15  can  be  summarized  in  the  following  paragraphs. 

151.  Th  most  critical  condition  in  a  rigid  pavement  occurs  when 
there  are  no  transverse  joints  or  cracks  in  the  pavement  and  also  when 
the  load  is  moving  along  its  edge.  In  a  jointed  pavement,  the  critical 
stress  occurs  when  the  load  is  halfway  between  the  joints  at  an  edge  and 
the  stress  has  *a  magnitude  close  to  that  of  a  pavement  with  no  transverse 
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joint.  When  the  load  is  placed  next  to  the  joint,  a  good  load  trans¬ 
fer  device  along  the  joint  can  reduce  pavement  stresses,  but  the  maxi¬ 
mum  stress  at  the  joint  with  a  zero  percent  load  transfer  efficiency  is 
nearly  equal  to  the  critical  stress  when  the  load  is  halfway  between 
the  joints  (and  when  the  distance  between  the  joints  is  long  enough). 

In  other  words,  the  pavement  thickness  cannot  be  reduced  because  of  the 
existence  of  efficient  dowel  bars  in  the  transverse  joint.  When  the 
load  is  moving  along  the  center  of  the  pavement,  the  stresses  are  smal¬ 
ler  and  are  nearly  independent  of  whether  the  load  is  at  the  center  or 

next  to  the  joint  of  the  point.  The  maximum  stresses  of  o  ,  o  ,  and 

x  y 

o  are  701.0,  697-6,  and  71^-9  psi  (column  b  of  Table  15),  respec- 
max 

tively,  when  the  load  is  moving  along  the  center  of  the  pavement,  as 
compared  with  the  corresponding  maximum  stresses  of  1271.8,  233.7,  and 
1271.8  psi  (column  c  of  Table  lU)  when  the  load  is  traveling  along  the 
edge  in  a  jointed  pavement.  Maximum  deflections  occur  when  the  load  is 
next  to  the  transverse  joint  for  both  interior  (0.1l6  in.;  column  a  of 
Table  15)  and  edge  (0.3^2  in.;  column  b  of  Table  lU )  loads.  The  corre¬ 
sponding  maximum  deflections  are  merly  0.068  in.  (column  b  of  Table  15) 
and  0.17^  in.  (column  c  of  Table  lH)  when  the  loads  are  U  ft  away  from 
the  joint.  Greater  pavement  deflections  mean  greater  subgrade  stresses 
and  subsequently  more  subgrade-related  damage  near  the  joint.  Based  on 
the  limited  amount  of  analysis,  it  seems  that  to  reduce  pavement  deflec¬ 
tions  along  the  joint,  it  is  more  beneficial  to  stabilize  the  subgrade 
soil  along  the  joint  than  to  use  butted  joints. 

152.  In  the  design  of  a  rigid  pavement  with  joints,  it  can  be  con¬ 
cluded  that  the  critical  stress  should  be  computed  by  placing  the  loads 
at  pavement  edge  halfway  between  the  joints.  (The  stresses  may  be  re¬ 
duced  by  decreasing  the  pavement  length. )  The  presence  of  an  efficient 
transverse  joint  next  to  the  loads  can  reduce  the  stresses  greatly  but 
increase  the  deflections  along  the  joint  drastically  and  consequently 
increase  the  subgrade  stresses  causing  pavement  damage  along  the  joint. 

153.  Computations  were  made  here  assuming  the  joint  has  a  100 
percent  efficiency  of  shear  transfer,  i.e.,  very  strong  dowel  bars  at 
the  joint.  In  reality,  most  pavements  transfer  less  than  100  percent 
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shear  across  the  joint.  Therefore,  deflections  along  the  loaded  slab's 
joint  should  be  greater  than  those  in  Tables  lb  and  15. 

Effect  of  Temperature  and  Gaps  Under  the  Slabs 

15b,  Temperature  differentials  cause  concrete  pavements  to  warp 
upward  (edge  up)  in  the  night  and  downward  (edge  down)  in  the  day.  Be¬ 
cause  of  the  concrete's  weight,  temperature  stresses  are  produced  in  the 
warped  pavement,  such  as  tensile  stresses  at  the  bottom  of  the  slab  in 
a  downward  warped  pavement.  Depending  upon  the  loading  position, 
stresses  in  the  pavement  can  be  magnified  when  temperature  is  considered. 

155.  Gaps  exist  under  pavements  for  a  number  of  reasons,  such  as 
pumping,  plastic  deformation  of  the  subgrade,  and  poor  construction 
practices.  The  presence  of  gaps  under  the  pavement  is  equivalent  to 

the  removal  of  subgrade  support  for  the  pavement  and  consequently  greater 
stresses  in  the  pavement.  Pavement  stresses  can  be  magnified  when  gaps 
are  combined  with  temperature  warping.  The  following  sections  present 
example  problems  computed  using  WESLIQID  to  show  the  effects  of  tempera¬ 
ture  warping  and  gaps  under  the  slab. 

A  single  slab 

156.  Figure  2b  shows  the  computed  stresses  and  deflections 

in  a  single  concrete  slab  subjected  to  temperature  warping.  A  large 
gap  exists  in  the  subgrade  at  the  slab's  center  under  the  applied  load. 
The  gap  has  a  magnitude  of  1  in.  as  shown  in  the  figure.  The  Young's 
modulus  and  Poisson's  ratio  of  the  concrete  are  6,000,000  psi  and  0.2, 
respectively.  The  weight  of  the  concrete  slab  is  considered  in  the 
computation . 

1 57-  Curve  A  shows  the  distribution  of  the  computed  results  sub¬ 
jected  to  the  applied  load  alone  without  consideration  of  temperature. 

The  slab  is  in  full  contact  with  the  subgrade.  Maximum  stress  occurs 
at  the  slab's  center,  with  tension  at  the  bottom.  Away  from  the  load 
near  the  slab  edge,  however,  small  tensile  stresses  occur  at  the  top 
of  the  slab.  Curve  B  shows  similar  results,  except  the  gap  under  the 
slab  is  considered  in  the  computation.  Distributions  of  curves  A  and  B 
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are  similar,  except  the  magnitude  of  stresses  and  deflections  are 
greater  in  curve  B  because  of  the  existence  of  the  gap  at  slab  center. 

158.  Curve  C  shows  the  distribution  of  stresses  and  deflections 
when  the  gap  and  temperature  are  considered.  The  temperature  differen¬ 
tial  is  -U5°F  between  the  top  and  the  bottom  of  the  slab,  i.e.,  3.75°F 
per  inch  of  the  concrete  slab,  and  the  slab  warps  downward.  It  should 
be  noted  that  a  U5°F  temperature  differential  is  an  extreme  case;  the 
use  of  such  a  large  value  is  for  the  purpose  of  illustration.  Figure  2h 
shows  that  in  spite  of  the  applied  load  at  slab  center,  the  center  por¬ 
tion  of  the  slab  curls  up  due  to  the  large  temperature  differential  and 
loses  contact  with  the  supporting  subgrade.  Along  the  slab  edge,  the 
slab  sinks  into  the  ground.  This  is  called  precompression  in  this  re¬ 
port.  In  the  figure  it  is  seen  that  two  thirds  of  the  slab  (at  the 
center)  has  lost  the  subgrade  contact  and  therefore  produces  very  large 
tensile  stress  under  the  load  and  slab  weight. 

15 9<  Curve  D  shows  the  distributions  of  stresses  and  deflections 
for  the  case  of  positive  temperature  differential  (+i*5^F).  Because  of 
the  large  temperature  differential  and  also  the  existence  of  gap  at 
pavement  center,  one  third  of  the  slab  along  the  edge  curls  above  the 
ground;  the  deflection  at  slab  center  is  also  the  largest  in  all  the 
cases.  However,  the  tensile  stresses  at  slab  center  in  this  case  are 
the  smallest,  because  under  a  positive  temperature  differential,  the 
slab  warps  upward  and  because  of  the  slab  weight,  compression  develops 
at  the  bottom  of  the  concrete  slab,  while  under  the  load  the  bottom  of 
the  slab  is  always  under  tension.  The  effect  of  the  applied  load  and 
positive  temperature  cancels  the  stresses  in  the  concrete  slab.  Fig¬ 
ure  2h  also  shows  that  fairly  large  tensile  stress  developed  at  the  por¬ 
tion  of  the  slab  away  from  the  center  load.  Because  of  the  slab  curling 
up  at  the  pavement  edge  and  because  of  the  weight  of  the  slab  pulling 
the  slab  edge  downward,  tensile  stresses  are  thus  developed  at  the  top 
of  the  concrete  slab. 

A  single  slab:  comparison 
with  the  Westergaard  solution 

l60.  A  complete  subgrade  contact  condition  was  assumed  in  the 
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Westergaard  solution  (1925).  The  slab  always  has  a  full  contact  with 
the  subgrade  soil,  and  gaps  are  not  allowed  between  the  slab  and  sub¬ 
grade,  no  matter  how  much  the  slab  has  warped  upward  due  to  temperature 
change  or  the  applied  load.  In  other  words,  the  slab  is  supported  by 
a  group  of  springs,  and  the  springs  are  always  connected  to  the  slab. 

In  reality,  the  pavement  can  lose  subgrade  support  at  some  points  due  to 
temperature  warping,  pumping,  and  plastic  deformation  of  the  subgrade. 
Computations  were  made  to  demonstrate  the  difference  between  the  solu¬ 
tions  obtained  from  the  Westergaard  solution  and  the  WESLIQID  program. 

In  the  latter  case,  the  slab  can  lose  subgrade  support  from  temperature 
warping. 

l6l.  The  same  slab  and  the  same  loading  used  in  the  computation 
shown  in  Figure  2k  were  used  in  this  computation,  except  that  the  gap 
under  the  pavement  was  removed  because  the  Westergaard  solution  is  not 
capable  of  considering  gap  under  the  pavement.  Figure  25  shows  the 
plotted  results  of  the  computations. 


162.  Curve  A  shows  the  deflections  and  stresses  computed  by  the 
Westergaard  solution  whem  temperature  is  not  considered.  Slight  tensile 
stress  (tension  at  the  top  of  the  slabs)  is  computed  near  the  slab 
edges.  Note  in  the  Westergaard  solution  that  the  slab  is  assumed  to  be 
weightless.  Cruves  B  and  C  show,  respectively,  the  plotted  results 
computed  by  the  Westergaard  solution  and  the  WESLIQID  program.  In  both 
cases,  the  positive  temperature  effect  causing  the  slab  to  warp  upward 
was  considered.  Because  a  complete  subgrade  contact  was  assumed  in  the 
Westergaard  solution,  the  negative  edge  deflections  computed  by  the 
Westergaard  solution  (curve  S)  are  much  smaller  than  those  by  the 
WESLIQID  program  (curve  C)  in  which  the  slab  is  free  to  curl  up  due  to 
the  temperature  change.  On  the  other  hand,  because  of  the  restraint 
imposed  near  the  slab  edges,  the  positive  deflections  (downward  direc¬ 
tion)  near  the  slab  center  computed  by  the  Westergaard  solution  are 
also  smaller  than  those  computed  by  the  WESLIQID  program  where  there  is 
no  restraint  in  the  subgrade  support. 

163.  Although  the  deflections  computed  by  the  WESLIQID  program  are 
much  greater  than  those  of  the  Westergaard  solution  because  of  the 
subgrade  restraint  in  the  latter  case,  the  negative  stresses  computed  by 
the  Westergaard  solution  are  greater  than  those  computed  by  the  WESLIQID 
program,  also  because  of  the  restraint  in  the  subgrade  support.  The 
negative  stresses  shown  in  curve  C  are  induced  by  the  weight  of  the  slab 
because  the  slab's  edges  are  curled  up.  When  the  slab  is  pulled  down 
along  the  edges  by  the  subgrade  reaction  forces  as  in  the  case  of  the 
Westergaard  solution,  greater  negative  stresses  are  developed  in  the 
slab,  as  shown  in  curve  E.  At  the  slab's  center,  stresses  computed  by 
the  Westergaard  solution  ar">  smaller  than  those  computed  by  the  WESLIQID 
program.  The  results  presented  in  Figure  25  demonstrate  the  signifi¬ 
cance  of  subgrade  support  in  a  rigid  pavement  and  the  inadequacies  of 
the  Westergaard  solution  in  some  cases.  Similar  computation  for  negative 
temperature  effect  causing  the  slab  to  warp  downward  was  made,  but  the 
difference  between  the  results  computed  by  the  WESLIQID  program  and  the 
Westergaard  solution  was  not  significant. 
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Stresses  and  deflections  due  to  the 
applied  load  alone  and  the  nonlinear 
response  of  the  pavement  to  the  load 

164.  In  Part  II  of  this  report,  it  was  stated  that  when  the  slab 
and  the  subgrade  are  in  partial  contact,  the  principle  of  superposition 
no  longer  applies;  i.e.,  the  computed  stresses  and  deflections  are  not 
linearly  proportioned  to  the  applied  load.  The  following  example  com¬ 
putation  shows  this  nonlinear  effect. 

165.  Figure  26  shows  the  finite  element  grid  pattern  of  a  square 


slab  subjected  to  a  corner  load.  The  temperature  differential  between 
the  top  and  the  bottom  of  the  slab  is  3°F,  which  causes  the  slab  to 
warp  upward.  Two  computations  were  made  separately.  The  first  one  was 
to  compute  the  stresses  and  deflections  due  to  the  temperature  and  the 
slab  weight,  and  the  second  was  to  compute  those  when  the  effect  of  the 
temperature  and  the  load  are  combined.  The  differences  between  the  com¬ 
puted  results  are  the  stresses  and  deflections  induced  by  the  applied 
load  alone.  It  should  be  noted  that  in  the  field  of  pavement  engineer¬ 
ing  research,  it  is  very  difficult  to  measure  thermal  stresses.  Pave¬ 
ment  engineers  tend  to  measure  stresses  due  to  the  load  only.  When  a 
stress  gage  is  embedded  in  a  test  pavement,  the  stress  induced  by  the 
applied  load  is  derived  from  the  difference  in  the  readings  taken  be¬ 
fore  and  after  the  application  of  the  load.  The  reading  taken  before 
the  application  of  the  load  takes  into  account  the  effect  of  temperature, 
slab  weight,  and  many  other  factors.  The  reading  taken  after  the  appli¬ 
cation  of  the  load  has  the  effect  of  load,  slab  weight,  temperature, 
and  other  factors.  The  difference  in  the  two  readings  is  thus  the  ef¬ 
fect  due  to  the  applied  load  alone. 

166.  Table  16  presents  the  distributions  of  stress  and  deflection 
computed  under  the  various  conditions.  Entry  1  shows  the  initial  curl¬ 
ing  due  to  temperature.  The  computation  is  made  under  the  assumption 
that  the  slab  is  weightless.  Entry  2  shows  the  stresses  and  deflections 
induced  by  temperature,  slab  weight,  and  applied  load.  Two  loads  were 
used  separately  in  computations.  One  had  a  unit  pressure  p  =  100  psi 
and  the  other  p  =  1000  psi  .  It  is  seen  that  the  stresses  and  deflec¬ 
tions  under  the  large  load  (p  =  1000  psi)  are  much  greater  but  not  10 
times  greater  than  those  computed  under  the  smaller  load  (p  =  100  psi); 
i.e.,  the  principle  of  superposition  (as  assumed  by  the  Westergaard 
solution)  does  not  hold.  Under  the  smaller  load,  the  nodes  along  the 
slab  edge  are  still  curled  up  due  to  the  temperature  warping,  but  many 
nodes  near  the  load  are  sunk  into  the  subgrade  under  the  load. 

167.  Entry  3  presents  the  computed  stresses  and  deflections  in¬ 
duced  by  temperature  and  the  slab  weight.  The  effect  of  applied  load 
was  excluded.  Because  the  effects  of  temperature  and  slab  weight  are 
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uniform  on  the  slab,  the  deflections  at  the  corner  nodes  9  and  99  are 
supposed  to  be  equal  because  of  symmetry.  The  slight  difference  in 
the  deflections  in  nodes  9  and  99  is  caused  by  the  nonsymmetrical  lay¬ 
out  of  the  finite  element  grid  and  is  considered  insignificant. 

168.  Entry  k  presents  the  stresses  and  deflections  induced  by  the 
applied  load  alone.  This  is  obtained  by  subtracting  the  values  in  entry 
3  from  entry  2.  Similar  to  entry  2,  the  stresses  and  deflections  in¬ 
duced  by  the  p  =  1000  psi  load  are  not  10  times  greater  than  those 
induced  by  the  p  =  100  psi  load. 

169.  It  should  be  reiterated  that  in  the  field  of  rigid  pavements 
the  design  engineers  are  interested  in  the  stresses  induced  by  the  com¬ 
bined  effect  of  the  applied  load,  temperature,  and  slab  weight.  But 
pavement  engineers  are  interested  in  the  stresses  induced  by  the  applied 
load  alone  because  it  is  the  only  stress  that  engineers  can  measure  with 
confidence.  These  stresses  can  be  computed  separately  by  use  of  the  pro¬ 
grams  developed  in  this  report. 

A  doweled  two- slab  pavement 
system  subjected  to  center 
and  corner  loads  at  the  joint 

170.  In  Part  III  of  this  report,  computations  were  made  to  com¬ 
pare  the  theoretical  computations  with  field  measurements  obtained  in  a 
number  of  airfields  in  the  United  States.  The  variables  compared  were 
the  strains  and  the  stress  transfer  across  the  joint.  The  computations 
were  made  assuming  that  (a)  the  pavements  have  an  initial  full  contact 
with  the  subgrade,  and  (b)  there  are  no  temperature  variations.  The 
agreement  between  the  measurements  and  the  computed  values  were  general¬ 
ly  good.  Computations  were  made  later  for  the  test  pavement  at  Lockbourne 
AFB  (Table  5),  assuming  that  temperature  variations  and  gaps  existed 
along  the  joints.  Besides  the  placement  of  loads  next  to  the  joint  at 
the  interior  of  the  concrete  slab  as  was  done  during  the  tests,  loads 
were  also  placed  at  a  corner  of  the  slab.  The  purpose  of  the  computa¬ 
tions  was  to  examine  the  combined  effect  of  temperature,  in  terns  of  day 
and  night  conditions,  and  gaps  under  the  joint  on  the  stress  and  deflec¬ 
tion  distribution  when  the  load  is  placed  at  different  locations  along 
the  joint. 
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171.  The  finite  element  grid  patterns  for  both  center  and  corner 
loads  are  shown  in  Figure  27.  The  grid  patterns  are  different  because 
of  the  difference  in  loading  positions.  The  deflections  computed  for 
seven  different  conditions  at  nodal  points  along  the  transverse  joint 
are  plotted.  Negative  deflection  indicates  upward  movements.  Curve  1 
presents  the  deflections  due  to  the  loads  only.  Maximum  deflections 
occur  under  the  load,  and  the  maximum  deflection  under  the  corner  load 
is  nearly  three  times  as  great  as  that  in  the  center  load  case. 

172.  Curve  2  shows  the  deflection  distribution  due  to  temperature 
warping  alone.  The  positive  temperature  differential  caused  the  slabs 
to  warp  upward,  leaving  the  pavement  along  the  joint  out  of  contact 
with  the  subgrade.  The  temperature  differential  between  the  top  and 
the  bottom  of  the  pavement  is  3.75°F  per  inch  of  the  pavement,  a  very 
extreme  condition.  Nodes  31  and  U5  in  the  plots  for  the  case  of  center 
and  corner  loads,  respectively,  are  located  at  the  same  location  next 
to  the  joint  along  the  center  line  of  the  pavement. 

173.  Curve  3  shows  the  deflection  curves  similar  to  curve  2  ex¬ 
cept  that  the  loads  are  added  on.  In  the  case  of  center  loads,  the 
center  of  the  pavement  along  the  joint  sinks  into  the  subgrade,  while 
the  pavement  edge  is  still  above  the  ground.  Under  the  corner  load, 
the  center  of  the  pavement  is  above  the  ground  because  of  the  positive 
temperature  differential.  However,  the  edge  of  the  pavement  sinks  into 
the  ground  due  to  the  applied  load. 

17^.  Curve  U  shows  the  deflection  curves  similar  to  curve  3  ex¬ 
cept  that  a  1-in.  gap  under  the  pavement  along  the  joint  is  considered. 
The  gap  has  a  width  of  12  in. ;  it  can  be  caused  by  either  pumping  or 
plastic  deformation  due  to  the  load.  It  is  interesting  to  note  that  in 
the  case  of  the  center  load,  the  deflection  curves  3  and  h  are  very 
nearly  the  same,  because  the  pavement  curls  up  along  the  joint;  thus  the 
gap  beneath  the  Joint  has  no  effect  on  the  pavement  deformation  along 
the  joint.  Under  the  corner  load,  the  shapes  of  curves  3  and  U  are  simi¬ 
lar  except  the  corner  deflection  is  greater  due  to  the  existence  of  the 
gap.  Note  that  the  shape  of  curve  k  would  be  different  if  the  magnitude 
of  the  gap  under  the  pavement  were  less  than  the  maximum  deflection  due 
to  the  load. 
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175-  Curve  5  is  similar  to  curve  2  except  that  the  negative  tem¬ 
perature  differential  caused  the  slabs  to  warp  downward,  leaving  the 
pavement  along  the  joint  to  sink  into  the  subgrade.  Note  that  the  pave¬ 
ment  along  the  edge  sinks  into  the  subgrade  more  than  that  at  the  center. 

176.  Curve  6  is  a  deflection  curve  similar  to  curve  5  except  that 
the  load  is  added  on.  Under  the  center  and  corner  loads,  the  pavement 
along  the  joint  sinks  further  into  the  subgrade,  with  maximum  deflection 
at  a  location  under  the  load. 

177-  Curve  7  is  similar  to  curve  6  except  that  a  1-in.  gap  along 
the  joint  is  considered.  Under  both  the  center  and  the  corner  loads, 
the  pavement  along  the  joint  sinks  into  the  subgrade  more  as  compared 
with  those  in  curve  6  because  of  the  existence  of  the  gap  along  the 
joint.  Also,  it  should  be  noted  that  the  shape  of  curve  7  would  be 
different  if  the  magnitude  of  the  gap  were  less  than  the  deflections 
computed. 

178.  The  stress  distributions  computed  for  the  seven  curves  for 
varied  conditions  shown  in  Figure  27  are  presented  in  Tables  17  and  18 
for  center  and  corner  loads,  respectively.  The  tabulated  stresses  are  at 
nodal  points  along  the  joint  and  near  the  loads.  The  locations  of  the 
numbered  nodes  are  given  in  Figure  27.  In  the  center  load  case  (Table 
17)  the  percent  stress  (or  load)  transfer  was  computed  only  at  nodes  -31 
and  36  under  the  load.  Under  the  corner  load,  however,  the  percent 
stress  (or  load)  transfer  is  difficult  to  determine  according  to  the 
adopted  definition  of  percent  stress  transfer,  because  the  stress  at 
most  nodal  points  changes  signs  across  the  joint. 

179*  In  curves  1-7  shown  in  Figure  27,  positive  stress  indicates 
that  the  tensile  stress  is  developed  at  the  bottom  of  the  pavement,  and 
negative  stress  indicates  that  the  tensile  stress  is  developed  at  the 
top  of  the  pavement.  Positive  temperature  differential  causes  the  pave¬ 
ment  to  curl  upward  (curves  2-M ,  and  because  of  the  weight  of  the  pave¬ 
ment,  negative  stress  is  developed  in  the  pavement.  Negative  temperature 
differential  causes  the  pavement  to  curl  downward  (curves  5-7)  and 
develops  positive  stress  in  the  pavement. 

180.  It  is  interesting  to  note  that  when  temperature  and  gaps  are 
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not  considered,  stresses  under  the  center  load  are  generally  greater 
than  those  under  the  corner  load,  while  the  deflections  are  reverse  in 
the  two  loading  cases. 

l8l.  In  the  case  of  center  load  in  Table  17,  the  stresses  due  to 
load  alone  are  positive  (curve  1  with  tension  at  the  bottom  of  the  pave¬ 
ment).  Since  the  thermal  stresses  are  negative  when  the  pavement  is 
warped  upward  (curve  2  compression  at  the  bottom  of  the  pavement)  and 
are  positive  when  the  pavement  is  warped  downward  (curve  5  with  tension 
at  the  bottom  of  the  pavement),  the  stresses  in  curves  3  and  U  in  which 
the  pavement  is  warped  upward  are  less  than  those  in  curve  1  because 
the  stresses  due  to  temperature  and  stresses  due  to  load  have  different 
signs.  However,  in  curves  6  and  7  in  which  the  pavement  is  warped  down¬ 
ward,  the  stresses  are  much  greater  than  those  in  curve  1,*  because  the 
stresses  due  to  temperature  and  stresses  due  to  load  have  same  sign.  It 
can  be  concluded  that  under  the  center  load,  the  positive  temperature 
differential  (nighttime  condition)  causing  the  pavement  to  warp  upward 
can  reduce  the  stresses,  but  the  negative  temperature  differential  (day¬ 
time  condition)  causing  the  pavement  to  warp  downward  increases  the 
stresses.  It  is  also  noted  in  Table  17  that  the  percent  stress  transfer 
in  conditions  (or  curves)  3  and  U  are  less  than  the  25  percent  under  the 
load  shown  in  curve  1,  but  the  percent  stress  transfer  in  curves  6  and  7 
are  more  than  the  25  percent.  This  is  true  because  when  the  pavement  is 
warping  upward,  the  stresses  on  the  unloaded  side  of  the  joint  reduce 
more  than  those  reduced  on  the  loaded  side.  The  reverse  is  true  when 
the  pavement  is  warped  downward. 

182.  Table  17  also  shows  that  stress  changes  very  slightly 

when  the  gap  under  the  pavement  is  present.  This  can  be  readily  ex¬ 
plained  for  the  case  when  the  pavement  is  warped  upward  (curve  U).  The 
major  portion  of  deflection  curve  *+  shown  in  Figure  27  is  still  above 
the  subgrade,  so  the  existence  of  the  gap  under  the  pavement  does  not 
change  the  stress  distribution  significantly.  In  the  case  of  the  pave¬ 
ment  warping  downward  sinking  into  the  subgrade  (curve  7),  however,  the 

*  The  maximum  edge  stress  in  the  pavement  when  the  same  loads  are  placed 
at  the  pavement  edge  halfway  between  the  joints  is  1151.5  psi. 
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lack  of  change  of  stresses  from  curve  6  in  which  the  gap  is  not  present 
is  difficult  to  explain.  When  the  pavement  is  warped  downward  and  the 
applied  load  is  forcing  the  pavement  into  the  subgrade  (curve  6),  the 
full  subgrade  reactive  forces  are  resisting  the  pavement's  deflection. 
When  the  subgrade  reactive  forces  are  removed  as  in  curve  7  in  which  the 
gap  exists  along  the  joint,  the  pavement  is  allowed  to  deflect  more 
freely  and  consequently  pavement  stresses  should  be  increased.  However, 
this  is  apparently  not  the  case  because  the  computed  stresses  listed  for 
curve  7  are  nearly  the  same  as  those  computed  in  curve  6.  The  only 
explanation  in  this  case  is  that  deflection  curves  6  and  7  shown  in 
Figure  27  (under  the  center  loads)  are  nearly  parallel  to  each  other. 
Since  stresses  in  a  pavement  are  induced  by  a  change  in  the  curvature 
of  the  pavement,  the  stresses  thus  remain  the  same.  It  is  also  noted  in 
Figure  27  that  deflection  curves  3  and  are  parallel  to  each  other, 
explaining  why  the  stresses  do  not  change  from  curve  3  to  curve  h . 

183.  It  should  be  emphasized  that  deflection  curves  6  and  7  may 
not  be  parallel  to  each  other  if  either  the  modulus  of  subgrade  reaction 
k  is  increased  or  the  pavement  thickness  is  decreased.  Under  such 
conditions,  the  shape  of  deflection  curve  6  may  change  more  than  that 

of  curve  7- 

184.  Table  l8  shows  the  stress  distribution  under  the  corner 
loads.  The  stresses  are  tabulated  at  nodal  points  around  the  load  and 
are  generally  smaller  than  those  computed  under  the  center  loads  shown 
in  Table  17-  Unlike  the  center  load  case,  the  effect  of  temperature 
warping  on  the  stress  distribution  in  the  corner  load  case  is  not  sig¬ 
nificant.  Under  temperature  warping,  the  stresses  increase  in  some 
nodal  points  but  decrease  in  some  others.  The  changes  were,  however, 
not  significant . 

185.  It  should  be  pointed  out  that  the  presence  of  the  load 
transfer  dowel  bars  can  reduce  the  stresses  in  the  pavement  near  the 
joint,  and  the  reduction  is  generally  more  in  the  corner  load  cases 
than  in  the  center  load.  To  clarify  this  point,  computations  were  made 
for  the  test  pavements  in  Lockbourne  AFB,  under  the  assumption  that 
there  was  no  load  transfer  across  the  joint,  and  the  results  are 


100 


presented  in  Table  19.  For  simplicity,  the  comparison  was  made  only 
for  the  case  when  the  temperature  effects  and  the  gaps  under  the  pave¬ 
ment  are  not  considered. 

186.  Table  19  shows  that  pavement  stresses  (o  )  are  drastically 

max 

increased  when  the  dowel  bars  at  the  joint  are  removed,  and  the  effect 
is  more  pronounced  in  the  corner  load  case.  Shear  stresses  t  play  a 
more  important  role  when  the  loads  are  applied  at  the  corner  of  the 
pavement  than  when  the  loads  are  applied  at  the  pavement  center  where 
the  shear  stress  is  zero  because  of  symmetry. 

187.  Results  presented  in  Figure  27  and  Tables  17,  18,  and  19 
for  a  two-slab  pavement  system  with  dowel  bars  subjected  to  corner  and 
center  loads  at  the  joint  are  summarized,  and  conclusions  are  made  in 
the  following  paragraphs. 

188.  Whether  temperature  and  gaps  under  the  joint  are  considered 
or  not,  aircraft  load  placed  at  the  corner  produces  much  larger  deflec¬ 
tions  than  when  the  load  is  placed  at  pavement  center  (curve  1  of 
Figure  27).  When  the  pavement  is  warped  downward  (daytime  condition), 
the  differences  in  deflections  between  the  corner  and  center  loads  be¬ 
come  even  greater. 

189.  When  temperature  and  gaps  are  not  considered,  stresses  under 
the  center  load  (curve  1  of  Table  17)  are  generally  greater  than  those 
under  the  corner  load  (curve  1  of  Table  18).  Under  the  center  load, 
when  the  pavement  is  warped  upward  (nighttime  conditions),  stresses  in 
the  concrete  pavement  can  be  reduced;  the  opposite  is  true  when  the 
pavement  is  warped  downward  (daytime  condition).  The  effect  of  tem¬ 
perature  is  not  so  significant  and  clear  when  the  load  is  placed  at  the 
corner. 

190.  The  presence  of  dowel  bars  in  the  joint  can  reduce  the 
stresses  in  the  pavement  near  the  joint.  The  reduction  is  more  signi¬ 
ficant  in  the  corner  load  case. 

Continuously  Reinforced  Concrete  Pavement  (CROP) 

191.  Results  presented  in  Table  lU  indicate  that  in  a  jointed 
pavement,  maximum  stresses  occur  when  the  edge  load  is  placed  halfway 
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between  the  joints,  and  the  stresses  can  reach  the  maximum  stresses 
that  occur  in  a  pavement  without  transverse  joints  if  the  pavement  is 
long  enough  in  the  longitudinal  direction.  The  question  arises  as  to 
whether  the  pavement  stress  would  be  reduced  with  a  reduction  in  pave¬ 
ment  length.  It  is  known  that  in  CRCP  closely  spaced  transverse  cracks 
developed  shortly  after  the  construction,  but  the  cracks  are  held  very 
tightly  by  the  steel  reinforcements.  It  is  reasonable  to  assume  that 
the  cracks  can  transfer  100  percent  shear  force  and  a  small  amount  of 
moment  when  the  pavement  is  in  sound  structural  condition.  Computations 
were  thus  made  in  a  CRCP  to  compare  the  computed  values  with  values  for 
a  regularly  jointed  pavement.  The  CRCP  used  in  the  computation  is  a 
five-slab  pavement  divided  by  four  cracks.  For  the  WESLIQID  program, 
the  cracks  were  treated  the  same  as  the  joints.  The  loading  condition 
and  pavement  dimensions  for  the  CRCP  are  the  same  as  those  used  in  the 
previous  examples  of  the  regular  jointed  pavement,  except  the  pavement 
length  is  reduced  to  b  ft.  The  joint  (or  crack)  condition  is  also 
identical,  i.e.,  100  percent  shear  transfer  and  zero  percent  moment 
transfer.  The  finite  element  layout  is  shown  in  Figure  28,  and  the  com¬ 
puted  results  are  presented  in  Table  20.  The  discussions  on  the  com¬ 
parisons  between  Table  20  (CRCP  with  closely  spaced  cracks)  and  Table 
lU  (regularly  jointed  pavement  with  edge  loads)  and  Table  15  (regularly 
jointed  pavement  with  center  loads)  are  presented  as  follows.  It  should 
be  pointed  out  that  the  maximum  stresses  may  change  signs  when  condi¬ 
tions  change;  the  comparisons  should  be  based  on  the  absolute  values. 


0^  ,  stresses  per- 
pendicular  to  the  joint 

192 .  Edge  loads  placed  next  to  the  joint  (loading  position  c). 

When  the  crack  spacing  is  reduced  (or  the  number  of  cracks  is  increased) 
as  in  a  CRCP,  the  maximum  stress  is  reduced  from  709. 0  psi  (column  b  of 
Table  1*0  to  1+1+5. 1+  psi  (column  c  of  Table  20).  It  is  noted  that  in 
column  d  of  Table  lU  when  the  front-axle  loads  are  placed  next  to  the 
joint  at  the  pavement's  edge,  the  maximum  occurs  at  a  point  far 

away  from  the  joint.  It  is  thus  reasonable  to  find  that  when  the  pave¬ 
ment  length  is  reduced,  such  as  in  a  CRCP,  the  stress  is  reduced. 
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When  the  crack  spacing  is  reduced,  the  maximum  stress  is  reduced  from 
1307.^  psi  (column  c  of  Table  lU )  to  735-6  psi  (column  a  of  Table  20). 
This  is  reasonable  because  the  results  presented  in  columns  b  and  c  of 
Table  2^  show  that  the  stresses  are  greatly  decreased  when  the 

edge  loads  are  moved  toward  the  joint. 

19U .  Center  loads  placed  next  to  the  joint  (loading  position  d). 
The  maximum  stress  is  reduced  from  338.^  psi  (column  a  of  Table  15)  to 
168.7  psi  (column  d  of  Table  20). 

195 -  Center  loads  placed  between  the  joints  (loading  position  b). 
The  maximum  stress  is  reduced  from  701  psi  (column  b  of  Table  15 )  to 
U2U.1+  psi  (column  b  of  Table  20). 

oy  ,  stresses  parallel  to  the  joint 

196.  Edge  loads  placed  next  to  the  joint  (loading  position  c). 
When  the  crack  spacing  is  reduced,  the  maximum  stress  has  changed  very 
little — from  -216.7  psi  (column  b  of  Table  1*+)  to  -218.9  psi  (column  c 
of  Table  20 ) . 

197.  Edge  loads;  placed  between  the  joints  (loading  position  a). 
The  maximum  stress  also  has  not  changed  much  from  -233.7  psi  (column  c 
of  Table  lU)  to  -22U.7  psi  (column  a  of  Table  20 ). 

198.  Center  loads  placed  next  to  the  joint  (loading  position  d). 
The  maximum  stress  has  changed  very  little — from  653. U  psi  (column  a  of 
Table  15)  to  6^6.5  psi  (column  d  of  Table  20). 

199 .  Center  loads  placed  between  the  joints  (loading  position  b). 
The  maximum  stress  is  also  slightly  reduced  from  697.6  psi  (column  b  of 
Table  lh)  to  6l8.3  psi  (column  b  of  Table  20). 


o  ,  maximum  principal  stress 
max  _ _ _ _ _ 

200.  Edge  loads  placed  next  to  the  joint  (loading  position  c). 

When  the  pavement  length  is  reduced,  the  maximum  principal  stress  is 
reduced  from  709.0  psi  (column  b  of  Table  lk)  to  -536.0  psi  (column  c  of 
Table  20). 


The  maximum  principal  stress  is  drastically  reduced  from  1307 . t  psi 
(column  c  of  Table  lU)  to  735*6  psi  (column  a  of  Table  20). 

202.  Center  loads  placed  next  to  the  .joint  (loading  position  d). 
The  maximum  principal  stress  has  changed  very  little — from  65^.2  psi 
(column  a  of  Table  15)  to  651. 1  psi  (column  d  of  Table  20). 

203.  Center  loads  placed  between  the  .joints  (loading  position  b). 
The  maximum  stress  is  reduced  from  71^.9  psi  (column  b  of  Table  15)  to 
619.^  psi  (column  b  of  Table  20). 

w  ,  vertical  deflections 

20U .  Edge  loads  placed  next  to  the  .joint  (loading  position  c). 

When  the  pavement  length  is  reduced,  the  maximum  deflection  is  reduced 
from  0.3^2  in.  (column  b  of  Table  1^)  to  0.185  in.  (column  c  of  Table 
20).  The  maximum  deflections  in  both  cases  occur  at  the  corner  node  of 
the  joint. 

205.  Edge  loads  placed  between  the  joints  (loading  position  a). 

The  maximum  deflection  is  increased  from  0.17^  in.  (column  c  of  Table 
1*0  to  0.183  in.  (column  a  of  Table  20 ).  It  is  to  be  noted  that  while 
the  deflection  is  increased,  the  0.l83-in.  maximum  deflection  is  much 
less  than  the  0.3^2  in.  (column  b  of  Table  1*0  produced  by  the  loads 
placed  ne^t  to  the  joint. 

206.  Center  load  placed  next  to  the  joint  (loading  position  d). 

The  maximum  deflection  is  reduced  from  0.116  in.  (column  a  of  Table  15) 
to  0.087  in.  (column  d  of  Table  20).  The  maximum  deflections  in  both 
cases  occur  at  nodes  along  the  joint. 

207.  Center  load  placed  between  the  joints  (loading  position  b). 
The  maximum  deflection  is  increased  from  0.069  in.  (column  b  of  Table 
15)  to  O.O85  in.  (column  b  of  Table  20).  Similar  to  the  maximum  deflec¬ 
tion  produced  by  the  edge  loads  placed  between  the  joints,  while  the 
deflection  is  increased  when  the  pavement  length  is  reduced,  the  0.085- 
in.  maximum  deflection  is  much  less  than  the  0.3^2  in.  (column  b  of 
Table  1*0  produced  by  the  edge  loads  placed  next  to  the  joint. 

208.  The  values  of  the  maximum  stress  and  deflections  computed  for 
different  loading  positions  for  both  the  regularly  jointed  pavement  and 
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the  CRCP  are  tabulated  in  Table  21.  It  is  seen  that  the  stress  o 

x 

(in  the  direction  perpendicular  to  the  transverse  joint)  is  drastically 
reduced  in  all  the  four  loading  positions,  but  the  stress  a  (in  the 

y 

direction  parallel  to  the  joint)  changes  very  little.  The  maximum 

principal  stress  o  is  reduced  greatly  in  the  edge  load  condition 

max 

but  is  changed  insignificantly  in  the  center  load  conditions.  The 
vertical  deflection  w  is  reduced  nearly  50  percent  when  the  edge  loads 
are  placed  next  to  the  joint  but  is  changed  insignificantly,  as  far  as 
the  pavement  design  is  concerned,  in  the  other  conditions. 

209-  Since  in  practice  the  thickness  of  the  CRCP  is  less  than 
that  of  regular  rigid  pavement,  computations  were  also  made  for  the 
same  CRCP,  but  the  thickness  of  the  rigid  pavement  was  reduced  from  8  to 
6  in.,  and  the  maximum  values  computed  are  presented  in  Table  21.  It  is 
seen  that  when  the  thickness  of  the  CRCP  is  reduced  by  2  in.,  while  the 
stresses  and  deflections  are  increased  greatly,  the  values  do  not  exceed 
those  in  an  8-in.  regularly  jointed  pavement. 

210.  It  should  be  pointed  out  that  the  computations  presented  in 
Tables  20  and  21  for  the  CRCP  are  based  on  the  condition  that  the  cracks 
(or  joints)  have  zero  efficiency  of  moment  transfer.  In  fact,  a  certain 
amount  of  moment  should  be  transferred  across  the  cracks  when  the  pave¬ 
ment  is  in  good  condition.  Therefore,  the  stresses  and  deflections  can 
actually  be  smaller  than  those  computed  values  shown  in  Tables  20  and  21. 

211.  From  the  discussions  presented  above,  it  can  be  concluded 
that  as  far  as  the  stresses  and  deflections  are  concerned,  it  is  advan¬ 
tageous  to  reduce  the  joint  spacing,  such  as  with  a  CRCP,  provided  that 
the  joints  (or  cracks)  are  kept  in  good  condition. 

212.  The  computations  for  the  CRCP  are  limited  to  l*-ft  crack 
spacing,  and  the  cracks  are  assumed  to  be  100  percent  efficient  in  the 
transfer  of  shear  forces.  In  reality,  when  the  traffic  repetition  in 
a  CRCP  is  increased,  the  number  of  cracks  would  be  increased,  and  the 
efficiency  of  load  transfer  across  the  cracks  would  also  be  reduced. 

When  punchout  distresses  become  imminent,  the  opening  of  the  cracks 
become  so  large  that  the  transfer  of  load  across  the  crack  becomes 
negligible.  To  investigate  this  case,  computations  were  made  for  the 
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same  CRCP  under  the  same  twin-tandem  truck  load  for  the  following  two 
cases:  (a)  the  crack  opening  is  reduced  from  4  to  2  ft,  but  the  effi¬ 

ciencies  of  shear  and  moment  transfer  are  still  assumed  to  be  100  and 
zero  percent,  respectively,  and  (b)  the  crack  spacing  is  2  ft,  but  the 
efficiencies  of  shear  and  moment  transfer  are  both  assumed  to  be  zero 
percent,  i.e.,  the  cracks  become  so  large  that  there  is  no  load  transfer 
across  the  cracks.  The  CRCP  used  in  the  computation  is  a  seven-slab 
pavement  divided  by  six  cracks.  The  computations  were  limited  to  two 
loading  positions:  (a)  the  edge  load  placed  next  to  the  joint  and 
(b)  the  center  load  placed  next  to  the  joint.  Computations  for  loads 
placed  between  the  joints  were  not  done  because  the  crack  spacing  is 
too  small  for  a  twin-tandem  load.  The  computed  maximum  stresses  and 
deflections  for  the  two  loading  positions  are  presented  in  Table  22. 

For  comparison,  the  computed  maximum  stresses  and  deflections  for  the 
other  two  pavements,  i.e.,  jointed  concrete  pavement  with  15-ft  slab 
lengths  and  the  CRCP  with  4-ft  crack  spacing,  were  included. 

21'j.  It  should  be  pointed  out  that  when  the  crack  spacing  is  re¬ 
duced  from  4  to  2  ft,  the  two  sets  of  dual -wheel  loads  of  the  twin- 
tandem  truck  load  are  placed  in  two  different  slabs.  When  the  crack 
opening  becomes  large  enough  that  there  is  no  load  transfer  across  the 
crack,  this  2-ft-long  slab  is  subjected  to  one  set  of  dual -wheel  load 
and  receives  no  effect  from  the  other  set  of  dual  wheels.  Mechanically, 
the  2-ft-long  slab  has  lost  the  advantage  of  the  slab  action  and  reacts 
to  the  load  as  a  2-ft-wide  beam.  Consequently,  the  stresses  in  the 
pavement  would  increase  drastically. 

214.  Table  22  shows  that  under  both  the  edge  and  center  loads, 
when  the  crack  spacing  is  reduced  from  4  to  2  ft  (cases  b  and  c),  the 
maximum  stress  in  the  direction  perpendicular  to  the  transverse 

joint  is  reduced,  but  the  maximum  stress  o  is  increased.  The  maxi- 

y 

mum  principal  stress  in  this  case  is  increased  partly  due  to  the  in¬ 
crease  of  shear  stress  (as  shown  in  the  computer  output).  As  expected, 
the  maximum  deflection  increases  drastically  as  the  crack  spacing  is 
reduced.  It  should  also  be  noted  that  in  case  a  the  maximum  deflection 
in  a  regular  jointed  pavement  is  0.342  in.,  which  is  greater  than  those 
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in  the  CRCP.  The  reason  for  this  is  that  in  a  regular  jointed  pavement, 
the  twin-tandem  load  is  placed  in  one  slab,  while  in  the  CRCP,  the  two 
sets  of  the  dual  wheels  are  placed  on  separate  slabs. 

215.  Computed  values  shown  in  columns  c  and  d  of  Table  22  show 

that  in  both  the  center  and  edge  load  cases,  the  maximum  principal  stress 

0  (and  0  )  and  deflections  are  more  than  nearly  double  when  the 
max  y 

efficiency  of  shear  transfer  across  the  crack  is  reduced  from  100 
(column  c)  to  zero  percent  (column  d).  The  significance  of  the  increase 
in  stress  is  further  discussed  in  the  following  paragraph. 

216.  It  is  a  well-known  fact  that  maximum  stresses  occur  at  the 
edge  of  a  rigid  pavement  at  a  direction  parallel  to  the  edge  of  the 
slab  (o  in  Figure  20a)  under  an  edge  load.  When  the  slab  length  is 
reduced,  i.e.,  more  cracks  have  developed  and  lost  the  load  transfer 
capability,  the  computed  results  indicate  that  the  location  of  the  maxi¬ 
mum  stress  is  no  longer  at  the  pavement  edge  but  has  moved  toward  the 
pavement  interior,  and  in  particular  the  magnitude  of  the  stress  has  in¬ 
creased  drastically.  In  this  case,  the  critical  stress  is  no  linger 

0  but  0  (or  a  ).  The  increase  in  the  number  of  cracks  also 
x  max  y 

has  a  profound  effect  on  pavement  stresses  in  the  center  load  case. 

217.  Figure  29  shows  the  distributions  of  stress  0^  across  the 
12- ft  pavement  cross  section  under  the  center  and  edge  loadings.  While 
Figure  29  shows  that  the  stress  0^  is  smallest  in  the  regular  jointed 
pavement  (case  a) ,  it  should  be  noted  that  the  most  critical  stress  in 
a  concrete  pavement  occurs  at  the  pavement's  edge  in  the  direction 
perpendicular  to  the  transverse  joint  (0  )  when  the  edge  load  is  placed 
between  the  joints  (1307.1*  psi  shown  in  Table  22  and  column  c  of  Table 
lU ) .  Also,  for  a  CRCP  with  U-ft  crack  spacing,  maximum  stress  occurs 
also  at  the  pavement's  edge  in  the  direction  perpendicular  to  the  joint 
(a  )  when  the  edge  load  is  placed  between  the  cracks  (735.6  psi  in 
Table  22  and  column  a  of  Table  20).  The  primary  purpose  of  presenting 
the  results  in  Figure  29  is  to  illustrate  the  important  fact  that  in  a 
CRCP  layer,  stresses  have  shifted  from  the  pavement's  edge  toward  the 
pavement's  center  (at  a  location  about  5  ft  away  from  the  edge  in  this 


218.  Figure  29  shows  that  for  a  2-ft  crack  spacing  in  the  pave¬ 
ment  when  the  efficiency  of  shear  transfer  across  the  cracks  is  reduced 
to  zero,  the  maximum  stress  increases  from  88!t.8  psi  (curve  C)  to  1878.7 
psi  (curve  D)  under  the  center  load  and  from  -537-9  psi  (curve  C)  to 
-1702.6  psi  (curve  D)  under  the  edge  load.  It  is  to  be  noted  that  in 
curve  D,  the  maximum  a  under  the  center  load  is  nearl .  the  same  as 

_ jl _ : _ 

that  under  the  edge  load,  which  is  contrary  to  the  well-established  fact 
that  the  edge  load  is  more  critical  than  the  center  load  in  a  regular 
jointed  concrete  pavement  (see  Tables  lit  and  15). 

219*  Computations  were  also  made  for  pavement  condition  D  of 
Table  22  and  it  was  assumed  that  the  pavement  had  lost  subgrade  contact 
due  to  pumping.  In  the  edge  load  case,  the  nodes  along  the  pavement 
edge  were  assumed  to  be  out  of  contact.  The  maximum  stress  a  was 

y 

computed  to  be  -3^98-7  psi,  which  is  twice  as  much  as  that  if  the  pave¬ 
ment  had  full  subgrade  contact.  The  situation  can  be  worse  if  the  tem¬ 
perature  differential  is  considered.  In  the  center  load  case,  the  nodes 
along  the  cracks  were  assumed  to  be  out  of  contact.  The  stresses  com¬ 
puted  were  not  significantly  different  from  those  computed  when  full 
contact  was  assumed. 

220.  It  is  interesting  to  note  in  Figure  29  that  under  the  edge 
load,  negative  stresses  are  developed  near  the  center  of  the  12-ft-wide 
concrete  pavement  and  that  under  the  center  load,  positive  stresses  are 
developed  directly  under  the  load.  The  question  arises  as  to  whether 
the  actual  stresses  in  a  concrete  pavement  are  smaller  than  those  indi¬ 
cated  in  Figure  29  because  the  stresses  induced  by  the  left  and  right 
truck  loads  may  cancel  each  other.  Computations  were  made  to  verify 
this  concept,  and  the  results  are  presented  in  Figure  30.  Two  series 
of  computations  were  made,  one  where  the  truck  travels  at  the  pavement 
center  and  the  other  where  the  truck  travels  at  the  pavement  edge.  The 
total  width  of  the  truck  from  the  edge  of  one  tire  to  the  edge  of  the 
other  is  assumed  to  be  8  ft  and  the  width  of  each  twin-tandem  load  is 
2  ft,  as  used  in  previous  computations.  The  CROP  is  assumed  to  have 
2-ft  crack  spacing. 


,  PARALLEL  TO  THE  TRANSVERSE  CRACK,  psi 


Center  load 


221.  Figure  30a  shows  the  distribution  of  stress  a  along  the 
transverse  crack  when  the  center  of  the  truck  is  in  line  with  the  center 
of  the  pavement.  Curve  A-A  shows  the  distribution  under  the  left  twin- 
tandem  load  alone.  A  maximum  positive  stress  of  1027*5  psi  was  computed 
under  the  inner  wheel  and  some  negative  stresses  computed  at  the  other 
half  of  the  12-ft  lane.  Curve  B-B  shows  the  distribution  under  the  full 
truck  load,  i.e.  ,  both  left  and  right  twin-tandem  loads.  The  distribu¬ 
tion  is  symmetrical  with  respect  to  the  pavement's  center,  and  the  max¬ 
imum  stress  is  reduced  from  the  1027-5  psi  (curve  A-A)  to  895*5  psi. 

The  reduction  of  stresses  is,  of  course,  caused  by  the  cancellation  of 
negative  stresses  shown  in  curve  A-A. 

Edge  load 


222.  Figure  30b  shows  the  distributions  when  the  edge  of  the  left 
twin-tandem  load  is  at  the  edge  of  the  pavement.  Curve  A-A  shows  the 
distribution  under  the  left  twin-tandem  load  placed  at  the  pavement's 
edge  that  causes  the  entire  pavement  along  the  crack  to  have  negative 
stresses  with  a  maximum  stress  of  1702.6  psi.  Curve  B-B  shows  the  dis¬ 
tribution  under  the  full  truck  load.  Because  the  right  side  twin-tandem 
load  produces  positive  stresses  in  most  parts  of  the  pavement  along  the 
crack,  the  combined  effect  due  to  both  left  and  right  loads  reduces  the 
pavement  stresses  as  shown  in  curve  B-B.  It  is  interesting  to  note  that 
in  curve  B-B  under  either  the  center  or  the  edge  load,  the  maximum 
stress  occurs  at  a  section  from  3  to  5  ft  away  from  the  pavement's 
edge  and  that  the  center  load  causes  the  pavement  to  have  positive 
stress  and  the  edge  load  causes  the  pavement  to  have  negative  stress. 
This  is  believed  to  be  the  primary  cause  of  the  punchout  failure 
frequently  observed  in  highway  CRCP. 

223*  This  study  found  that  the  maximum  stress  in  a  highway  CRCP 
does  not  occur  at  the  pavement’s  edge,  but  occurs  at  the  interior  part 
of  the  slab.  This  was  also  concluded  by  Darter,  LaCoursiere,  and 
Smiley  (1979)  in  their  computation  of  a  12-ft-wide  highway  CRCP  sub¬ 
jected  to  an  l8-kip  axle  load.  They  found  that  when  the  crack  spacing 
has  reduced  to  1  or  2  ft  and  the  crack  opening  has  become  so  large  that 
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no  load  is  transferred  across  the  cracks,  the  most  critical  tensile 
stress  developed  in  the  pavement  at  a  location  about  U— 6  ft  away  from 
the  pavement's  edge.  They  concluded  that  this  critical  stress  can  be 
the  main  cause  for  the  punchout  failures  frequently  observed  in  CRCP's 
in  the  State  of  Indiana. 

22b.  Results  presented  in  Figures  29  and  30  and  in  Tables  13, 

15,  20,  21,  and  22  are  summarized  and  conclusions  are  made  in  the  fol¬ 
lowing  paragraphs. 

225.  In  a  regular  concrete  pavement  with  joints  spaced  at  15  ft 
or  more,  the  most  critical  condition  occurs  when  the  single  twin-tandem 
highway  track  load  is  placed  at  the  pavement's  edge  away  from  the  joints 
(l307.**-psi  edge  stress  shown  in  column  c  of  Table  lU).  Compared  with 
a  regular  jointed  concrete  pavement,  the  stresses  and  deflection  in  a 
CRCP  are  reduced  if  the  pavement  is  in  good  condition.  However,  when 
the  number  of  cracks  is  increased  and  consequently  the  crack  spacing 
is  reduced,  such  as  the  2-ft  case  presented  in  the  computations  of  this 
report,  not  only  the  stresses  start  to  increase,  but  also  the  location 
of  the  maximum  stress  (under  the  edge  load)  shifts  from  the  pavement's 
edge  to  the  interior  part  of  the  pavement.  As  the  load  transfer  capa¬ 
bility  across  the  crack  is  reduced  to  zero  due  to  the  increased  load 
application,  the  stresses  in  the  pavement  increase  drastically.  It  is 
important  to  realize  that  under  a  single  twin-tandem  load,  the  maximum 
stress  is  no  longer  at  the  pavement's  edge  due  to  the  edge  load  as 
generally  known,  but  at  locations  directly  under  the  center  loads. 

This  drastic  increase  in  stresses  at  locations  away  from  the  pavement’s 
edge  induced  by  either  the  edge  load  or  the  center  load  may  be  the 
effect  responsible  for  the  punchout  failures  frequently  observed  in  a 


highway  CRCP  when  the  number  of  cracks  becomes  larger  and  when  the 


twin-tandem  load  placed  at  either  the  pavement's  edge  or  the  pave¬ 
ment's  center  due  to  the  cancelling  effect,  since  the  edge  load  causes 
the  pavement  to  have  negative  stresses  and  the  center  load  to  have 
positive  stresses.  The  maximum  stress  occurs  at  a  section  about  3  to 
5  ft  away  from  the  pavements *s  edge. 


llU 


PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


227.  The  computer  programs  WESLIQID  and  WESLAYER  have  the  capa¬ 
bilities  of  calculating  the  stresses  and  deflections  in  a  rigid  pavement 
with  cracks  and  joints.  WESLIQID  is  designed  for  pavements  on  a  liquid 
foundation,  and  WESLAYER  is  for  pavements  on  layered  linear  elastic 
solids. 


228.  Comparisons  were  good  between  the  computed  results  of 
WESLIQID  and  those  of  available  solutions,  such  as  the  Westergaard 
solution,  Pickett  and  Ray's  influence  charts,  and  the  discrete  element 
method.  The  comparisons  between  the  percent  stress  (or  load)  transfer 
across  the  joint  computed  by  the  WESLIQID  and  WESLAYER  programs  and 
those  measured  in  a  series  of  full-scale  test  sites  were  also  good. 

229.  The  following  conclusions  were  derived  based  on  the  results 


of  the  analysis  of  many  pavements  with  the  WESLIQID  program: 


a. .  The  dowel  is  not  a  very  efficient  load  transfer  device 

across  the  joint  in  a  rigid  pavement.  The  current  Corps 
of  Engineers  practice  of  assuming  only  a  25-percent 
stress  (or  load)  transfer  in  the  design  and  evaluation 
of  rigid  pavements  is  warranted. 

b.  The  stress  transfer  across  a  joint,  as  modeled  in  this 
code,  is  influenced  by  many  factors,  most  drastically  by 
the  modulus  of  dowel  support  (or  subgrade  elastic  modu¬ 
lus)  and  the  number  and  size  of  steel  bars  in  the  joint 
and  to  a  lesser  degree  by  the  modulus  of  concrete,  joint 
spacing,  pavement  thickness,  and  modulus  of  subgrade 
reaction. 


While  vertical  deflections  are  easy  to  measure  in  a  rigid 
pavement,  one  should  be  cautious  in  the  use  of  deflec¬ 
tions  across  a  joint  to  estimate  the  efficiency  of  stress 
transfer  because  unless  the  pavement  is  failed,  the  de¬ 
flections  across  a  joint  do  not  differ  as  much  as  the 
stress  Oy  along  the  joint.  Stress  0  is  in  the  di¬ 
rection  parallel  to  the  joint.  ^ 


d.  The  efficiency  of  stress  transfer  across  a  joint  has  an 
insignificant  effect  on  the  stresses  and  deflections  in 
the  slab  when  the  load  is  placed  at  the  center  of  the 
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slab,  but  has  a  significant  effect  when  the  load  is 
placed  next  to  the  joint. 

e_.  When  the  slab  is  in  partial  contact  with  the  subgrade 
due  to  temperature  warping  or  other  causes,  the  princi¬ 
ple  of  superposition,  as  assumed  by  the  Westergaard  solu¬ 
tion,  is  no  longer  valid;  i.e.,  the  pavement  stresses  and 
displacements  are  not  linearly  proportional  to  the  load, 
even  though  the  slab  is  still  within  the  elastic  range. 

_f.  The  most  critical  condition  in  a  rigid  pavement  occurs 
when  there  is  no  transverse  joint  or  crack  in  the  pave¬ 
ment  and  also  when  the  load  is  moving  along  its  edge. 
Under  the  edge  load,  the  presence  of  joints  and  cracks 
can  reduce  pavement  stresses  near  the  joint  and  the 
cracks  but  increase  the  deflections  in  the  same  area. 

In  a  jointed  pavement,  the  critical  stress  occurs  when 
the  load  is  halfway  between  the  joints  and  the  stress 
can  have  a  magnitude  close  to  that  of  a  pavement  with  no 
transverse  joint.  Therefore,  the  presence  of  a  joint 
does  not  reduce  the  maximum  stress  in  a  pavement.  When 
the  load  is  moving  along  the  center  of  the  pavement,  the 
stresses  are  smaller  and  are  nearly  independent  of 
whether  the  load  is  at  the  center  or  next  to  the  pavement 
joint.  Under  both  center  and  edge  loads,  maximum  deflec¬ 
tion  occurs  when  the  load  is  next  to  the  transverse 
joint.  Greater  pavement  deflections  induce  greater  sub¬ 
grade  stresses  and  consequently  more  severe  plastic  de¬ 
formation  in  the  subgrade  soil,  which  may  lead  to  the 
creation  of  voids  in  the  subgrade  soil  along  the  joint 
and  cause  earlier  pavement  failure.  To  reduce  pavement 
deflections  along  the  joint,  stabilizing  the  subgrade 
soil  was  found  more  beneficial  than  using  the  butted 
joints. 

£.  Compared  with  a  regular  jointed  concrete  pavement,  the 
stresses  in  a  CRCP  are  reduced  if  the  pavement  is  kept 
in  good  condition.  However,  when  the  number  of  cracks 
is  increased  and  consequently  the  crack  spacing  is  re¬ 
duced,  not  only  the  stresses  start  to  increase,  but  also 
the  location  of  the  maximum  stress  under  the  edge  load 
shifts  from  the  pavement  edge  to  the  interior  part  of 
the  pavement.  As  the  load  transfer  capability  across 
the  crack  is  reduced  to  zero  due  to  the  increased  load 
applications,  the  stresses  increase  drastically  in  the 
center  load  case  reaching  that  under  the  edge  load. 

This  drastic  increase  in  stresses  at  locations  away  from 
the  pavement  edge  induced  by  either  the  edge  load  or  the 
center  load  can  be  attributed  to  the  punchout  failures 
frequently  observed  in  a  highway  CRCP  when  the  number  of 
cracks  is  increased  and  when  the  crack  openings  become 
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visible.  The  stresses  in  the  edge  load  case  can  be 
doubled  if  the  pavement  has  lost  the  subgrade  support 
along  the  pavement  edge  due  to  pumping,  temperature 
warping,  or  other  reasons. 

h.  In  a  CRCP,  stresses  under  a  full  truck  load,  i.e.,  both 
left  and  right  twin-tandem  loads,  are  less  than  those 
under  a  single  twin-tandem  load  placed  either  at  the 
pavement  edge  or  pavement  center.  The  edge  load  causes 
the  pavement  to  have  negative  stresses,  and  the  center 
load  causes  the  pavement  to  have  positive  stresses.  The 
combined  effect  of  edge  and  center  loads  reduces  the 
pavement  stresses.  The  maximum  stress  occurs  at  a 
section  about  3  to  5  ft  away  from  the  pavement  edge. 

Recommendations 


230.  The  computer  program  WESLIQID,  developed  for  rigid  pavements 
on  liquid  foundations,  is  versatile  since  it  has  many  options  to  deal 
with  problems  of  different  natures.  The  program  is  economical  to  oper¬ 
ate  and  requires  only  reasonable  computer  core  space.  It  is  recommended 
therefore,  that  the  WESLIQID  program  be  used  for  routine  pavement  design 
analysis,  and  research  purposes.  Because  the  computer  program  WESLAYER, 
developed  for  rigid  pavements  on  layered  elastic  solids,  is  limited  only 
to  two-slab  systems  and  requires  more  computer  core  space  and  computer 
time  to  operate,  it  is  recommended,  therefore,  that  the  program  be  used 
for  research  and  analysis  purposes. 
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Table  1 

Range  of  Modulus  of  Dowel  Reactions  from 
Various  Sources  (Finney  1956) 


Range  of  Modulus  of 
Dowel  Reaction 


10^  pci 

Source 

Remarks 

0.3  to  1.5 

Grinter 

Estimation 

Max .  2.5 

Friberg 

Tests  on  embedded 
dowels ,  1938 

0.71  to  1.17 

MS  HD* 

Load-deflection 
test,  1947 

0.78  to  5.89 

MS  HD* 

Tests  on  embedded 
dowels ,  1947 

0.89  to  8.3 

Marcus 

Dowels  with  uniform 
bearing  pressure 

2.45 

Loe 

Load-deflect ion 
tests,  1952 

0.9  to  8.6 

MS  HD* 

Tests  on  embedded 
dowels,  1954 

# 


Michigan  State  Highway  Department. 


Table  2 

Comparison  of  Stresses  Computed  from  the  Westergaard  Solution 
and  the  Finite  Element  Method  for  Three  Different  Loading 
Conditions  at  Different  Suugrade  Soil  Conditions 
P  =  10,000  lb,  E  =  3,000,000  psi,  u  =  0.15 


Maximum  Stress 

in  Slab  with 

Loaded  Radius 

a  =  0  ,  psi 

Modulus 

of  Westergaard': 

3  Finite 

Element  1 

Method 

Thickness 

Subgrade 

Solution* 

WESLIQID 

of  Slab, 

h  Reaction, 

k  Corner 

Interior 

Edge  Corner 

Interior 

Edge 

in. 

pci 

Load 

Load 

Load  Load 

Load 

Load 

50 

300 

181 

287  272.6 

200.9 

301.1 

10 

100 

300 

172 

270  262.0 

191.1 

295.0 

200 

300 

162 

253  247.7 

181.3 

271.4 

*  Computed  values  are 

available 

in  Tables  II,  III,  and  IV 

of 

Westergaard  (1925). 

Table  3 

Distribution  of  Stresses  Around  the  Concentrated  Corner  Load  in 

a 

Concrete  Pavement, 

P  =  10,000 

lb,  h  *  10 

in. ,  E  =  3,000,000  psi 

1. 

u  =  0.15 

,  k  =  100  pci 

Node 

0 

X 

0 

y 

T 

0 

max 

0  . 
min 

T 

max 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

0.0 

-154.8 

0.0 

-154.8 

0.0 

77.4 

3 

0.0 

-174.3 

0.0 

-174.3 

0.0 

87.1 

10 

-154.8 

0.0 

0.0 

-154.8 

0.0 

77.3 

11 

-97.1 

-97-1 

164.9 

262.0* 

67.8 

164.9 

12 

69.8 

-134.3 

123.2 

-229.4 

25.3 

127.4 

19 

-174.3 

0.0 

0.0 

-174.3 

0.0 

87.1 

20 

-134.3 

-69.8 

123.2 

-229.4 

25.3 

127.4 

21 

-96.9 

-96.9 

97.8 

-194.7 

91.5 

97.8 

*  Maximum  stress. 


Table  4 

Comparison  of  Stresses  and  Deflections  Between  th 


Finite  Element  and  the  Discrete  Element  Methods 


Maximum  Principal  Stresses,  psi 


Finite  Discrete 
Element  Element  Difference 


t  691. 


1101.8 

77U.3 

0.0 

1U68.7 

1583.2 

11+27.8 


1091.2 


-989.5 
1021.2 
10i*6. 6 


753.i 
675.' 
1*56.' 
-55. 
1550.0 


1090.0 


1270.0 

1100.0 


1510.0 
0. 
1500.0 


1310.0 


1370.0 


65. 


-1000.0 

-1030.0 

-1060.0 


— 

0.398  0 

— 

0.1*37  0 

— 

0.1*73  0 

— 

0.505  0 

— 

0.535  0 

5.25 

0.358  0 

2.87 

0.39**  0 

2.21 

0.1*27  0 

0.0 

0.1*58  0 

— 

0.1*86  0 

0.53 

0.278  0 

1.82 

0.307  0 

1.86 

0.335  0 

0.1+1* 

0.361  0 

— 

0.387  0 

1.00 

0.172  C 

0.33 

0.192  0 

1.02 

0.211  0 

1.23 

0.230  0 

— 

0.21+9  0 

1.60 

0.093  0 

0.96 

0.106  c 

0.76 

0.119  c 

0.90 

0.131  0 

n  i  n 

1.37 

0.01+2  0 

1.05 

0.0l*9  0 

0.85 

0.057  0 

1.26 

0.065  0 

— 

0.072  0 

Vertical  Deflections,  in. 


Finite  Discrete 
Element  Element  Difference 


3.1b 

3.10 

2.87 

2.88 
2.73 
3. 21* 
2.96 
2.95 
2.76 
2.80 
2.80 
2.85 
2.62 
2.70 

2.52 
2.82 
2. 51* 
2.76 
2.5l* 
2.35 
3.13 
2.75 
2.U6 
2.2l* 
2.05 
2.33 
3.92 
3.39 

1.52 


*  Nodes  along  the  edge. 

**  Node  at  the  corner. 

t  Positive  stress  indicates  tension  at  the  bottom  of  the  slab 
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Strain  measurements  are  indiated  in  microinches. 

Total  strain  is  the  sum  of  the  strain  measured  at  the  loaded  and  unloaded  sides  of  the  joint 


Table  7 

Summary  of  Average  Load  Transfer  for  Doweled  Joints 
(Ohio  River  Division  Laboratories  1939) 


Average  Measured  Strain 
_ *  10  ^  in. /in. _ 


Location 

Method 
of  Dowel 
Installation 

Loaded 
Side  of 
Joint 

Unloaded 
Side  of 
Joint 

Total 

Strain 

Load 

Transfer 

% 

Lockbourne  AFB, 
Ohio 

Cast  in  place 

127 

3*+ 

l6l 

21.1 

Lincoln  AFB, 
Nebraska 

Cast  in  place 

55 

32 

87 

36.8 

Hunter  AFB, 
Georgia 

Remove  and 
replace 

70 

26 

96 

27.1 

McCoy  AFB, 
Florida 

Remove  and 
replace 

6o 

18 

78 

23.1 

ORDL  Test 

Track  "A" 

Remove  and 
replace 

82 

32 

111* 

28.1 

Ellsworth  AFB, 
South  Dakota 

Dummy 

half-dowel 

6U 

U3 

107 

L0.2 

Beale  AFB, 
California 

Oversize  dummy 
half-dowel 

67 

32 

99 

32.3 

March  AFB, 
California 

Split  dowel 

iou 

50 

15b 

32.5 

Dow  AFB,  Maine 

Split  dowel 

5b 

7 

6l 
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is  greater  than  0.01  in.,  a  K  value  of  1,500,000  psi  is  used. 
Computed  by  the  WESLIQID  finite  element  program. 


Table  9 
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Solutions  and  Experimental  Measurements 
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Table  15 

Cocparisons  cf  Stresses  and  Deflections  fcr  Different  Loading  Positions  Shown  in  Figure  H  , 
A  Twc-Slab  Pavement  System  with  Interior  Loads;  tfESLIQID  Program 


r*"’ 


The  stress  ox  is  expressed  in  psi.  Positive  stress  indicates  that  the  tensile  stress  is  at  the  bottom  of  the  slab  and  negative  stress  indicates  that 
the  tensile  stress  is  at  the  top  of  the  slab.  The  deflection  is  expressed  in  inches.  Positive  deflection  represents  upward  gap  and  negative  deflection 
represents  the  deformation  sinking  into  the  ground. 


Center  Loads,  WESLIQID  Program 
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of  the  pavement. 

The  locations  of  nodal  points  are  shown  in  Figure  27. 
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*  The  locations  of  nodal  points  are  shown  in  Figure  27. 

'*  Maximum  stressses.  The  maximum  free  edge  stress  is  1151.5  psi  (o 
t  The  coordinates  of  x-  and  y-axis  are  shown  in  Figure  27. 


Letters  refer  to  loading  positions  shovn  in  Figure  28. 
Maximum  stresses  or  deflections. 


-* 

CVJ 

OJ 

o 


o 

o 


co 

on 

o 

o 


C  P-« 
■H  O 
.  I  « 
No  o 


C!  Pk 
•H  O 
,  I  « 
CO  o 


31 


o 

CO 

cm 


CO 

i—4 

o 


CO 

co 


CO 

o 


On 

O 

l— 


-if 

o 


CO 

00 

00 


LA 

CO 

o 


*M 

4- 

•  cd 

OJ 

-if 

NO 

On 

d  H 

-=r 

t— 

rH 

VO 

•h 

ro 

H 

rH 

O 

l  M 

• 

CO  0) 

o 

O 

O 

O 

« 

ro 

O 

NO 

LA 

C  Pk 

• 

• 

• 

•H  O 

ON 

o 

CO 

00 

1  « 

co 

-S' 

O 

LA 

NO  O 

CM 

O 

ON 

1 

rH 

rH 

. 

o 

VO 

rH 

-S 

C  P* 

• 

• 

• 

• 

•H  O 

NO 

LA 

rH 

On 

1  K 

00 

00 

IA 

rH 

co  o 

LA 

1 

t— 

NO 

NO 

+~ 

•  3 

o 

-S 

O 

On 

C  H 

• 

• 

• 

• 

•H  d 

On 

t- 

-if 

-S 

0&  § 

O 

o 

LA 

rH 

t— 

CO 

NO 

r — 

K 

rH 

H- 

• 

<M 

NO 

ON 

NO 

C  Ck 

• 

• 

• 

• 

•H  o 

00 

LA 

LA 

rH 

vAg 

ON 

CO 

ON 

LA 

ro 

00 

1 

On 

On 

On 

f- 

LA 

00 

d  Pm 

• 

• 

• 

• 

•H  O 

CO 

-S 

NO 

CO 

cAg 

i — 1 

00 

CM 

CM 

1 

-S 

NO 

rH 

VO 

u 

c- 

0- 

-=t 

NO 

•  CO 

d  H 

VO 

00 

OO 

t— 

•h  d 

rH 

oo 

LA 

ON 

cA  8 

CM 

CM 

1 

VO 

NO 

K 

• 

CO 

o 

ON 

-S 

d  Ck 

• 

• 

• 

• 

•H  O 

H 

o 

LA 

o 

1  K 

-S 

-if 

t— 

CO 

NO  O 

t— 

CM 

CM 

- f 

rH 

1 

. 

-S 

VO 

0- 

-S 

d  p< 

• 

• 

• 

•rH  O 

LA 

LA 

00 

-if 

1  03 

-S 

co 

VO 

C\J 

co  O 

-s 

f- 

rH 

-if 

o 

rH 

O 


o 

d 

0 

a 

p  _ 

4J 

P 

4J 

4J 

'd 

4) 

P  C0 

> 

aJ 

P 

>  ^ 

><  o 

P 

0 

0 

H 

3 

P  (0 

4J  rH 

4> 

rH 

0 

V 

0 

4)  O 

c 

P 

•H 

a 

rH 

43  rH 

u 

b 

P 

4) 

»d 

4) 

•H 

0> 

fd  a> 

4)  P 

4> 

P 

CO 

0> 

bO 

(V  M 

o  a 

O 

d 

o 

o 

'd 

o 'd 

(0  OJ 

a) 

4> 

pk 

(0 

4) 

C0  4) 

H  O 

rH 

O 

H 

rH 

04 

P< 

P4 

r> 

d  «. 

* 

r 

d 

P 

P 

■d  p 

'S 

P 

•H 

CO 

c 

co  a 

a)  d 

a 

d 

t* 

'd 

•rH 

*d  »h 

O  -H 

0 

•H 

oJ 

S 

o 

C0  o 

rH  O 

H 

O 

Q 

0 

•o 

o  ^ 

•o 

rH 

rH 

U 

V 

<V 

a 

<V  4 ) 

4/ 

V 

o; 

43 

4)  43 

P  J3 

p 

43 

M 

P 

5) 

d  p 

d 

P 

*3 

0) 

4> 

w 

w 

o 

O 

M 

U 

■H 

43 

P 

I 


MOV. 
c  as  O 

O  CO 

S.  ^  ^ 


CO  QJ  CO 
V. 

-  M  w  O 
d  o  c  » 
*h  a)  aJ  p 
I  U  u 
CO  O  P  4> 

—  -A  B  oT 

CM  O  O 

Be 

**  fH  P 
P  >  C  4) 
C  0)  43 
d,  o  P 
LA  O  V 

"  J  c 

O.  *H 


O  0) 
V  H 


V. 

O  44 
O  - 

43  *H  4>  .  . 
P  X  N  £ 
M  P  > 
C  I  'O 
0)  •  c  P 
H  G  «J  «J 

•P  rH 
CO 


P  I  V 
d  VO  <U 

<U  <M 

a  ca  “ 


<D 

as* 


~  d 
P  *H 


cd  xl  ^  *d 

43 
P 

s 


SO)  4) 

4;  o 


4= 

CO 


r*> 

O  *d 

c  a) 
o  o 


Ctf  ^  _ 

p«  a  o  h 

aJ  -h  a  m 
tJ  JH  y  4> 

4>  O  V.  *d  4) 

*  -  *  sts 

V.  P  V 

o  i  3 

d 


■SC 


0k 


APPENDIX  A:  EQUATIONS  FOR  STRESSES  AND  DEFLECTIONS 
UNDER  A  POINT  LOAD  AND  UNDER  A  CIRCULAR  LOAD 
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1.  The  stresses  and  vertical  deflection  in  a  linear  elastic 
medium  under  a  point  load  can  be  computed  from  the  following  equations 
(Harr  1966). 
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normal  stresses  in  the  r  ,  0  ,  and  z 

directions 

a  point  load 

see  Figure  Al 

shear  stress  in  the  rz  direction 
'vertical  deflection 
Poisson's  ratio 
elastic  modulus 


2.  Under  a  circular  load,  the  vertical  stress  and  deflec¬ 

tion  w  under  the  center  of  the  circular  load  can  be  computed  by  the 
following  two  equations  (Harr  1966) 
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where 

q  =  unit  pressure 
a  =  radius  of  the  loaded  area 
n  =  z/a 


APPENDIX  B :  NOTATION 
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[H] 


Hij 
i * j 


I 
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[K] 

[K] 
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p' 


Dimensions 

Area  of  the  dowel  bar 

Diagonal  matrix  representing  the  area  over  which  subgrade 
reaction  is  distributed 

Dowel  bar  diameter 

Initial  curling  of  a  weightless  and  unrestrained  slab  due 
to  a  temperature  differential  between  the  top  and  the  bottom 

Joint  spacing  or  one  half  of  the  slab  length 

Modulus  of  concrete  or  modulus  of  dowel 

Vertical  force 

Nodal  forces 

Composite  nodal  forces 

Shear  modulus  of  steel 

Thickness  of  the  concrete  slab 

Vertical  force  at  node  i  due  to  vertical  displacement  at 
node  J 

Stiffness  matrix  of  the  subgrade 

Vertical  force  at  node  i  due  to  displacement  at  node  j 
Node 

Moment  of  inertia  of  the  slab  section 

Modulus  of  subgrade  reaction 

Modulus  of  dowel  support 

Stiffness  matrix  of  the  concrete  slab 

Composite  stiffness  matrix  of  the  system 

Vertical  force  at  node  i  due  to  vertical  displacement 
at  node  j 

Radius  of  relative  stiffness 
Length  of  the  concrete  slab 
Moments 

Moment  about  the  x-axis 
Moment  about  the  y-axis 
Unit  pressure 

Vertical  shear  force  or  a  point  load 

Shear  force  on  dowel  to  affect  A' 
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B1 
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R 


A  composite  matrix  and  is  equal  to  [K]  +  [ H ] 

Radius  of  the  spherical  surface 
s  Precomposition  or  gap 

t  Pavement  thickness 

w  Deflection  in  the  z-direction 

a  Coefficient  of  thermal  expansion 

8  A  parameter 

{6}  Nodal  displacements,  each  consisting  of  a  vertical  deflec¬ 
tion  and  two  rotations 

{6*}  Subgrade  displacement 

A  Difference  in  deflection  across  the  joint 
Ac  Deformation  of  concrete  due  to  shear  force  on  the  dowel 

,  Input  parameter  specifying  the  deformation  of  concrete  below 

c  which  modulus  K-^  is  used  and  above  which  modulus  Kp  is 

used 

As  Shear  deformation  of  the  dowel  bar 

AT  Temperature  differential 

0x  Rotation  about  the  x-axis 

0y  Rotation  about  the  y-axis 

v  Poisson's  ratio  of  the  concrete 

a  „  Maximum  principal  stress 

max  ^ 

ax  Stress  in  the  x-direction 

oy  Stress  in  the  y-direction 

t  Shear  stress 

Shear  stress  in  the  x-y  direction 
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In  accordance  with  letter  from  DAr.M - RDC ,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Chou,  Yu  T. 

Structural  analysis  computer  programs  for  rigid 
multicomponent  pavement  structures  with  dis¬ 
continuities — WESLIQID  and  WESLAYER  :  Report  1  : 
Program  development  and  numerical,  presentations  /  by 
Yu  T.  Chou  (Geotechnical  Laboratory,  U.S.  Army 
Engineer  Waterways  Experiment  Station).  —  Vicksburg, 
Miss.  :  The  Station  ;  Springfield,  Va.  :  available 
from  NTIS ,  [ 1981 ] . 

I1 * * IV.*?  p.  in  various  pagings  :  ill.  ;  27  cm.  — 
(Technical  report  /U.S.  Army  Engineer  Waterways 
Experiment  Station  ;  GL-81-6,  Report  l) 

Cover  title. 

"May  1981." 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army 
under  RDT&E  Project  No.  UA762719AT1iO,  Work  Units  001 
and  003." 
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